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Abstract 

Over the past one hundred years Brownian motion theory has contributed substan- 
tially to our understanding of various microscopic phenomena. Originally proposed 
as a phenomenological paradigm for atomistic matter interactions, the theory has 
since evolved into a broad and vivid research area, with an ever increasing num- 
ber of applications in biology, chemistry, finance, and physics. The mathematical 
description of stochastic processes has led to new approaches in other fields, cul- 
minating in the path integral formulation of modern quantum theory. Stimulated 
by experimental progress in high energy physics and astrophysics, the unification of 
relativistic and stochastic concepts has re-attracted considerable interest during the 
past decade. Focusing on the framework of special relativity, we review here recent 
progress in the phenomenological description of relativistic diffusion processes. Af- 
ter a brief historical overview, we will summarize basic concepts from the Langevin 
theory of nonrelativistic Brownian motions and discuss relevant aspects of relativis- 
tic equilibrium thermostatistics. The introductory parts are followed by a detailed 
discussion of relativistic Langevin equations in phase space. We address the choice of 
time parameters, discretization rules, relativistic fluctuation-dissipation theorems, 
and Lorentz transformations of stochastic differential equations. The general theory 
is illustrated through analytical and numerical results for the diffusion of free rela- 
tivistic Brownian particles. Subsequently, we discuss how Langevin- type equations 
can be obtained as approximations to microscopic models. The final part of the arti- 
cle is dedicated to relativistic diffusion processes in Minkowski spacetime. Since the 
velocities of relativistic particles are bounded by the speed of light, nontrivial rela- 
tivistic Markov processes in spacetime do not exist; i.e., relativistic generalizations 
of the nonrelativistic diffusion equation and its Gaussian solutions must necessarily 
be non-Markovian. We compare different proposals that were made in the liter- 
ature and discuss their respective benefits and drawbacks. The review concludes 
with a summary of open questions, which may serve as a starting point for future 
investigations and extensions of the theory. 

Key words: Brownian motion, special relativity, stochastic processes, relativistic 
Langevin equations, Fokker-Planck equations, diffusion processes, relativistic 



Preprint submitted to Elsevier 



26 January 2009 



thermodynamics 

PACS: 02.50.Ey. 05.40.-8. 05.40.Jc. 47.75.+f 



Email addresses: jorn.dunkelOphysics.ox.ac.uk (Jorn Dunkel), 
peter.hanggi@physik.uni-augsburg.de (Peter Hanggi). 

URLs: http : //www-thphys . physics . ox . ac . uk/people/ JornDunkel/ (Jorn 
Dunkel), http://www.physik.uni-augsburg.de/theol/hanggi/ (Peter Hanggi). 



2 



Contents 

1 Introduction 

1.1 Historical background 

1.2 Relativistic diffusion processes: problems and general strateg 

1.3 Structure of the review 

2 Nonrelativistic Brownian motion 

2.1 Langevin and Fokker-Planck equations 

2.2 Microscopic models 

2.3 Remarks and generalizations 

3 Relativistic equilibrium thermostatistics 

3.1 Preliminaries 

3.2 Stationary systems with confinement 

4 Relativistic Brownian motion processes in phase space 

4.1 Relativistic Langevin and Fokker-Planck equations 

4.2 One-dimensional examples and mean square displacement 

4.3 Proper-time reparameterization 

4.4 Moving observers 

4.5 Relativistic binary collision model 

5 Non-Markovian diffusion processes in Minkowski spacetime 

5.1 Reminder: nonrelativistic diffusion equation 

5.2 Telegraph equation 

5.3 Relativistic diffusion propagator 

6 Outlook 
Symbols 

A Stochastic integrals and calculus 



A.l Ito integral [l06 

A. 2 Stratonovich-Fisk integral IIO8 

A. 3 Backward Ito integral IllC 

A. 4 Comparison of stochastic integrals Ill2 

A. 5 Numerical integration Ill3 

B Surface integrals in Minkowski spacetime Ill5 

C Relativistic thermodynamics ll2C 

C.l Reminder: nonrelativistic thermodynamics |12C 

C.2 Relativistic case 121 



4 



1 Introduction 




In his annus mirabilis 1905 Albert Einstein published four manuscripts [H, 
that would forever change the world of physics. Two of those pa- 
[sl laid the foundations for the special theory of relativity, while an- 
other one [3] solved the longstanding problem of classical (nonrelativistic) 
Brownian motionj^ Barring gravitational effects j^, 0], special relativity has 
proven to be the correct framework for describing physical processes on all 
terrestrial scales 0, @] • Accordingly, during the past century extensive efforts 
have been made to adapt established nonrelativistic theories such as, e.g., 
thermodynamics, quantum mechanics or field theories to the requirements 
of special relativity. Following this tradition, the present review focuses on 
recent progress in the theory of special relativistic Brownian motion and dif- 



26l, HQ, 0, [SBHlii, IS 



fusion processes M [nL MJ, [m [JJ, ll5|, ll6|, ll7|, ll8|, ll9|, l20|, M, l22, l23|, M, 125 



1.1 Historical background 



Historically, the term 'Brownian motion' refers to the irregular dynamics 
exhibited by a test particle (e.g., dust or pollen) in a liquid environment. 
This phenomenon, already mentioned in 1784 by the Dutch physician Ingen- 
Housz 



35 



361 ] ■ was first analyzed in detail by the Scottish botanist Robert 



Brown 37] in 1827. About eighty years later, Sutherland 38], Einstein and 
von Smoluchowski 39|] were able to theoretically explain these observations. 
They proposed that Brownian motion is caused by quasi-random, microscopic 
interactions with molecules forming the liquid. In 1909 their theory was con- 



firmed experimentally by Perrin [40|, providing additional evidence for the 
atomistic structure of matter. During the first half of the 20th century the 
probabilistic description of Brownian motion processes was further elaborated 
in seminal papers by Langevin 4ll . |42| |. Fokker |43ll, Planck 14J1, Klein |45 



Uhlenbeck and Ornstein 46 and Kramers 



13 IIP 



In parallel with these early theoretical studies in the field ofphysics, out 



standin; 
roff 



mathematicians like Bachelier 
Feller 



54, 55, 56 



57[, and Levy 



Wiener [5l|, l52|, |53||, Kolmogo 



58l . l59l | developed a rigorous basis 



for the theory of Brownian motions and stochastic processes. Between 1944 



and 1968 their groundbreaking work was com plem ented by Ito [60|, l6l|], Gih- 
man 62, 63, 64], Fisk 65, 66l] and Stratonovich 67, 68, 69], who introduced and 
characterized different types of stochastic integrals or, equivalently, stochastic 



^ Einstein's first paper [l| provided the explanation for the photoelectric effect. 
^ Excellent early reviews are given by Chandrasekhar 48|], and Wang and Uhlen- 
beck \E 
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differential equations (SDEs). SDEs present a very efficient tool for modeling 
random processes and their analysis has attracted an ever-growing interest 

over the past decades 0, H Q H H H 17610 

Nowadays, the modern 

theory of stochastic processes goes far beyond the original problem considered 
by Einstein and his contemporaries, and the applications coyer a wide ran ge o f 
different areas includirig physics [73, [tS, 79, 80; 
economy and finance [8 



81|, 182|, 183|, l8j , biology [85|, l86 



A central topic of this review concerns the question how SDE-based Brownian 
motion models can be generalized within the framework of special relativity. 
In the physics literature 



75 



tions 



41, 42 



SDEs are often referred to as Langevin equa- 
and we shall use both terms synonymously here. From a math- 
ematical perspective, SDEs [tS] determine well-defined models of stochastic 
processes; from the physicist's point of view, their usefulness for the descrip- 
tion of a real system is a priori an open issue. Therefore, the derivation of 
nonrelativistic Langevin equations from microscopic models has attracted con- 
siderable interest over the past sixty years 2l|, |90|, |9l|, l92|, |93|, [OJ, l95|, l96| . Efforts 
in this direction not only helped to clarify the applicability of SDEs to physi- 
cal problems but le d, among others, also to the coiicept of quantuna Browni an 



motion 0, H H S, ffl, lioil, [loa, E, EqI, |Toi, Em, [To3, [TosI, [Toi, [noj 



If one aims at generalizing the classical Brownian motion concepts to special 
relativity, then several elements from relativistic equilibrium thermodynamics 
and relativistic statistical mechanics play an important role. More precisely, 
thermostatistic principles govern the stationary behavior of Brownian parti- 
cles and, thus, impose constraints on the structure of relativistic Langevin 
equations. The first pape rs on relativis tic t hermodynamics were published in 



1907 by Planck [l2ol . 



12ll | and Einstein [122|. A main objective of theirs was to 



identify the Lorentz transformation laws of thermodynamic variables (temper- 
ature, pressure, etc.)0 In 1963 the results of Einstein and Planck became ques- 



tioned by Ott |127) . whose work initiated an intense debate about the correct 



rela t ivist i c transformat i on behayior of thermodynamic quant i ties 



130 . 
146, 



131 



147 



132 



148 



133 



149 



135 



151 



136 



152 



137 



153 



138 



154 



139 



155 



140 



156 



141 



157, 



142 



158 




^ The history of the mathematical literature on Brownian motions and stochastic 
processes is discussed extensively in Section 2.11 of Ref. [74,]; see also Sections 2-4 
in Nelson [71 1. 

^ The vast literature on classical Brownian motion processes and their vari- 
ous app lications in non relat i vistic physics is discussed in several survey arti- 
cles (m ItqI. liol. [ill, lii. lii. [ml . fni . fm . Tni . fni l . Nonrelativistic generalizations of 
the st andard th eory as, e.g., anomalous diffusion processes have been summarized 
in 83|, IllO . Ill7l | , whi l e review articl es on nonrelativistic quantum Brownian motion 
can be found in [iOtI. lOsl . E^ . flli | . 

^ See also Pauli 123 |. Eddington 124], Tolman 125 ] and van Dantzig 1261 ] for 



early discussions of this problem. 
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Tei . [Tei . fiei . [Tesl . fieel . [TgtI . [igsI . [T69l .[T7ol . [ml . fiTi . (TtI [iTi, [TtsI, [Tt^ 



However, as clarified by van Kampen [150[ | and Yuen [178j, the controversy sur 



rounding relativistic thermodynamics can be resolved by realizing that ther- 
modynamic quantities can be defined in different, equally consistent ways. 

While some authors considered relativistic thermodynamics as a purely macro- 
scopic theory, others tried to adopt a more fundamental approach by focussing 
on relativistic equilibrium statistical mechanics. Early pioneering w ork i n the 



latter direction was provided by Planck's students von Mosengeil |184 | and 



von Laue 185 



186 , who derived in 1911 the 



^ and his collaborator Jiittner ^ 

relativistic generalization of Maxwell's velocity distribution 187| F1 Research 
on relativistic equilibriu m thermostati s tics experienced its mos t inte nse phase 



between 1960 and 1970 137, 138, 146 



tics experiencea its mos t intense pnast 
153,Ji3,j20i, ffl, ffl, H, I2q3 [208 



209I . [210I . [ml . [2li . iSl . mi UM UM- An excellent exposit 



tion 



on the conceptual foundatio ns and difficul ties of relativistic statistical me- 



221 



chanics was given by Hakim 
field has conti nued to attrac t interest 

223 . liisl . [22^ . 



201 



237 




22^, 



222 iJn 1967. Over the past 



m ligi. ITggl. [i97l. 



240, 241, 242 



2271 . 12281 . 12291 . 12301 . 12311 . 12321 . 1233 . 



years the 



234 



1991. 



2351. 



23^, 



The recurring debate on relativistic thermostatistics can be traced back to 
the difficulty of treating many-particle interactions in a relativistically consis- 
tent manner!^ In nonrelativistic physics interactions may propagate at infi- 



183, 



181, 18 



The pre-1970 literature o n thi s disputed issue has been reviewed by Yuen [ITi 
and Ter Haar and Wegland 179l ]. see also Israel 177l |: more recent surveys can be 



found in 

A gen eral introduction to relativistic gases is given in Synge's textbook clas- 
sic [188*]. Relativis tic g eneralizations o f equ ipartition and virial theorems 189] are 



discussed by Pauli 123 ] and Einbinder 19d| : rn ore recent investigations of these two 
specific topics can be found in, e.g., Refs. [l9l|, SEgj ]. Equilibri um d istributions 
for ideal relativistic quantum g a ses were also deriv ed by Jiittner 1941 ] in 1928. A 

Eai, Egg [l97l. Iigi . Iligl . liool ] have raised doubts about the 



few recent papers [22 

correctness of Jiittner 's results |l86l . Il94il. but relativistic molecular dynamics sim 



20ll.l20i]. cf. discussion in Sections 13.2.21 and 



ulations confirm Jiittner's prediction 
Ol below. 

^ Seminal contributions to the theory of relativist ic niany-p article i nter actions were 



provided by Fokker [243i] . Wheeler and Feynman 
Goldberg 247, 248t . and Van Dam and Wigner 



244 



250 



iiil . Pryce (m^], Havas and 



249l. l250t] . Oyer t he past decades 
25ll. I252I. lisij . These forbid. 



several no-interaction theorems were proven 
within their respective qualifications, certain types of interaction models within the 
framework of special relativity. The mathematical structure of relat ivist ic many- 
particle interactions was analyzed in detail by Arens and Babbitt 254]: v ariou s 
semi-relativistic approximations have been discussed, e.g., in Refs. i255l. bsdlisTj. 
Kerner 2581 ] has edited a reprint collection covering large parts of the pre-1972 
literature on relativistic action-at-a-distance models; more recent contributions in- 



clude ^5^, lieS, lisH, liei, liei, liei, H, lie^, H- An alternative 
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nite speed, i.e., they can be modeled via instantaneous interaction potentials 
which enter additively in the Hamilton function; from that point on, non- 



relativistic statistical mechanics emerges without much difficulty [2871 . |288 



Unfortunately, the situation becomes significantly more complicated in the 
relativistic case: Due to their finite propagation speed, relativistic interactions 
should be modeled by means of fields that can exchange energy with the par- 
ticles These fields add an infinite number of degrees of freedom to the 
particle system. Eliminating the field variables from the dynamical equations 
may be possible in some cases but this procedure typically leads to retar- 
datio n eff e cts, i .e., t he particles' e quations of motions become non-local in 



220l . [iill . H, O, H, [iij. Thus, in special relativity it is usually 



time 

very difficult or even impossible to develop a consistent field-free Hamilton 
formalism of interacting particles. 

In spite of the difficulties impeding a rigorous treatment of classical relativis- 
tic many-particle systems, considerable progress was made during the second 
half of the past century in construc t ing an approximat e relativis t ic k i netic 



300l.[301 



302, 303 



304 



305 



306 



theo rv Il67l. l225l. l255l. l289l. 1291 l29ll. l292l. l293l. l294 l295l . l29d l297l . l298l . l299 



307 



308 



309 based on relativistic Boltz- 



mann equations for the one-particle phase space probability density functions 
(PDFs)l£j From such a kinetic theory it is only a relatively small step to formu- 
lating a theory of relativistic Brownian motion processes in terms of Fokker- 
Planck e quations a nd Langevin equations. While the relativistic Boltzmann 
equation [3 111 . |312| ] is a nonlinear partial integro-differential equation for the 
PDF, Fokker-Planck equations are linear partial differential equations and. 



therefore, can be more easily solved or analyzed [73 



The present article focuses primarily on relativistic stochastic processes that 
are characterized by linear evolution equations for their respective one-particle 
(transition) PDFs. The corresponding phenomenological theory of relativis- 
tic Brownian motion and diffusion processes has experienced considerable 
progress during the past decade, with applications in various areas of high- 



ener gy physics liii . mdi . [3T7I . liisl . I3T9I . [320I . I321I . |322| | and astrophysics |323l . 



eiiergy_pi 




From a general perspective, relativistic stochastic pro- 



intensely studied method f or describing relatiyistic interact i ons is based on the so- 
called constraint formalism fl96l. &i. liofll . [27ll . 127^ . Itt^ . liri . [27rI . [27fil . liT?! . [27^. 
279l . [280I . [281I . [isi . lisi . [284, 28^ The foundations of this approach were worked 
out by Dirac 269f |. who aimed at cons t ruct i ng a relat ivistic quantum theory for 
interacting many- particle systems ^5^, [25l|, ES, lis^ . For a detailed discussion 
of relativistic many-particle theory, we re fer to the insightful c onsiderations in the 
original papers of Van Dam and Wigner and Hakim [22l|, [22^.[224| 

as well as to the recent review by Hakim and Sivak 286l |. 

^ Comprehensive introductio ns t o relativistic Bolt zman n equations can be found 
in th e te xtbooks by Stewart [310(], de Grootet al. |31ll ]. and Cercignani and Kre- 



mer 



3I4], or also in the reviews by Ehlers [3131 1 and Andreasson [314] . 



8 



cesses provide a useful approach whenever one has to model the quasi-random 
behavior of relativistic particles in a complex environment. Therefore, it may 
be expected that relativistic Brownian motion and diffusion concepts will play 
an increasingly important role in future i nvest i gations of , e.g. , thermalization 
and relaxation proces ses i n as t rophysics | 323l. l324j . l325l . |326| or high-energy 



collision experiments [315|, l316l . l318l . l319l . l328l . |329 



1.2 Relativistic diffusion processes: problems and general strategies 



According to our knowledge, the first detailed mathematical studies on rela- 
tivistic di ffusio n processes were performed independently by Lopuszahski 330 



Rudberg |33l| . and Schay |332l | between 1953 and 1961. In the 1960s and 70s 
their pion eering work was fu rther elaborated by Dudley who published a series 
of papers 3331 . l334j . l335l . |336[ | that aimed at providing an axiomatic approach to 
Lorentz invariant Markov processes 74 1 i n phase space. Inde pendently, a sim- 
ilar program was pursued by Hakim 220, 221 . 222 . 337, 338 1. whose insightful 
analysis h elped to elucidate the conceptual subt leties of relativistic stochastic 

) and Hakim (Proposition 2 



processes 
in 



j338 |. Dudley (Theorem 11.3 in 333 



338| ) proved the non-existence of nontrivia^ 



10 



Lorentz invariant Markov 



processes in Minkowski spacetime, as already suggested by Lopuszahski [330 



This fundamental result implies that it is difficult to find acceptable r elativistic 
generalizations of the well-known nonrelativistic diffusion equation 2871 . 1339 



|. = BVV (1) 

where T) > is the diffusion constant and g{t, x) > the PDF for the particle 
positions as G M'^ at time t. In order to circumvent this 'no-go' theorem for 
relativistic Markov processes in spacetime, one usually adopts either of the 
following two strategies 

• One consider s non-Markoy ian diffusion processes X{t) in Minkowski space- 
time 0, 0, H, H, [iej . 



'nontrivial' if a typical path has a non- 



A diffusion process is considered as 
constant, non- vanishing velocity. 

The mathematical interest in r elativistic diffusion processes increased in the 1980s 
and 1990s, when several auth ors l34nl. |34ll . 'Hii. IHi^ . [34^ . [3491 . [Hiol . 
35ll . [35^ . B, H H H I357I. Issalssa . .360 . .36L ,36Z] considered the possibility 
of extending Nelson's stochastic quantization approach [363] to the framework of 
special relativity. These studies, although interesting from a mathematical point of 
view, app ear t o have relatively little physical relevance because Nelson's stochastic 
dynamics 363l | fails to r epro duce the correct quantum correlation functions even in 



the nonrelativistic case [364]. 
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One constructs relativistically acceptable Markov processes in phase space 
by considering not only the position coordinate X(t) of the diffusing par ticle 



but also its momentum c oordinate P(t) fill - 12, [isl, 'iS, 

M H H H M, H H H, H, B, S Ei,l£i, 



17, 



337 



18, 19 



338 



1.2.1 Non-Markovian diffusion models in Minkowski spacetime 



A commonly cons i dered 're lativistic' generalization of Eq. ([T]) is the telegraph 



equation [lO, 365 . 366 , 367 
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(2) 



where > denotes a finite relaxation time scale. Unlike the classical diffu- 
sion equation ([T]), which is recovered for r„ = 0, the telegraph equation ([2]) con- 
tains a second order time- derivative and, therefore, describes a non-Markovian 
process. While the classical diffusion equation ([1]) permits superluminal propa- 
gation speeds, the diffusion fronts described by Eq. ([2]) travel at finite absolute 
velocity v = {Ti/t^Y/'^; cf. the discussion in Section below. 



Historically, Eq. ([2]) was first obtained by Thomson (Lord Kelvin) in 1854 365 



while studying the signal transduction for the first transatlantic cable. The 
first probabili stic d erivation of Eq. ([2]) for the one-dimensional case was given 
by Goldstein |366l | in 1950. His approach was based o n a s o-ca lled persistent 
random walk model originally introduced by Fiirth 3681 . |369| in 1917 as a 



paradigm f or d iffusive motion in biological systems and later also considered 
by Taylor 370[ | in an attempt to treat turbulent diffusion[2H In contrast to 
standard non-directed random walk models, which lead to the classical diffu- 
sion equation ([T]) when performing an appropriate continuum limit [t^], the 
random jumps of a persistent walk take into account the history of a path 
by assigning a larger probab ility to th ose jumps that point in the direction of 
the motion before the jump 3661 . l367l |. Persistent random walk models ca n be 
used to descri be the transra ission of light i n rnultip le scattering media 372 
such as foams |373l . l374j . l375l | and thin slabs |37ll . |376| . Similarly, the telegraph 
equation ([2]) has been applied in various areas of physics over the past decades 
e.g., to model the propagation of electric signals and heat waves 



14 



Moreover, an interesting connection between the free particle Dirac equa- 



Masoliver and Weiss |lO|| discuss several possibilities of deriving Eq. 
different underl ying models. 



(121) from 



2i.l371| 



See also Kac 3671] and Boguna et al. 

A detailed review of the pre-1990 research on heat waves was provided by Joseph 
and Preziosi |377l. 13781] . wh ile more r ecent discussi ons and applications of Eq. 1^ 
can be found hTflJlii, 0, H H H H IsSlj . 



10 



tion 



382l | and the telegraph equation ([2]) was pointed out by Gaveau et al. |383 



in 1984: The solutions of both equations may be linked by means of an ana- 
lytic continuation quite similar to the relation between the classical diffusion 
equation ([T]) and the free-particle Schrodinger equation in the nonrelativistic 



case 



15 



On the other hand, the telegraph equation ([2]) is not the only pos- 
sible generalization of Eq. ([1]) and a rather critical discussion of Eq. ((21) in 
the context of relativistic heat transport was given by van Kampen [1531] in 
1970. Section [5] below takes a closer look at the properties of Eq. ([2]) and also 
addresses potential alternatives 23, 39d . 



1.2.2 Relativistic Markov processes in phase space 



As an alternative to non-Markovian diffusion models in spacetime, one can 
consider relativistic Markov processes in p hase spac e [11 



17 



18Ul9|j20L 121|, 122|, 12J, [26, 



31, 32, 33, 34, 220, 221,^ 




336l . l337l . |338| . Typical examples are processes described by Fokker-Planck 



equations (FPEs) or Langevin equations 



11 



12, 13, 14 



)ed by tok 



20|,l2l|,l24 



25|, m, iMl, 132, l33|, l3J, [3911, [392, l393|, l39J, l395| . Similar to the relativistic Boltz- 
mann equation, relativistic FPEs in phase space can be used to describe non- 
equilibrium and relaxation phenomena in relativistic many-particle systems. 
FPEs can be derived from Lange vin eq uations, as approximations to more gen- 



eral linear master equations |332( | or by approximating the collision inte- 



grals in nonlinear Boltzmann equation by differ entia l expressions that contains 
effective friction and diffusion coefficients 16|, l299l |. In particular, the latter 
method was successfully applied i n different areas of physics over the past 
decades, incl uding plas ma physics |299|, H, H, H,|399|,j4o3,j40ll, fioil, C 



1404 14051 . 14061 . 14071 . |408| hi g h-ene rgy physics Imlili . 13171 . IsTl IsTft . ,328, . .329lj , 
and astrophysics [sii, H, H, Efl [411 the 1980s and 90s this 

approach was further elaborated 13961. 1397l.l398l.l399l. liool . lioil. lioi. lioi. 



404 



418 



and several numerical methods for solving FPEs were developed |399l . l417 

Recent applications include the modeling of diffusion and thermaliza- 
tion processes in qua rk-gluon plasmas, as produce d in relativistic heavy ion 



collision experiments 315 . 316, 318, 319 . 328, 329l|, as wel l as the description 



of c omplex high-energy processes in astrophysics 
iH. 



323, 324, 325, 409, 410, 411 



15 Pqj. fupi^her reading about path integral r epresentations of the Dirac propagator 
we refer to the papers of Ichinose 384 . 385], Jacobson and Schulman 38d |. Barut 



and Duru [383] , and Gaveau and Schulman 3881; see also footnote 7 in Gaveau et 



al. 3831 ] and pp. 34-36 in Feynman and Hibbs [389i] . 

Relativistic Fokker-Planck-type equations also played a role in the debate about 
whether or not the black body ra diation spectrum is coin patible with Jiittner's 
relativistic equilibrium distribution 19ll . |413| . 1414 . |415| . 14161 ] . 

A main reason for the lively interest in relativistic FPEs at that time was the 
prospect of building plasma fusion reactors. 
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A complementary approach towards relativistic stochastic processes in phase 



space starts fro m Lange v in e q uati o ns [111 . Il2l . Il3l . Il4l . Il5l . llTl . Il8l . l20l . l2ll . |24 



25l . l26l . l3ll . l32l . l39ll . l392l . l393l . 1394 |395| . Stochastic differential equations of 



the Langevin type yield explicit sample trajectories for the stochastic motion 
of a relativistic Brownian particle. Relativistic Langevin eq uations may ei- 
ther be postulated as phenomenological model equations ll|, [iTI] or obtained 



from more precise microscopic models by imposing a sequence of approxima- 
Compared with the nonrelativistic case, the latter task becomes con- 



21 



tions 

siderably more complicated due to the aforementioned conceptual and tech- 
nical difficulties in formulating consistent relativistic many-particle theories. 
The phenomenological Langevin approach to relativistic Brownian motion was 



initiated by Debbasch et al. [Uj, who in 1997 proposed a simple relativistic 



generalization of the classical Ornstein-Uhlenbeck process j46(], representing a 



special limit cases of a larger class of relativistic Langevin processes |15l. 1 191.126 



From a practical point of view, relativistic Langevin equations provide a useful 
tool for modeling the dynamics of relativistic particles in a random environ- 
ment, because these SDEs may be simulated by using well-established Monte- 
Carlo techniques that are numerically robust and efficient 76|, |89|, 14201 ] . Recent 



appli cations include the analysis of thermalization effects in quark-gl uon p las- 



mas 



318 , 319, 329, 421] and ultrarelativistic plasma beam collisions 3261 



1.3 Structure of the review 



The present article intends to provide a comprehensive introduction to the 
theory of relativistic Brownian motions with a particular emphasis on rela- 
tivistic Langevin equations. For this purpose, the subsequent parts are orga- 
nized as follows. Section [2] summarizes the Langevin theory of nonrelativistic 
Brownian motions in phase space. Section [3] discusses aspects of relativistic 
equilibrium thermostatistics as far as relevant for the subsequent discussion. 
Relativistic Langevin equations in phase space and their associated FPEs are 
considered in SectionUl SectionOis dedicated to relativistic diffusion processes 
in Minkowski spacetime; as outlined above, such processes must necessarily be 
non-Markovian. Our review concludes with a summary of open questions in 
Section [6], which may serve as a starting point for future investigations and ex- 
tensions of the theory. Since the topic of this article resides in the intersection 
between statistical and high-energy physics, we hope that the presentation is 
accessible for members of both communities. To keep the discussion as self- 
contained as possible, the article includes three appendices which summarize a 
few essentials about stochastic integrals, surface integrals in Minkowski space- 
time and relativistic thermodynamics. 
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2 Nonrelativistic Brownian motion 



This section summarizes basic definitions and mathematical tools as well as 
relevant results from the nonrelativistic theory of nonrelativistic Brownian mo- 
tions. In particular, the Langevin and Fokker-Planck equations considered in 
this part define the nonrelativistic limit case of the relativistic theory described 
Section m In Section [TT] we shall briefly recall the general mathematical struc- 
ture of Langevin and Fokker-Planck equations, the relevance of discretization 
rules and the choice of fluctuation-dissipation relations. Section 12.21 focusses 
on the question how stochastic differential equations (SDEs) can be motivated 
and/or derived from microscopic models. As typical examples, the well-known 



90 



9ll . |92| . l93l . 194 l95l . |96| and a recently proposed 



21 will be considered. In contrast to the oscillator 



harmonic oscillator model 
binary collision model 
model, the collision model can be generalized to the framework of special 
relativity, and its relativistic version will be discussed in Section 14.51 



2.1 Langevin and Fokker-Planck equations 



112 



The condensed discussion in this part is primaril y b ased on the papers o f Uh 
lenbeck and Ornstein [46] , Wang and Uhlenbeck j49| 



and Klimontovich 

For further reading about nonrelativistic stochastic processes and their numer- 
ous applications in physics a nd m athematics, we refer to the review articles 



of Chandrasekhar 
Georges 



48| . Fox 111| . Hang gi and Thomas [77], 



116| . Metzler and Klafter |117l |. H anggi a nd Marchesoni [81 



and Kroy [82 



Bouchaud and 
Frey 



7J, l75|, l76|, |422|. For conceptual 



or the textbooks references 
clarity, we restrict ourselves to the simplest case where motions are confined 
to one space dimension {d = 1). The generalization to higher space dimensions 
is obvious. 



2.1.1 Langevin equations and discretization rules 

As a standard paradigm for Brownian motion, we consider the one- dimensional 
motion of a point-like Brownian particle (mass M), which is surrounded by a 
stationary homogeneous heat bath consisting, e.g., of smaller liquid particles 
(mass m <C M) at constant temperature T. The inertial rest frame I S of 
the heat bath will be referred to as lab frame hereafter. The position of the 
Brownian particle in S at time t is denoted by X{t) and its velocity is given by 
V{t) := dX{t)/dt. The associated nonrelativistic momentum of the Brownian 
particle is defined by P{t) := MV(t). 



18 



By definition, the mean velocity of the heat bath particles vanishes in S. 
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According to the Langevin picture of Brownian motion, the stochastic dynam- 
ics of the Brownian particle due to the interaction with the bath and in the 
presence of a conservative external force field 3^ (t, x ) can be described by the 
differential equations H H H S, H H 113, EsS 



^ = 9^(t,X)-a(P)P + £(P,t), (3b) 



complemented by the initial conditions I X(0) = xq and P(0) = po- The 
second term on the rhs. of Eq. ( ]3bl) is a friction force with the shape of the 
friction coefficient function a{p) > depending on the microscopic details of 
the particle-bath interaction. The stochastic Langevin force 



L{P{t),t) = [2DiP)f'Qat) (3c) 

reflects fluctuations in the surrounding heat bath. The symbol © signals the 
choice of a specific discretization rule to be discussed in more detail below. 
The amplitude of the fluctuating force £ is tuned by function D{p) > 0. For a 
spatially inhomogeneous heat bath the functions a and D would also depend 
on the position coordinate x. The stochastic 'driving' function ({t) is often 
taken to be a Gaussian white noise, i.e., ({t) is characterized by: 



(CW) = o, 
(C(t)C(^)) = 5(t 



(4a) 
(4b) 



with all higher cumulants being zero. In Eqs. (jlj), the bracket ( ■ ) symbolizes 
an average over all possible realizations of the noise process C(^)- 



In the mathematical literature T4, 1761], SDKs like the Langevin Eq. ([3]) are 
usually written in the differential notation 



dX{t) = (P/M) dt, 

dP{t) = 3^(t, X) dt -a{P)Pdt 



[2D{P)Y/^QdB{ty, 



(5a) 
(5b) 



Here, dX(t) := X(t + dt) —X(t) denotes the position increment and dP(t) : = 
P{t + dt) — P{t) the momentum change. The random function B{t) is a 
standardized one-dimensional Brownian motion or, equivalently, a standard 
Wiener process jHH, 74, 76, 77], whose increments 
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Without loss of generality we fix the initial time to = 0. 
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dB{t) := Bit + dt) - Bit) 



(5c) 



are defined to be stochastically independently and characterized by the Gaus- 
sian probability distribution 



P{dfi(t) G [y^y + dy]} = (27rdt)-^/'exp[-yV(2dt)l dy 



(5d) 



i.e., the increments dB{t) are independent random numbers drawn from a 
normal distribution with variance dt. From Eq. fl5dl) and the independence of 
the increments at different times s 7^ t, it follows that 



(dS(t)) = 0, {dB{t)dB{s)) = { : J- (6) 




where now the expectation ( ■ ) is taken with respect to the probability measure 
P of the Wiener process Bif). The two different representations ([3]) and ([5]) 
may be connected by formally identifying 

d5(t) = C(t) dt 

In the remainder, we will write SDEs primarily in the differential notation of 
Eq. ([5]), which m ay a lso be viewed as a simple numerical integration scheme. 



cf., e.g., Ref. [89|, |420| and Appendix |Al 



It is worthwhile to summarize the physical assumptions, implicitly underlying 
Eqs. (E]) and dSD: 

• The heat bath is spatially homogeneous and stationary; i.e., relaxation pro- 
cesses within the heat bath occur on time scales much shorter than the 
relevant dynamical time scales associated with the motion of the heavy 
Brownian particleJ^ 

• Stochastic impacts between the Brownian particle and the constituents of 
the heat bath occur virtually uncorrelated. 

• On a macroscopic level, the interaction between Brownian particle and 
heat bath is sufficiently well described by the friction coefficient a and the 
stochastic Langevin force L. 
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This means that the joint probabiUty density of an arbitrary collection of subse- 



quent increments dB{ti) is a product of the Gaussian distributions see, e.g., [7J, |76| 
for a precise mathematical definition. 

Interaction with a spatially inhomogeneous non-stationary heat bath can be mod- 
eled, e.g., by considering friction and noise amplitude functions of the form a(t, x,p) 
and D{t, x,p). 
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• Equations ([3]) and ([5]) hold in the lab frame S, corresponding to the spe- 
cific inertial system, where the average velocity of the heat bath particles 
vanishes for all times t. 

In Section 12.21 we shall review how stochastic dynamical equations similar to 
Eqs. ([3]) and ([5]) can be derived or, at least, motivated by means of specific 
microscopic models. 



Discretization rules A stochastic force with a momentum dependent noise 
amplitude function D{p) as in Eq. fl5bl) is usually referred to as 'multiplica- 
tive' noise. When considering SDEs that contain multiplicative noise terms, 
the specification of the discretization rule © is necessary since, for fixed func- 
tions a{p) and D[p), different discretization sche mes in general lead to non- 
equivalent stochastic processes [tI, 76, 77, 422 1; put differently, the values 
of the stochastic integral P(t) defined by Eq. (l5bl) depend on the choice of 
discretization rule. This is the most essential difference compared with ordi- 
nary differential equations, whose integral curves (i.e., solutions) are indepen- 
dent of the underlying discretization scheme when taking the continuum limit 



dt — > Olfil The three most commonly considered discretization rules are the 



following ones |424 
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*", is defined 



The pre-point discretization of Ito [60, |61[, denoted by "0 
by computing the function C(P) = [D{p)Y^'^ at P{t), i.e., 

C{P) * dB{t) := C{P{t))dB{t). 

The post-point rule "0 = •", sometimes also referred to as kinetic |ll2l 



(7a) 



or backward Ito rule [42J], is defined by evaluating the function C{P) 
[P)(p)]V2 at P{t + dt), i.e.. 



C(P) • dP(t) := C(P(t + dt))dP(t). 



(7b) 



The mid-point rule "0 = o" of Stratonovich 67, 68, 69] and Fisk 6^, 66] is 
defined by taking the mean value of the Ito and the backward Ito stochastic 
integral, i.e.. 



C(P) o dB{t) := ^[C(P) * dB{t) + C{P) • dB{t)]. 



(7c) 



22 For example, the discretization rule is irrelevant when the driving process B{t) is 
a regular function, e.g., if B{t) = sin(ti;t) such that dBit) = cos{ujt) di. If, however, 
the driving process is a strongly fluctuating function (of unbounded variation) as, 
e.g., a Wiener (white noise) process, then different integration rules may yield non- 
equivalent trajectories; cf. Problem 2.29 in 74 1. 

A brief summary of the different discretization rules and their consequences with 
regard to stochastic differential calculus is given in Appendix Rl 
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From the mathematical point of view, the choice of the discretization rule 
reduces more or less to a matter of convenience 7^, |76| . To briefly illustrate 
this, consider Eqs. (jSD with post-point rule "0 = •", reading 



dX{t) = (P/M) dt, (8a) 
dP{t) = [S'it, X) - a{P) P] dt + [2D(P)]i/2 . (8b) 

Then, for each pair of sufficiently smooth functions {a{p), D{p)), one can de- 
termine new friction coefficients ao|*(p) such that the pairs (ao|*(p), -D(p)) 
describe exactly the same stochastic dynamics as Eq. ([H]) when combined with 
the corresponding discretization rule o and *, respectively. Specifically, one 
can replace Eq. (]8bl) by the equivalent Stratonovich-Fisk SDE 



dP{t) = [3^(t, X) - ao{P) P] dt + [2D(P)]i/2 o (9g^) 

ao{p):=a{p)-D'ip)/{2p), (9b) 

or, alternatively, by the equivalent Ito SDE 

dP{t) = [3^(t, X) - a^P) P] dt + [2D(P)]i/2 * dB{t), (10a) 

a.{p):=a{p)-D'{p)/p, (10b) 

where D'{p) := dD{p)/dp. The modified friction coefficients in Eqs. (l9b|) and 

(llObp account for the fact that the three discretization rules are characterized 
by different conditional expectations, respectively: 

([2P)(P)]^/2 . ^^(^) I p(^) =p) = D\p) dt, (11a) 

([2P)(P)]i/2 o ^p(^) I p(^) =p) = D'{p) dt/2, (lib) 

([2P)(P)]^/2 * dB{t) I P{t) =p) = 0. (11c) 



From a practical point of view, each of the three above mentioned discretiza- 
tion methods has its own merits and drawbacks: Ito's pre-point rule (*) is 
particularly convenient for numerical simulations, but care is required when 
considering nonlinear transformations G{P) of the momentum coordinate due 
to modifications of the differential calculus, cf. Appendix \M By contrast, if 
one adopts the Stratonovich-Fisk mid-point rule (o), then the transformation 
rules from ordinary differential calculus carry over, but it becomes more dif- 
ficult to implement this mid-point rule in numerical simulations. The latter 
disadvantage also applies to the post-point rule employed in Eq. (l5b|) . How- 
ever, as we shall see next, the post-point rule (•) leads to a particularly simple 
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form of the fluctuation-dissipation relation 
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Fokker-Planck equation When studying SDEs of the type ([5]), one is 
typically interested in the probabihty 

f(t,x,p) dxdp 

of finding the Brownian particle at time t in the infinitesimal phase space 
interval [x,x + dx] x [p,p + dp]. The non-negative phase space probability 
density function (PDF) f{t,x,p) > of the Brownian particle is normalized 
at all times, i.e. 



dxdpf{t,x,p), Vt>0; (12) 

where, here and below, unspecified integrals range over the full phase space, 
position space, or momentum space, respectively. Given the phase space PDF 
f{t, X, p), the marginal momentum PDF (j){t, p) and the marginal position PDF 
g{t,x) are defined by 



(f){t,p) = J dx/(t,x,p), (13a) 
g{t,x) = / dpf{t,x,p). (13b) 



Deterministic initial data X{0) = xq and -P(O) = po translate into the initial 
conditions 



f{0,x,p) 
g{0,x) 



■ 6{x — Xq) 6{p 

■■6{x - xo). 



Po), 



(14a) 
(14b) 
(14c) 



Adopting the post-point rule (•), the FPE describing the phase space density 
f{t,x,p) of the stochastic process ([5]) is given bju^ 



dt M dx ' dp dp 



a{p)pf + D{p) 



dl 
dp 



(15) 



The general conversion formulae for the different discretization rules are summa- 
rized in Appendix [a! 

If we had considered Eq. ()5bp with another stochastic integral interpretation 
(e.g., pre-point or mid-point discretization), then the corresponding FPE would be 
different from Eq. ()16p . 



18 



Equation (fT5|) is a linear partial differential equation in /, i.e., more gen- 
eral general solutions can be obtained by integrating the special solution 
with deterministic initial condition fll4l) over some arbitrary initial distribu- 
tion /o(a^O)Po)- Moreover, the FPE (IT^ is of first order in time, reflecting the 
fact that the Langevin equations ([5]) describe a Markovian process. 



If there are no external forces present, i.e., 3^{t,x) = 0, then Eq. 
the following FPE for the momentum PDF (f)(t,p) 



T5|) yields 



90 d 
dt dp 



a{p)p(f) + D{p) 



dp 



(16) 



The stationary solution of Eq. (1161) is given by 



^(p) = !Nexp 



P* D{p') 



(17) 



where !N" is a normalization constant, and some arbitrary constant such that 
the integral in the exponential exists. 

The general form (ITTj) of the stationary solution implies that one may generate 
arbitrary momentum distributions (e.g.. Maxwell, Bose, Fermi or power law 
distributions) by choosing the f rictio n and noise amplitude functions a{p) and 
D{p) in a suitable manner [l^, 425 |. To briefly illustrate this, consider some 
normalized target PDF 0(p) > 0. We would like to fix the relation between a 
and D such that the stationary solution 0oo(p) coincides with </)(p). Equating 
(j){p) with 0oo(p) from Eq. f|T71) . taking the logarithm and differentiating with 
respect to p we find the conditioE 
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«(P) _ d 
^p=--log0(p). 



(1^ 



For instance, if a classical nonrelativistic Brownian particle is in thermal equi- 
librium with the surrounding heat bath, then is a Maxwell distribution 



Mp) = (27rM//5)-^/%xp[-/5p7(2M)] 



P := (A^bT)- 



(19) 



where 7 is the temperature of the heat bath and k-Q the Boltzmann con- 
stant. In this case. E g. (flSl) reduc es to the the generalized Einstein fluctuation- 
dissipation relation [25|, Ill2l . Ill3 



If reexpressed in terms of the corresponding Ito or Stratonovich-Fisk friction 
coefficients a^,\o{p) from Eqs. (|10bp and (|9bp. then the derivative D'{p) enters the 
Ihs. of Eq. (fT8|) . which thus would take the form of a differential equation. 
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D(p) = a{p) M^bT = a{p) M/(3. 



(20) 



The fluctuation-dissipation relations (ITS]) and (|2UI) do fix only one of the two 
coefficients a{p) and D{p). Put differently, one is still free to adapt, e.g., the 
function a{p) such that the stochastic process fl5b|) exhibits the correct relax- 
ation behavior. This freedom is a main reason why the Langevin approach 



is successfully applicable to a wide range of thermalization processes |112 



Phy sical l y reasonable expression s for a{p) may be deduced from kinetic the- 
ory 315l . ISld . l426l . 14271 . l428l . l429l | or microscopic Hamiltonian models that take 
into account the interactio ns as well as the statistical properties of the heat 



bath |2l|, S, H H H H liij. Examples will be discussed in Section [221 



With regard to our subsequent discussion of relativistic Brownian motions, 
it will be important to keep in mind that the nonlinear Langevin equations 
like Eq. fl5bl) provide a tool for constructing Brownian motion processes with 



arbitrary stationary velocity and momentum distributions |19l . 1425 



2.1.2 Nonrelativistic Ornstein- Uhlenbeck process 

The standard paradigm for a nonrelativistic Brownian motion process in the 



absence of external forces is the classical Ornstein-Uhlenbeck process |46 
corresponding to constant coefficients 



D{P) = Do , a{p) = ao, 

yielding, e.g., the Ito SDE 



(21) 



dX = (P/M) dt, 

dP = S'it, X) dt-aoPdt + (2Do)^/^ * dB{t). 



(22a) 
(22b) 



Note that in this particular case the choice of the discretization rule is not 
relevant when integrating the momentum equation fl22bl) . but a rule must 
be specified when considering nonlinear transformations G{P)] cf. remarks in 
Appendix |Al 



Free motion Considering free Brownian motions with 5'(t, x) = ffist, the 
solutions of Eqs. read explicitly 



X(t)=X{0)+ f'dsP{s)/M, 
Jo 

Pit) = P(0) e-°o* + (2P)o)^^^e-°o* T e"°' * dB{t 

Jo 



(23a) 
(23b) 
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Combining the solution (l23l) with Eq. (15dp . one finds for the first two moments 
of the momentum coordinate 46|, |48 



(P(t))=P(0) e-"°*, 

P(t)2\ = p(o)2 e-2"ot + (1 _ e-2"o*)^ 



(24a) 
(24b) 



while the first centered moments of the position coordinate are obtained as 



(X(t)-X(O)) 

[x{t)-x{o)r 



m 

aoM 

_2Do_ 

{aoMy 



(1 - e-"°*) 



1-e 



aoM 

Do ( 
^r-^ (-3 + 4e" 



+ 



(25a) 



). (25b) 



The asymptotic spatial diffusion constant Dqo, not to be confused with the 
noise amplitude Dq; is usually defined by 



21) 



oo.- lim- ([X(t)-X 



(26) 



From Eq. ( ]25bl) we find for the Ornstein-Uhlenbeck process the classical re- 
sultf^l 



(27) 



The FPE governing the momentum PDF (f){t, p) of the free Ornstein-Uhlenbeck 
process reads |287 



at \ op 



Adopting the deterministic initial condition (114b[l , 0(0 , p) 
time-dependent solution of Eq. fl28|) is given by [461. 287| 



(28) 

Sip - Po), the 



0(t,p) 



1/2 



ao[p-Po exp(-aot)]^ ] 



27rP>o[l - exp( 



-2aot)] j ''''^1 2P)o[l - exp(-2aot)] j' ^^^^ 



A useful integral formula for the diffusion constant for nonlinear one-dimensional 
Brownian motion processes was derived by Lindner [25| recently, see Eq. (|149p below. 
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In the limit t ^ oo this solution reduces to the stationary Gaussian distribu- 
tion 



For a given momentum distribution (f){t,p) of the Brownian particle, the cor- 
responding velocity PDF ipij:, v) is defined by 



ilj{t,v):=^<l){t,p{v)), (31) 

where p = Mv in the nonrelativistic case. Hence, by imposing again the Ein- 
stein relation from Eq. (!20|) . which now reduces to Dq = aoMk^^, the sta- 
tionary momentum PDF (l30l) is seen to be equivalent to Maxwell's velocity 
distribution 



^„M.|^V'%.pf_^V (32, 



Moreover, by virtue of the Einstein relation Dq = aoMfceT, the asymptotic 
diffusion constant from Eq. flTTI) can be written in the form 



Doo := kB7/{aoM). (33) 



External force fields The FPE describing the phase space density f{t, x, p) 
of the stochastic process ( !22l) in an external force field 3'(t, x) reads 



Models of this type have been intensely studied during the past century, cover- 



ing a wide range of application (see, e.g., Ref. [78| for a review). However, for 
arbitrary time and position dependent force fields 5'(t, x) it is generally very 
difficult, and in many cases even impossible, to find exact time-dependent 
solutions of the Fokker-Planck equation 0341) . In the simpler case of a time- 



independent, conservative force field 3'(t, x) = F{x) with confining I potential 
<^{x), i.e.. 



Conventionally, a potential $(x) is called 'confining' if it increases sufficiently fast 
for |x| ^ oo so that the phase space PDF / is normalizable. 
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F{x) = -^^x) 



(35) 



one can determine the stationary solution attained in the hmit t —>■ oo. Im- 
posing as above the Einstein relation Dq = aoMfce^; the st ationary solution 
of Eq. flMl) is given by the Maxwell-Boltzmann distribution 2871 . 1288 



f{x,p) = Z exp<-f3 



2M 



+ $(x) 



f3 := {k^r)-\ 



(36) 



where the normalization constant Z is determined by Eq. ( fT2l) . 



Another important class of applications concerns time periodic force fields, sat- 
isfying 3^{t, x) = 3^{t + At, x) for some fixed period At. In this case it is some- 
times possible to derive approximate asymptotic solutions of the FPE ( !34l) 
by considering the limit t — >■ cx). These asymptotic solutions are usually also 
time periodic and can exhibit phase shifts. They may give rise to a number of 
inte r esting phenomena such as, e .g., stochastic resonance [43ll . |432| . |433| . 1434 

43i . H, mi H, H, H, H 



From the purely mathematical perspective, SDEs define well-defined models 
of stochastic processes 7^; from the physicist's point of view, their usefulness 
for the description of real systems is a priori an open issue. Hence, before 
directing our attention to the relativistic case, it is worthwhile to recall how 
nonrelativistic Langevin equations can be justified by means of microscopic 



models [21|, 190, 191|, M, M, M, M, M 



2.2 Microscopic models 



When considering Langevin equations of the type ([5]), one may in principle 
distinguish between the two following tasks: 

(a) One can postulate Langevin equations as phenomenological model equa- 
tions, study the mathematical consequences and compare the resulting pre- 
dictions with experiments in order to (in)validate the theory. Adopting this 
approach, the parameters and the explicit functional form of the friction 
and noise amplitude functions have to be determined from experimental 
data |442) . 

(b) Alternatively, one can try to motivate and derive Langevin equations from 
microscopic models. If successful, this approach yields explicit expressions 
for the friction and noise functions in terms of the microscopic model pa- 
rameters. 
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The remainder of this section addresses the latter problem, which has attracted 



co nside r able intereste d over the past decades [16|, |21, 1901, i9ii, m m m m 



96l. l4li . l428l . l430l . 443|. Langevin equations provide an approximate stochastic 



description of the 'exact' microscopic dynamics. Hence, in order to derive SDEs 
from, e.g., Hamilton mechanics one has to impose certain approximations. 
These approximations determine the range of applicability of the Langevin 
approach. Generally, one can pursue at least two different routes for deriving 
SDEs of the type ([5]) from more precise models: 



(1) Starting from a Boltzmann-type equation 2891 . 13121 . l444j or master equa- 



tion 



75l | for the one-particle probability density of the Brownian particle, 
one can try to reduce these integro-differential equation s to a Fokker-P l anck 
equ a tion by perfor r ning suitable approximations [16|, l315l . l316l . l415l . |427 . 



4281 . 14291 . 14431 . 14451 . 14461 . 1447j . Once the Fokker-Planck equation has been 
found, it is straightforward to write down a corresponding Langevin dy- 
namics 75|, |76|]. The microscopic collision dynamics is then encoded in the 



scatterin g cross-sec tions appearing in the collision integral of the Boltzmann 
equation j312l . |444 . 

(2) Alternatively, one may start from a microscopic (e.g., Hamiltonian) model 
describing the interaction between Brownian particle and heat bath. After 
eliminating the heat bath deg rees of freedom from the equa tions of motion 



for the Brownian particle M, H H S, M, M, H 



one obtains a 

generalized Langevin equation which may be reduced to the form ([5]) in cer- 
tain limit cases. As a byproduct, the fluctuation-dissipation relations arise 
quite naturally within this approach upon assuming a probability distribu- 
tion for the (initial) bath configuration. 



To briefly illustrate the second procedure by example, we next consider the 
oscillator model 
model 
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90l . l9ll . |92| . l93l . l94j . l95l . l96l | and the elastic binary collision 



In contrast to the more frequently studied harmonic oscillator 
model, the collision model from Section [2.2.2l can be extended to the relativistic 
case; cf. discussion in Section 1^751 



2.2.1 Harmonic oscillator model 



The harmonic oscillator model presents the classical paradigm for constructing 
a generalized Langevin equation from a Hamiltonian model BHHHH 



94, l95j96 

reads 96 



The Hamiltonian function upon which the derivation is based 



H 



2M 



Pr 



2m, 



+ 



nirUt. 



(37) 
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Here M, X and P are the mass, position and momentum of the Brownian par- 
ticle and is an external potential field; Xr and Pr denote the position and 
momentum of a heat bath particle with mass rrir, oscillator frequency and 
coupling constant c^. Equation (157|) yields the following Hamilton equations 
of motions: 



X = {P/M), P = F{X)+J2cr(^ 

Xr = {Pr/mr), pr = 




rrirUJj.Xr + CrX 



(38a) 
(38b) 



where F{x) = — d$(x)/dx is the conservative external force acting on the 
Brownian particle. As evident from Eqs. (155]) . Brownian particle and heat bath 
are coupled via linear forces in this model. By formally integrating Eqs. fl38bp 
and inserting the solutions into Eq. fl38al) . one may eliminate the heat bath 
coordinates from Eqs. fl38ap . yielding the exact generalized Langevin equa- 



tions 



93, 96 



X- 
P- 



■ {P/M), 
-F{X)- 



ds u{t - s)P{s) + L{t), 



(39a) 
(39b) 



where, for given initial values X(0),P(0), {0:^(0), ^^(0)}, the memory friction 



kernel i^(t — s) and the Langevin noise force L[t) are given by 96 



1 



COs[uJrit — s)], 



Xr{0) 



rrirUJ^ 



X(0) 



COS{UJrtj H Sm{u!rt) 



nirUOr 



(39c) 
.. (39d) 



In order to be able to characterize the properties of the noise force L{t) by 
means of an fiuctuation-dissipation relation, one still needs to impose a distri- 
bution for the initial conditions {xr.(0),pr(0)} of the bath variables. In prin- 
ciple, this initial distribution can be chosen arbitrarily. Of particular interest 
in canonical thermostatistics are equilibrium distributions of the Maxwell- 
Boltzmann type. In the case of the generalized Langevin equation fl39|) . a 
plausible choice for the initial bath distribution corresponds to the PDF 



/b(K(0),p,(0)}|X(0) = xo: 

\p.,W 



exp< 



V, , ,2 



2mr 



x,(0) 



Xq 



(40) 



where j3 = {k^7) ^ is the inverse thermal energy, T the temperature, and Z 
the normalization constant. The initial position X[0) = Xq of the Brownian 
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particle enters in Eq. (HOj) as an independent parameter, i.e., averages with re- 
spect to fh{{xr{0),Pr{0)} I -^(0) = Xq) are conditional on the initial Brownian 
particle position X(0) = Xq. Averaging the stochastic force L(t) with respect 
to /b from Eq. (HQ]) , one finds 



(L(t)), = 0, (41a) 
{L{t)L{s))^ = MkB7iy{t-s). (41b) 

Equation fl41bp represents the fluctuation-dissipation relation for the gener- 
alized Langevin equation (l39l) given the initial bath distribution (HOl) . The 
generalized Langevin equation (139!) differs from Eqs. ([5]) and fl22l) through the 
memory friction z/(t — s). The SDE (!22|) . describing the Ornstein-Uhlenbeck 



process in an external force field, is recovered from Eqs. (1391) in the limit casd 



u{t- s) = 2aoS{t- s), (42) 



where a is a constant friction coefficient. The limit case fl42|) can be illustrated 
by rewriting the friction kernel (I39cl) in the more general form 



p{t-s)= duj Ciuj) cos\uj{t - s)]. (43) 
Jo 

By fixing the amplitude function C{(jj) as 



C{uj) 



5{UJ — UJr 



(44) 



one recovers the memory friction fl39cl) as a special case of Eq. fB3l) . In order to 
obtain the limit case fH2l) from Eq. fj43l) . one can use the cosine-decomposition 
of the Dirac ^-function 



I poo rco 

S{t-s) = — duj e*^(*-^) = - / dco cos[uj{t - s)]. (45) 

27T J-oo TX Jo 

Hence, upon comparing Eqs. (1451) and fH51) . the white noise limit (1421) corre- 
sponds to the particular choice 



C{uj) = (2«o)/vr. (46) 



The prefactor '2' is required in Eq. (j42|) because of the convention 
/o*dsd(t-s)P(s) = P(t)/2. 
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The harmonic oscillator model provides a useful microscopic justification 
for the Langevin equations ([5]) and (!22l) of the Ornstein-Uhlenbeck process. 
Unfortunately, this model cannot be transferred to special relativity, as it is 
based on instantaneous harmonic interactions-at-distance which violate fun- 
damental relativistic principles. Therefore, in the last part of this section we 
shall consider another microscopic model which is based on strictly localized 
elastic binary collisions and, thus, can be extended to special relativity. 



2.2.2 Elastic binary collision model 



The one-dimensional elastic binary collision model |2l|, |448[ | is based on the 
idea that the stochastic motion of a Brownian particle (mass M) is caused _b 
frequent elastic collisions with smaller heat bath particles (mass m <^ M) 
As before, we denote the coordinates and momenta of the heat bath particles 
by {xr,Pr}, where r = 1, . . . , N and ^ 1. 



Collision kinematics An elastic collision of the Brownian particle [velocity 
V, momentum P = MV, kinetic energy E = P^/ (2M)] with a heat bath 
particle [velocity Vr, momentum pr = mvr, kinetic energy = p^/(2m] is 
governed by the energy-momentum conservation laws 



E + er = E + er 



P+Pr = P + Pr 



(47) 



Here, hat-symbols refer to the state after the collision. Taking into account 
the kinematic conservation laws (H7|) . we find that the momentum gain APr 
of the Brownian particle per single collision is given by 



APr =P-P = 



2m 



■P + 



2M 



M + m M + m 



Pr- 



(4J 



To construct a Langevin-like equation from Eqs. fj47|) and fj48|) . one considers 
the total momentum change SP{t) = P{t + 6t) —P{t) of the Brownian particle 
within the 'mesoscopic' time interval [t,t + 6t], assuming that: 



The quantum mechanical generalization of the harmonic oscillator model rep- 
resents a_paradigm for quan t um Brow nian motions and has been studied, e. g., in 



3 uan t 
103, [loi, [loij; see also the reviews by Grabert et al.lllij 
and Hanggi and Ingold [108l |. 

Similar approaches are known from unimolecular rate theory, see, e.g.. Section V 
in [t^. In the context of quantum Brownian motions, a quantun i-me chanical version 
of the collisio n mo del was proposed and studied by Pechukas 1051 ]. and Tsonchev 



and Pechukas lod j. 
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• collisions occurring within [t, t + St] can be viewed as independent events; 

• the time step 6t is sufficiently small, so that there occurs at most only one 
collision between the Brownian particle and a specific heat bath particle r 
and that \6P{t)/P{t) \ < 1 holds true; 

• 6t is still sufficiently large, so that the total number of collisions within 6t 
is larger than 1. 

These requirements can be fulfilled simultaneously only if m <^ M holds. With 
the above assumptions, one can approximate 



N 

6Pit) ^Y^^^r ^r{t,St), (49) 
r=l 

where Ir{t,St) G {0, 1} is the indicator function for a collision with the heat 
bath particle r during the interval [t,t + 6t]] i.e. 



^ , , , 1 1 if a collision has occurred in [t, t + 5t\, , , 

Ir{t,dt) = \ L' J' (50) 

I U otherwise. 

In the one-dimensional case, /^(t, 6t) can be expressed in the fornf^ 



Ir{t, 6t) = e(X - Xr) e(x; - X') + 0(x, - X) &{X' - x'^), (Sla) 

where X = X(t), Xr = Xr{t) are the 'initial' positions of the colliding particles 
at time t, and 



X' = X + V 5t , x'^ = Xr + Vr5t (51b) 



their projected positions I at time t + 6t. The collision indicator from Eq. flSTl) 
is characterized by 



/,(t,0) = 0, [Ir{t,6t)y = Ir{t,6t), J = 1,2,.... (52) 

The Taylor-expansion of Ir at 5t = reads [2l| 

The Heaviside- function Q{x) is defined as the integral over the Dirac 5- function, 
i.e., Q{x) := 0, x < 0; 0(0) := 1/2; @{x) := l,a: >0. When considering higher 
space dimensions, the expression (|5ip for the indicator function has to be modi- 
fied accordingly, e.g., by taking into account the geometric shape of the Brownian 
particle. 

Of course, in the case of a coUision the position of the particles at time t + 6t wiU 
be different from the 'projected' positions X' and x'. 
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Ir{t, 6t) ^^\Vr- V\ 6{Xr - X). 

Combining Eqs. (gH]), (gQ]) and flSSaj) yields 



(53a) 



m 



■ N 

y 

M + m 



Ir{t,6t) 



N 



p{t)+2Y: 



M 



{M + m 



Prlr{t,6t), (53b) 



where, additionally, it was assumed that for each collision occurring within 
[t,t + 6t], the momentum of the Brownian particle before the collision is ap- 
proximately equal to the 'initial' value P{t). In view of m <C M, Eq. (I53bp 
can be simplified further to give 



N 

P{t) +2j2pr Irit,St). (53c) 
r=l 

A comparison with the Langevin equation fl5bp suggests that, heuristically, the 
first term on the rhs. of Eq. ( 153cp can be interpreted as a 'friction' term, while 
the second term represents 'noise'. However, although looking quite similar to a 
Langevin equation, Eq. f l53cl) is still considerably more complicated than, e.g., 
the Langevin equation (]5b|) . This is due to the fact that the collision indicators 
Ir{t,St) from Eq. fl53al) depend not only on the Brownian particle's position 
and velocity but also on the stochastic bath variables {xr, Vr}- Nevertheless, it 
is possible to calculate the statistical properties of the momentum increments 
SP{t) from Eqs. provided one specifies a distribution for the heat bath 
particles. 



dPit) 



^ m 



Bath distribution In principle, one can use Eqs. (!53|) to calculate the sta- 
tistical moments ((^-P)"')^ for an arbitrarily given heat bath PDF {{xr,Pr})- 
Here, we shall focus on the situation where the (infinitely large) heat bath is 
given by a quasi-ideal gas which is in thermal equilibrium with its environ- 
ment. In this case, the one-particle PDF f^{xr,Pr) is given by the spatially 
homogeneous Maxwell distribution 



f^i^r,Pr) = (27rmA;BT)-'/' L'^ exp[-pl/{2mkj,7)], (54) 

where Xr G [0, L] with L being the one-dimensional container volume. More- 
over, we will assume that: 

• the heat bath particles are independently and identically distributed; 

• the distribution of the bath particles is not affected by the collisions with 
the Brownian particle. 
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The above assumptions can be justified for a sufficiently large bath, if collisions 
between the bath particles rapidly reestablish a spatially homogeneous bath 
distribution. 



Mean drift force We define the mean (momentum) drift as the average 
momentum chan ge ( SP)^ over the interval [t,t + St], given th e momentum 



value P at time In the case of Eq. fl53cl) . one finds [2l|, |448 



(SPit))^ = -2N j {Ir{t, St))^ P + 2N {pr Ir{t, St))^ . 



(55) 



To calculate the averages on the rhs. we note that, for a spatially uniform 
bath distribution as in Eq. fl5^ . one-particle expectation values of the form 



{G{xr, Vr) Ir{t, St))^ cslu be Calculated to first order in 6t as [21 



6t 

{G{Xr, Vr) Ir{t, St))^ =7^7 ^""r G{X, Vr) \Vr - V\ 1ph{Vr) ■ 



(56) 



Here, ^b(t'r) denotes the one-particle velocity PDF of the heat bath particles, 
which in the case of Eq. (15^ is given by the Maxwellian 



^Jhivr) = iylTi) ^'"^ e^^{-vl/vl) 



v-Q := {2k-QJ/m 



,1/2 



(57) 



By making use of Eq. (156!) . one obtains for the mean drift of the collision 
model: 



{5P{t))^^-2n^ A;bT TT- 




(5J 



where nb = N/L is the number density of the heat bath particles, 
Pb '■= Mv-Q = M(2/cBT/m)^/^ a characteristic thermal momentum scale, and 
the error function erf(z) is defined by 



erf (2;) := 




More precisely, one should write the mean momentum drift in the form of a 
conditional expectation {SP{t) \ P{t) = p)^; however, for ease of notation we write 
{SP)^ at this stage. 

Higher moments and correlation functions may be calculated in a similar manner. 
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-4 -2 2 4 



P Ipb 

Fig. 1. Nonrelativistic binary collision model. Mean drift force 'X{P) = {5P{t)/5t)^ 
from Eq. (|58p with nb = N/L denoting the number density of the heat bath particles 
and pb '■= M{2k-QT /mY/"^ the characteristic momentum of a Brownian particle 
(mass M), surrounded by heat bath particles of mass m and temperature T. 

Figure [1] depicts the mean drift force 



%{P):={5P{t)/5t)^, (59) 

obtained from Eq. ( !58l) . The absolute value of this drift force grows linearly 
for small momentum values (Ornstein-Uhlenbeck regime) and quadratically 
for large momentum values. 



In the remainder of this section, we are going to illustrate how one can use 
the result (158|) to obtain a systematic procedure for approximating the model 
equations (!53l) by a nonlinear SDE of the type (l5b|l . 



Langevin approximation Similar to the oscillator equation fl39|) . the 
Eq. fl53cp for the momentum increments SP{t) in the binary collision model 
is not yet a Langevin equation. SDEs of the type 



dP{t) = -a{P)Pdt + [2D{P)]^/^ • dB{t) (60) 

are phenomenological model equations that provide a simplified description 
of the microscopic dynamics, which in the case of the collision model is more 
precisely described by Eq. f l53cp . Therefore, to obtain a useful Langevin model, 
the coefficients a{p) and D{p) have to be chosen such that they yield the best 
possible approximation within this class of SDEs defined by Eq. (l60l) . Plausible 
general criteria for the 'best approximation' can be formulated as follows: 

The stochastic process described by Eq. (J60l) should 
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• approach the correct stationary momentum distribution; 

• exhibit the correct mean relaxation behavior. 

The first criterion is equivalent to imposing the appropriate fluctuation-dissipation 
relation on the functions a and D. For the elastic collision model considered 
here, the expected stationary momentum PDF of the Brownian particle is 
given by the Maxwell distribution 



(p) = (27rMA;BT)"^/^exp -p'/{2Mk^7) 



(61) 



According to the discussion in Section I2.1.H this implies that a and D must 
be coupled by the Einstein condition 



D{P) = a{P) Mk^7. (62) 
The second (drift) criterion can be expressed mathematically as 
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/ dP(t) 
\ dt 



P{t) = p 



(63) 



Taking into account the Einstein relation and Eq. filial) , the Ihs. of Eq. fl63l) 
is given by 



/dm 

\ dt 



P{t)=p) = -[a{p)p-a'{p) MA;bT], 



(64a) 



where a'{p) := da{p)/dp. The rhs. of Eq. (!63l) is obtained by substituting 
P = p on the rhs. of Eq. (!58|) . yielding, i.e.. 



P{t)=p) =X{p), (64b) 
/ b 

where the mean drift force % was defined in Eq. f l59|) . Hence, by virtue of 
Eqs. flMl) . we see that the drift criterion fl^ is equivalent to the following 
ordinary differential equation (ODE) for a{p): 




-a{p)p +a'{p) MkB7=X{p). (65) 

With respect to the two criteria formulated above, the solution of this ODE 
gives the friction function a that provides the 'best' Langevin approximation 

( I P{t) = p) denotes by the conditional expectation with respect to the Wiener 
measure of the Brownian motion B(t). 
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Fig. 2. Nonrelativistic binary collision model, (a) Friction coefficient 
Q^oo(p) = —^{p)/p from Eq. (I66p with nb = N/L denoting the number den- 
sity of the heat bath particles and ps ■= M{2kB7/m)^^'^ the characteristic 
momentum of a Brownian particle (mass M), surrounded by heat bath particles of 
mass m and temperature 7. (b) The ratio x{p) from Eq. (|67p as a measure for the 
quality of the approximation Uoo (p) ■ 

to the binary collision model. The initial condition fo r a(p ) must be specified 
such that the correct asymptotic behavior is obtained 448j |. Information about 
the collision model and the bath distribution is encoded in the mean drift 
force Evidently, the procedure leading to Eq. fl65l) can be generalized to 

other interaction models/bath distributions as well - provided the stationary 
distribution of the Brownian particle is known. Other types of interactions 
(e.g., nonelastic) would result in another function %{p). A non-Maxwellian 
bath distribution would affect not only the rhs. of Eq. (l65l) but also its Ihs. 
due to a modified fluctuation-dissipation relation. 

Unfortunately, it is usually very difficult or even impossible to find the exact 
analytical solution of the ODE ( |65l) for a realistic drift function 3C(p). For 
practical purposes, one can obtain useful approximations, e.g., by considering 
the asymptotic behavior for p ^ oo and p ^ 0, respectively. In the case of 
the collision model, one finds that the approximation 



a{p) ~ -X{p)/p =: aoo(p), (66) 

which becomes exact for p oo, is also applicable at small |p|-values, if 
m <C M. This is illustrated in Fig. [2] (b), which depicts the dimensionless 
ratio 



X{P) ■■= [-aoo{p)p + a'oo(p) MkB7]/%{p). (67) 

The function x is an indicator for the quality of the approximative solu- 
tion aaoip), which is plotted in Fig. [2] (a). As evident from the dotted curve 
in Fig. [2] (b), for m -C M the function x is close to unity even for small 
values \p\. This means that a^o is a good approximation to the exact solu- 
tion of Eq. (I65ll . Thus, a Langevin equation of the type (l60l) with aoo{p) from 
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Eq. (166|) and D{p) = aoo{p)MkB7 yields the correct stationary momentum dis- 
tribution (16T1) and exhibits almost exactly the same mean relaxation behavior 
as Eq. fl53p . In particular, such a nonlinear Langevin equation provides a con- 
siderably more accurate description of the Brownian motion in a gaseous heat 
bath than, e.g., a classical Ornstein-Uhlenbeck process with constant friction 
and noise coefficients, cf. Eq. fl22bl) . For instance, an Ornstein-Uhlenbeck (or 
Stokes-like) approximation is obtained by replacing aoo{p) through its mini- 
mum value 

ao := «oo(0) = -%'{0) = nb ^ (^] ^'^ (68) 

M \ nm J 

and fixing the Einstein relation D{p) = a^M k-B7 =: Dq. Adopting these 
additional simplifications, the corresponding Ornstein-Uhlenbeck process can 
be expected to provide a useful description for slow Brownian particles at 
sufficiently low temperature values T. 

Similar to the harmonic oscillator model from Section [2.2.1l the collision model 
provides useful insights into the approximations that must be made in order 
to obtain a Langevin equation from a microscopic model. Compared with the 
oscillator model, an advantage of the collision model is given by the fact that 
it can be extended to special relativity, cf. Section 14.51 below. 

From a more general perspective, the above examples illustrate which objec- 
tives can(not) be achieved by phenomenological Langevin models that are 
based on Brownian motion processes. Langevin equations of this type, and 
their corresponding Fokker-Planck equations, provide a simplified descrip- 
tion of the underlying microscopic dynamics. The coefficient functions in the 
Langevin/Fokker-Planck equations allow one to construct stochastic processes 
that exhibit the same asymptotic relaxation behavior and approach the same 
stationary or asymptotic distribution like the actual physical process. Station- 
ary distributions can often be inferred from thermostatistical (e.g., maximum 
entropy) principles, while the relaxation behavior must be deduced from the 
exact microscopic dynamics. In many cases, the resulting stochastic models 
are sufficient for comparing with experimentally accessible data, but they may 
become inaccurate for describing higher order correlations and/or the relax- 
ation dynamics far from the asymptotic state. 

2.3 Remarks and generalizations 

In this section we have focused on the most commonly considered examples of 
nonrelativistic Brownian motion processes, described by SDEs of the type ([5]). 
Such nonrelativistic Langevin equations have been studied extensively in var- 
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ious contexts over the past decades (see, e.g., the reviews |78l. Ill2[ |). The hst 
of succe ssful apphcatio ns covers a wide range of d ifferent areas inclu ding laser 
phy sics 



454, 455 



n^ . liigf , optical lattic es |450l . lisit , plasma physics |452l . |453 



high energy physics 315l . l3ld|, bi ologically and chemically moti- 



vat ed population and reaction dyna mics [456|, active Brownian motion mod- 



els |429l . |457| . l458l . l459l . l460l . l46ll . |462[ | , or theoretical and exp erimenta l studies o f 



excitation and transition phenomena in nonlinear systems |463l . |464| . |465| . |466 



The stochastic processes defined by Eqs. ([5]) share as a common feature that 
the underlying noise source is modeled by a standard Wiener process B{t). 
In general, one can also consider othe r driving proce sses such as Poisson pro- 



cesses 



74 |76| or Levy processes |467l . l468l . [469 



470| . which may give rise to 



so-called anomalous sup er- o r sub-diffusion effects; see , e.g., th e reviews by 



Bouchaud and Georges |116 



117 Furthermore, 



and Metzler and Klafter 
one can abandon the assumption (14bl) of (5-correlated 'white' noise by consid- 
ering stochastic processes that are driven by 'colored' noise, for example, by 
replacing Eq. (]4bl) with an exponential correlation 



(C(t)C(s)) = -exp(-|t-s|/rn), (69) 

where the parameter is the relaxation time of the driving noise (. The 
mathematical analysis of processes driven by colored noise is considerably 



more complicated than in the case of 5-correlated white noise (see Ref. |114l | for 
a review). The study of non-Gaussian and/or colored driving processes within 
a relativistic framework still represents an open problem for the future. By 
contrast, the relativistic generalization of stochastic processes of the type ([5]) 
has been subject of intense research in recent years, leading to the relativistic 
Brownian motion theory described in Section |H 
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Anomalous diffusion processes exhibit a nonlinear growth of the 

mean square displacement, e.g., of the for m to t'^ w ith k > 1 and k < 1 correspond- 
ing super- and sub-diffusion, respectively jllTI . |467 |. 
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3 Relativistic equilibrium thermostatistics 



The brief summary of nomelativistic Brownian motion theory in the preced- 
ing section shows that equilibrium thermostatistics plays an important role in 
constraining the relation between friction and noise coefficients in Langevin 
equations by means of suitably chosen fluctuation-dissipation relations. More- 
over, a priori knowledge about the equilibrium properties of the heat bath is 
required if one wishes to derive (generalized) Langevin equations from micro- 
scopic models. Therefore, the present section intends to summarize relevant 
aspects of relativistic equilibrium thermostatistics in order to provide for the 
subsequent discussion of relativistic Langevin equations in Section |H This part 
is structured as follows. 



Section 13.11 introduces some notation and discusses general properties of rela- 
tivistic probability density functions. Section [3^ focuses on the thermostatis- 
tics of stationary systems, since these will play the role of a heat bath later 
on. In this context, particular emphasis will be placed on the relativistic gen- 
eralization of Maxwell's distribution for the following reason: Nonrelativistic 
Brownian motion models such as the classical Ornstein-Uhlenbeck process are 
in obvious conflict with special relativity because they permit particles to 
move faster than the speed of light c. Most directly, this can be seen from the 
stationary velocity distribution, which is a Maxwell velocity distribution and 
thus non-zero for velocities \ v\ > c. The recent literature has seen some debate 
about the correct generalizat ion of Maxwell's distribution in special relativ- 



ity [22, ll96|, [1981, |199|, l232l.l474 



dynamic s sim ulations 
Jiittner 
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In Section [3^ we shall discuss recent molecular 
which favor a distribution that was proposed by 
1861 1 in 1911, i.e., six years after Einstein had formulated his theory of 



special relativity |3| . In addition to its relevance with regard to relativistic 
Brownian motion theory, the Jiittner gas also provides a useful model for illus- 
trating the subtleties of relativistic thermodynamics. The latter topic is dis- 
cussed separately in Appendix [Cl also addressing the controversy that has sur- 
rounded the Lorentz transformation law s of temperature and other thern i ody 
narnic quantit i es over the past 100 



Tii-fiii. 

149l.[l56l. 



year s 



T35I . [Tiel . [T37l . [i3i . fiiol . [moI . [mi 



121, 122 



142, 



ml. [lisl. E29I. [i3ol . IT31I . 



143, 144 



isil . [Tiol . [Teol . [Tesl . [TggI . [Teol . [ml . liTi . ITtsI . I2i3 . Eii . I2T7I . H, Sii 



1461. 147 



132, 



148, 



3. 1 Preliminaries 



Section 13.1.11 summarizes definitions and notational conventions]^^! Subse 



quently, the 'peculiar' Lorentz transformation behavior of one-particle phase 



For a more detailed introduction to special relativity we refer to Weinberg |6y or 
Sexl and Urbantke 0]. 
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space probability density functions (PDFs) will be addressed. 



3.1.1 Notation and conventions 



In special relativity, an inertial frame S corresponds to a global Cartesian 
spacetime coordinate system. A spacetime event is labeled by a (1 + d)- 
dimensional coordinate tuple x = = {ct,x) = {t,x^, . . . ,x'^) in E, where 
d is the number of space dimensions and, adopting natural units from now 
on, the speed of light c = 1. Upper and lower Greek indices a,/?,... take 
values 0,1, ... ,d, and we use Latin indices i,k, . . . G {1, . . . ,d} for the spatial 
components. Vectors with upper indices are called contravariant. 

With respect to the Cartesian coordinate frame S, the cornponents rjajs of the 
metric tensor of flat Minkowski spacetime are defined by |6| 



a = (3 = 
a = P = 1, 
a ^13. 



(70) 



By definition, the components of the covariant vector ( obtained by 

contracting the contravariant vector (x"") with rjap, i.e. 



d 
(3=0 

The tuples (x") and (xa) will be called four-vectors, regardless of the number 
of space dimensions. The corresponding four-vector gradients are defined by 



The components 7]°'^ of the inverse metric tensor are determined by the con- 
dition 



= = V (x°), (73a) 

or, equivalent ly, by 
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Throughout, we will use the Einstein summation convention defined in Eq. (j7T]) . 
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= S%, (73b) 
where 5"^ is the Kronecker 5-symbol, yielding 



rj-f" = Va/3. (74) 

The Minkowski spacetime distance between two events xa — (xX) — {tA, Xa) 
and xb — — {tB, Xb) is defined by 



d{xA:XBf ■■^rjapix'X - X%){x'^ - X^s) 

= -{tA - tsf + {xa - XBf. (75) 
By definition, the separation of two events is 

• time-hke, if d{xA^XBY < 0^ 

• hght-hke, if d{xA,XB)^ = 0; 



• space-hke, if d{xA:XB) > 0- 

Events with time-hke separation can be causally connected by (a series of) 
signals travelling slower or equal to the speed of light. Events with light-like 
separation can be causally related only by undisturbed signals travelling at 
the speed of light. Events with space-like separation are causally disconnected. 

The classical motion of a massive point-like particle through spacetime corre- 
sponds to a sufiicicntly smooth, time-like curve in S, referred to as world-line. 
Consider a stationary observer 0, who is at rest in S. It is natural that 
parameterizes the particle motion using the coordinate time t of E, i.e., 
describes the world-line as a curve {x^{t),x'{t)) with x°{t) = t. In the vicinity 
of any point (event) on the particle's world-line, an infinitesimal proper time 
differential can be defined by 

dr := {-riaf^dx^'dx^y/^ = {df - dx^^^ ^dt{l- v^^^, (76) 

where v{t) :— dx{t)/dt is the particle velocity in E. According to the pos- 
tulates of special relativity, dr is the time interval measured by an intrinsic 
clock, comoving with the particle, whereas dt is the coordinate time interval 
measured by a clock at rest in S. The four-velocity {u") of a massive particle 
is defined as the derivative of the world-line with respect to its proper time. 



d T*^ 

:= -— Uau'^ = -1. (77) 

dr ^ ^ 

For a point-like particle with rest mass m > 0, the energy-momentum four- 
vector (p") = {p^,p^, . . . ,p'^) = (e, p) is defined by 
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:= mw" =^ paP" = —m^. (7J 
Upon comparing with fl76l) . one finds for a particle with velocity v in S 



p° = e = 7717(1;), p = m^{v)v, 7(v) := (1 — v^) (79) 



Lorentz-Poincare transformations In special relativity affine-linear Lorentz- 
Poincare transformations (LPTs) of the form 

x=Ax + a ^ x'" = AV^ + «" (80a) 

describe the transition from an inertial frame S to another inertial frame S'. 
The constant four-vector a" shifts the origins of time and space, while the 
constant Lorentz matrix (A"^) may account for a spatial rotation, a chaiige 
of orientation and/or a relative velocity between the two frames S and S' [Gj. 
The matrix components A°/3 are determined by the condition 



d{x'j^,x'j^Y = d{xA,XBy ^ A''^A'^sVaf3 = VjS, (80b) 
which means that causal relations remain preserved during transitions between 



inertial systems. The LPTs (IHOj) form a grouper Of particular interest for our 
purpose, is the subgroup of proper Lorentz transformations (LTs), defined by 
a" = and the additional constraints 



A% > 1 , det(A"^) = +L 



The requirements (IHTj) exclude time reversal and space inversion. Examples 
are pure rotations 

A% = 1, A^^o = A°, = 0, ^'j = Rij, (82) 

where {Rij) is a rotation matrix [i.e., det(-Rjj) = 1 and RijRkj = Sij], and 
Lorentz boosts 



A% = 7, A\ = A^ = -^w\ A^. = 5^.+ (^-1) (83) 

with velocity w = (w^ , . . . , w^) and Lorentz factor 7 := (1 — w"^)^^^"^. To 
briefly illustrate the effect of a boost, consider a particle at rest in the spatial 

For group theoretical aspects of Lorentz transformations see, e.g., Refs. QMM- 
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origin of S and, therefore, being described by the world-hne = (t,0) 

in S. By applying the Lorentz boost (!83|) to (x") = (t, 0), we find 



x"' = A%x° = -ft = t' , x" = A'ox" = -^wH = ~wH', (84) 

which means that the particle travels at constant velocity w' = —w through E'; 
i.e., S' moves with velocity w relative to S. The inverse of the transformation 
matrix (183|) is obtained by replacing w with —w. 



From Eq. flHOj) and the definition f[78|) of the four-momentum, one finds the 
relativistic energy-momentum transformation law 



= A"^/. (85) 



Combining Eqs. (|85|1 and (I80bp . one can verify the well-known mass-shell 
condition 



2 2 2 /2 /2 2 /'an\ 

m=e— p=e— p =m, (86) 

which means that the rest mass m is a Lorentz invariant. In particular, the 
mass shell condition (186|) implies that Eq. (185|) is equivalent to a restricted 
nonlinear transformation p ^ p' = p{p'), given by 



/(p) = A*o(m2 + p2)i/2 + (87) 



5.1.^ Probability densities in special relativity 

With regard to the subsequent discussion, it is worthwhile to address a few 
subtleties concerning the definition and transfo rma t ion b ehavior of probability 



density functions (PDFs) in special relativity 15ll . 1230 



Relativistic one-particle phase space distributions To start with, we 
consider the one-particle phase space PDF f(t,x,p) > 0, where (t,x,p) are 
measured with respect to the inertial 'lab' frame S. For a relativistic many- 
particle system with conserved particle number N (e.g., a gas of i dentical 



particles), the function / can be defined operationally as follows [151| : 
An observer 0, who is at rest in S and observes the system at S-time t, finds 

Nf{t,x,p) d'^xd'^p 
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particles in the 2(i-dimensional phase space interval [x, x + da;] x [p,p + dp]. 
Assuming that the dynamics of each particle is described by functions Xr{t) 
and Pr{t) in S, the fine-grained one-particle PDF / of the A^-particle system 
is defined bjcll 



N 



fit, X,p) = iV-1 6{X - Xr{t)) S{p - Pn{t)). 



r=l 



From this definition it is evident that / satisfies the t- simultaneous normal- 
ization condition 

1 = J d'^xd'^p fit,x,p). (88b) 

Note that this integral is taken along the hyperplane "t=constant" in S. 

Alternatively, when considering, e.g., the random motion of a single Brown- 
ian particle in a fluctuating medium, the quantity f{t,x,p) d'^xd'^p can be 
interpreted as the probability of finding the Brownian particle at lab time t 
in [a;, a? -|- da;] x [p,p + dp]. In the latter case, it is usually assumed that a 
potential trajectory is realized with a certain a priori probability. Mathemat- 
ically, this idea is implemented by introducing latent variables uj in order to 



label the potential trajectories '^^1 by writing X{t] uj) and P{t; uj). The assign- 
ment of a priori probabilities is equivalent to specifying a PDF <I>(a;) on the 
set of the l atent vari ables {u}. In this case, the phase space density / in S is 



defined by [15l|, l23Q 



f{t, x,p)= J dcu $(cj) 5{x - X{t; Lu)) 5{p - P{t; to)), (88c) 

and this / is again subject to the normalization condition (I88b[) . 

Equations (IHHj) refer explicitly to the inertial rest frame S of the observer 0. 
Now consider a second observer 0' at rest in another inertial frame S' that 
moves with constant velocity w ^ relative to S. Employing an analogous 



Since the kinetic momentum P,. is uniquely linked to the velocity Vr{t) := 
dXr{t)/dt, the definition (j88ap formalizes the idea of classifying the particle curves 
Xr{t) according to their positions and time derivatives at time t. 
^■^ For example, if the particle dynamics in S is described by differential equations of 
the form dX{t)/dt = V[X{t), P{t)], dP{t)/dt = K[X{t),P(t)] with given (well- 
behaved) functions F : M'^ x M*^ ^ M*^ and K : M"* x R'' ^ W^, then a trajectory is 
uniquely determined by specifying the values X(to) = xq and P(to) = Po at some 
instant to inT,. In this case, one could choose to = {xq,Pq). More generally, uj could 
also label different realizations of some background field which affects the particle 
dynamics. 
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t 



t = const // /' 

t' = const ^ 

X 

Fig. 3. The solid curves represent two world-lines starting at the same spacetime 
point (dashed lines indicate the forward lightcone) The one-particle phase space 
PDF f{t,x,p) measures the number of world-lines that pierce through the hyper- 
plane "t=constant" within the interval [a;, a;-|-da;] while having momentum values in 
the range [p,p+dp]. Similarly, the PDF f'{t', x' ,p') measures the number of world- 
lines piercing through the hyperplane "t'=constant" within the interval [x' , x' + dx'] 
while having momenta in the range [p',p' + dp']. 

operational definition as 0, the moving observer 0' will measure another distri- 
bution f'(t', x', p') and it arises the question how the two functions f'(t', x', p') 
and f{t,x,p) are related to each other. In the nonrelativistic theory, the 
change from one inertial system to another does not affect the time coor- 
dinate; hence, one can use the standard transformation laws for PDFs in that 
case [see, e.g., Eq. ( l3TI) ]. By contrast, the situation becomes more complicated 
in the relativistic theory, because now the definition of / and /' is based on 
an observer-dependent notion of simultaneity: The measurements of and 
0' refer to the two different hyperplanes "t=constant" and "t'=constant" in 
Minkowski space, respectively. This is illustrated in Fig. [31 

In an insightful paper [151 ] published in 1969, van Kampen proved that the 
one-particle phase space PDF / transforms as a Lorentz scalar, i.e.^1 



fit,x,p) = f'it',x',p'), 



^9) 



where {t,x,p) and {t',x',p') are connected by Lorentz transformations. His 
(first) proof of Eq. (1891) merely uses an assumption about the uniqueness of 
particle trajectories and a reparameterization of the par ticles trajectories in 
terms of their invariant proper times (cf. Section 3 in 15l| ). As a consequence, 
Eq. fl89|) represents a generic kinematical result and applies to a broad class 
of interaction models. Furthermore, van Kampen showed thai 
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'^^ One can fi nd s e vera l insufficient 'proofs' of Eq. (j89p in the literature, cf. the 
discussion in [l5ll. [230|. In this context it is sometimes claimed, erroneously, that 
the phase space element d'^xd'^p is a Lorentz scalar; in Section 2 of their paper, 
Debbasch et al. [230] demonstrate that this is not true in general. 

The result (j90ap is nontrivial due the fact that the integrals refer to different 
hyperplanes "t=constant" and "t'=constant" in Minkowski spacetime, respectively. 
The proof of Eq. (j90ap follows from Gauss' theorem, by making use of the fact that 
the absolute particle velocities are bounded by the speed of light c = 1, cf. Eqs. (34) 
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J dVdV/(t',^',p') = / d'xd'pf{t,x,p), 



(90a) 



implying that the function /' satisfies the t'-simultaneous normahzation con- 
dition 



1 = J dVdV f{t',x',p'). 



(90b) 



We next summarize a few consequences of Eq. fl89l) . 



Density-current vector Equation (15^ imphes that the quantities 



git,x) = J d'^pf{t,x,p) 
jit,x) = J d'^p fit,x,p) V, 



(91a) 
(91b) 



where v = p/e = p/p^ is the particle velocity, can be combined into a density- 
current four- vector field j{x) = (j") = {g,j). As any covariant vector field, j 
transforms under a Lorentz transformation x ^— x' = Ax as 



J'(x')=Aj(A-V) = Aj(x). 



(92) 



The four-vector character of (j") becomes particularly evident by rewriting 
Eq. (I9T1) in the form 



r{t,x) 

where 



l^fit,x,p)p'^, 



(93) 



d'^p d'^p' 



p 



P 



/O 



(94) 



is invariant under Lorentz transformations 151. 312 . Furthermore, as shown 



in Section 4 of van Kampen's paper 15l[ |. satisfies the continuity equation 



(95) 



and (35) in [1511 • 



43 



stating the conservation of particle number or probability, in agreement with 
Eq. (I90al) . The kinematical proof of Eq. (!95l) as given in [151j does not require 
knowledge about the equations of motions, but uses only the existence of the 
particle trajectories. 



Energy-momentum density tensor Equation (p3|) can be generalize d to 



define a symmetric energy-momentum (density) tensor field by [l5ll . |312 



e'^^{t,x)= f(t,x,p)py^. (96) 



In particular, the "00" -component gives the (kinetic) energy density per parti- 
cle. In contrast to the particle number conservation law (!95l) . the conservation 
of the energy-momentum tensor is a dynamical property which requires knowl- 
edge about the equations of motions (see example in the next section). 

In principle, "higher moment tensors" can be defined in a similar manner: 
r^^-{t,x) = J^f{t,x,p)py^p^.... (97) 



How ever, conventional thermodynamic and hydrodynamic analysis [16J, Il77 



476l | usually focuses on relations involving and 6'"^. 



Marginal momentum distributions The zero-component g{t, x) of the 
density-current vector j(x), defined in Eq. fl91al) . can be interpreted as the 
marginal distribution of the particle positions. Similarly, one may define marginal 
momentum distributions with respect to S and S' by 

0(t,p)= I d^x/(t,a;,p), (98a) 

0'(t',p')=/dV/'(t',a;',p')- (98b) 

Again, and 0' refer to different hyperplanes in Minkowski spacetime, re- 
spectively. Thus, in general, it is impossible to calculate (p' from (p alone or 
vice versa. There exist, however, a few exceptions such as, e.g., a system of 
freely moving, identical particles (i.e., no interactions, no external fields, no 
walls). In any inertial frame, such a collection of particles is desc ribed by a 



time-dependentLzil one-particle phase space PDF / satisfying |151 
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We assume that at least two particles have different velocities. 
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f{t,x,p) = f{0,x-vt,p), (99) 

where v = p/^m^ +p^)^/^ is the velocity. Moreover, since the individual par- 
ticle momenta do not change, the marginal momentum distribution must be 
time-independent in any inertial frame. In particular, in this case - and only in 
this case - it is true that each particle observed by as having momentum p 



is observed by 0' as having momentum p', i.e., 151 



(f,{p)d'^p = 0'(p')dV, (100) 

where p'{p) is the restricted Lorentz transformation from Eq. flHTl) . Taking 
into account that d'^p' = [p'^ /p^)d'^p, cf. Eq. f lM|) . one thus obtains 



It should be stressed again that this formula holds true only in the case of an 
unconfined, non-interacting systems; it is not valid anymore in the presence 
of a confinement (see example in Section [3^ . 



Multi-particle distributions The discussion in the remainder mostly con- 
cerns one-particle distributions. Nonetheless, we mention that Eq. fl89l) _can be 



generalized to the case of A^-particle phase space PDFs /at, yielding |151 



fNituXi,p^;...;tN,XN,PN) = fNi^'v^'vPv ■ ■ ■'^'^'n^^'n^p'n)^ (102) 

where for = 1, . . . , the coordinates (t^,a;^,p'„) and (t„,a;„,p„) are con- 
nected by a Lorentz-transformation, and J'n is the multiple-time probability 
density for the lab observer to observe particle 1 at time ti near (a;i,p^), 
and particle 2 at time t2 near {x2,P2), etc.. 

The above results clarify the transformation behavior of PDFs in special rel- 
ativity, but they do not yet answer the question as to which PDF provides 
the correct description for a given physical system as, e.g., a relativistic gas 
in equilibrium. The latter problem will be addressed in the next part. 



3.2 Stationary systems with confinement 



Of particular relevance in equilibrium thermodynamics are confined systems 
that can be described by an isotropic, stationary PDF f(t,x,p) = ip{x,p) in 
a specific inertial frame E. A typical example is an equilibrated gas, enclosed 
in a container box which is at rest in the lab frame S. In standard Brownian 
motion theory such systems often play the role of the heat bath. 
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3.2.1 General remarks 



Lab frame In the lab frame S, a spatially homogeneous gas is described 
by a PDF of the form 



^{x, p) = J{x; V) (103a) 

where V = Ad(V) is the rest volume number (i.e., the Lebesgue measure |477| ) 
of the spatial container region V C M"^ in E, and J(a;; V) denotes the indicator 
function of the box, i.e.. 




" " ^' (103b) 

To be more specific, we consider a cubic container of length L such that 
V = [-1/2,1/2^. In this case, V = Ad(V) = L'^ and the indicator function 
can be expressed as 



d 

J{x; V) := n 0(^2 + x') 0(L/2 - x'). (104) 

i=l 

Isotropy and stationarity of the gas in the lab frame S require the marginal 
momentum distribution to be rotationally invariant, 0(p) = 0(|p|), yielding 
for the four-current (j°) = in S 



g{t,x) = V~^ J{x;Y), (105a) 
j{t,x) = 0. (105b) 

It is obvious that this current satisfies the continuity equation (195!) . daj°' = 0. 
Furthermore, due to the isotropy of the momentum distribution in E, the 
energy-momentum tensor (196|) becomes diagonal with components given b} 
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r (e) , a = P = 0, 

e^'^it^x) =V-^3{x;Y)l{p-v) /d, a = p = l,...,d, (106a) 

[o, ay^(3. 

Here, we have defined the one-particle mean values 

As usual, we denote hy p ■ v the ordinary Euclidean scalar product of two d- 
dimensional vectors p = (p*) and v = (v*), i.e., p ■ v := Yli=iP^''^^ ■ Moreover, we 
abbreviate p^ : = p ■ p = pip^ = Yli=iP^P^- 
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{e) := J d'^p <f>{p) e, 
{p ■ v) := J d'^p (j){p) p ■ V. 



(106b) 
(106c) 



with e(p) = (p^ + m^)^/^ denoting the energy of a gas particle, and v(p) = p/e 
the velocity. 

It is worthwhile to calculate the four-divergence of the energy momentum 
tensor: 



0, 



p = o, 

P = z, 



9,J(a;;V), 

where, in the case of a cubic box V := A(i(V) = [—L/2, L/2Y, we find 



(107a) 



di3{x; V) = S{L/2 + x') - 5{L/2 - x' 



X 



\{Q{L/2 + x^)Q{L/2- 



X- 



(107b) 



i.e., the boundaries are sources of stress [152 
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To illustrate the meaning of the energy-momentum tensor consider the 
mean (integrated) energy-momentum vector 



:= j d'^xd'^p f{t,x,p) /. 



This quantity can be rewritten in terms of 6'"^ as follows: 



(108a) 



J d'^x ^°^(t, x) = jda^ r^(t, x). 



(108b) 



The directed surface element normal to the hyperplane "t=constant" in S- 
coordinates is given by {daa) = (d'^x, 0); cf. Appendix[Bl Thus, a (i-dimensional 
spatial integration in S is equivalent to a surface integral over the hyper- 
plane "i:= const ant "in (1 + (i)-dimensional Minkowski space. For the energy- 
momentum tensor (11061) one finds explicitly 



Usually, it is assumed that the momentum conservation violation of the gas, 
f^a^gas 7^ 0, is compensated for by the energy-momentum tensor ^"((^f of the confine- 
ment (environment), i.e., dcAO^ + O'^^^^) = 0; a similar problem occurs in continuum 
models of the electron, cf. [478 1. 
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^ = »• (108c) 
/3 / 0. 

It is important to note that {p^)t is a non-local quantity, as it represents a sum 
over components of the energ y-mo ment um tensor at different spacetime points. 



As emphasized by Gamba |475j and discussed below, this aspect becomes 
relevant if one considers the question how a moving observer 0' could, in 
principle, determine {p^)t- 



Moving frame Consider an observer 0' who is at rest in a frame S' that 
moves at velocity w = (w,0, . . . ,0) along the x^-axis of S. Denoting the 
spacetime coordinates and four-momentum in S by (x") = (t,x^, . . . ,x'^) and 
(p") = {e,p^, . . . ,p'^) = (e, p) and those in S' by {x'") = {t',x'^, . . . ,x"^) and 
\p'^) = le',p'\ . . . = (e',p')> we obtain from Eqs. dgUa]) . and ([HZD the 
explicit Lorentz transformations 



X 



t' = -f{t - wx^) t = -f{t' + wx'^) 



7(-wt + a;^) x^ = -f{wt' + x'^) (109) 



e' = 7(e — wp^) e = 7(e' + wp''^) 



p'^ = 'j{—we + p'^) p^ = 'j{we' + p'^) 

x'^ = x^ pr' = p'^ j = 2, . . . ,d, 



where 7 = 'j{w) = (1 — w"^)^^^"^ is the Lorentz factor; e = (m^ and 
e' = (m^+p'^)^/^ denotes the relativistic energy, respectively. From Eqs. (11091) 
one can obtain an explicit expression for the phase space PDF /' measured by 
the moving observer. If the observer moves at three-velocity w = {w, 0, . . . , 0) 
through the lab frame, then from her point of view the container box moves 
at velocity 

w' = (w', 0, . . . , 0) = {-w, 0,...,0) = -w (110) 

through S'. By inserting Eqs. (I109p and (11101) into Eq. fl5^ for a stationary 
lab distribution of the form (I103al) . one obtains the PDF /' in the moving 
frame S' as 

f{t',x',p') = V-^ Ji^i-w't' + x"),x'\..,x"';Y) X 

<P{^{-w'e' + p"),p'\...,p"'). (Ill) 

While the phase space PDF / was stationary in S, the associated PDF /' is 



not stationary in E' due to the motion of the container boxc^ Considering 
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In fact, this also happens in the nonrelativistic case. 
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the representation fll04p of the indicator function "J in E, one can introduce a 
corresponding quantity J' in E' by 



a'(t', x'; V) ■.= J{-f{-w't' + x"),x''', x""- V) 

= 0(L/2 + 7(a;'^ - w't')) Q{L/2 - ^{x'' - w't')) x 

d 

\{Q{L/2 + x") 0(1/2 -x") 

i=2 

= Q{L'/2 + {x" - w't')) Q{L'/2 - {x'^ - w't')) x 

d 

l[eiL/2 + x")eiL/2-x"), (112) 

i=2 

where L' = L/'-f is the Lorentz-contracted length in x'^-direction. By means 
of Eq. (11121) . one can express the components of the conserved four-current 
if") = (q'J') = (7^, -7«^f?) in S' as [cf. Eq ([lOS])] 



g'{t',x') = V'-^'J'{t',x';Y) (113a) 
j'{t', x') = V'-^ J'it', x'; V) w' = g'w', (113b) 

with V := V/'y denoting the Lorentz-contracted box volume in the moving 
frame S'. By inserting the expression (11121) into Eq. (Illll) . one may perform a 
volume integration along the hyperplane "t'=constant" to obtain the marginal 



momentum density in the moving frame as |151 



</.'(p') := / dV fit', x',p') = mP'))- (114a) 

This result was already obtained by Dirac in 1924; it is sometimes written 
in the equivalent form [151.JJ 



0'(p')/^' = m/v. (114b) 

As anticipated above, Eqs. (11141) differ from the corresponding result (110 II) for 
an unconfined system of non-interacting particles. At first sight, it may seem 
surprising that the presence of the box alters the transformation properties of 
the momentum distribution. The underlying physical explanation is that the 
observations by and 0' are not synchronous and that in the time between 
their observations some particles collide with the container walls 



50 



Experimentally, the momentum PDF 4>{p) can be determined by collecting the 
momentum values p^it) of the N gas particles, measured at the same instant of 
time t in S, into a histogram. Similarly, the function 4>'{p') can be determined from 
a histogram of the momentum values p[.{t')^ measure d at the same time t' in S'. 
Cf. discussion in Section 6 of van Kampen's paper 15ll |. 
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The marginal momentum distributions present another example of non-locally 
defined quantities. In the history of relativity, such quantities have been a 
source of consi derable confusion] I Pitfalls may be avoided by keeping in 
mind that 475 : 



In general, non-locally defined quantities are not connected by Lorentz trans- 
formations if they refer to different subsets (e.g., hyperplanes) in Minkowski 
spacetime. Therefore, a complete measuring instruction for a nonlocal observ- 
able must include a statement about the set of events over which the quantity 
is defined. 

To illustrate this, we may consider the non-locally defined vector 



: / d'x'd'p' nt', x', p') p" = / d< e""{t'. 



X 



'115a) 



representing a surface integral along the hyperplane "t'=constant"l£H Al- 
though the average four- momentum vector {p'^)t' from Eq. flll5al) looks quite 
similar to {p^)t from Eq. (llOSp . the crucial difference is given by the fact that 
these two objects refer to two distinct hyperplanes, respectively, as signaled 
by the different labels t and t'. For the en ergy- momentum tensor f ll06p . the 



rhs. of Eq. f lll5al) can be evaluated to give [218 



h{{e)+w^p-v)/d), /3 = 0, 

{P^)t' = + {p-v)/d), (3 = 1, (115b) 

[o, (3>1. 

Upon comparing this result with Eq. fllOScI) . we note that {p^)t and {p'^)t' are 
indeed not related by a Lorentz transformation, i.e.. 



Gamba 4751] gives a detailed list of examples, including the Lorentz contraction 
of length: "Only rarely is it pointed out clearly that the measurements of the two 
observers ...do not refer to the same set of events. The 'ends' of the rod being taken 
as contemporary for both observer at rest and moving observer, are in fact different 
points in the four-dimensional (absolute) spacetime. Once this is clearly indicated, 
the accepted definition [of length] is not particularly harmful." 

As discussed in more detail in Appendix [Bl for a surface integral along the hy- 
perplane "t'=constant" in T,' one finds in S'-coordinates 

J^da'J""^-{t',x') = J dV0'O'3-(t>'). 
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This inequality is, in fact, caused by the non- vanishing divergence of the 
energy-momentum vector in the presence of a container, cf. Eq. (11071) . The 
actual Lorentz transformation {p'^)t of { p^)t is obt ained by keeping the un- 



derlying hypersurface "t=constant" fixed }l78l . |475 



I.e., 



l^p'f^^^ .= ^^1 d^x ^°^(t, x) = k'^f.jda^ = jda'^ (116a) 
producing the expected result 




/? = 0, 
13 = 1, 
(3>l. 



(116b) 



With regard to experimental observations this means that, in order to deter- 
mine {p'^)t, a moving observer 0' must first reconstruct the momentum values 
on the hyperplane "t=constant" from her data before being able to compare 
with the averages of a lab observer 0. 

This gives rise to a bit of a dilemma when attempting to formulate a relativistic 
thermodynamic theory that is based on global system averages such as, e.g., 
the total internal energy, momentum, etc.: 



Following Planck [121| and Einstein [122| , one can build a theory based on 
observer- simultaneously defined quantities, as those in Eq. ( lllSbj) . but this 
leads to a loss of covariance because different observers would use different 
hyperplanes depending on their velocities. 

Alternatively, one can define thermodynamic quantities with respect to a 



specific, fixed hyperplane to obtain a manifestly covariant formalism [178 



In this case, however, the question arises as to which hyperplane is the most 
appropriate one (cf. discussion below). 

These difficulties c an, i i i pri i icipl e , be avoided by fo r mula ting a thermody- 



namic field theory [l6i, [T3, mi H H, H, EtsI, EsJ in terms of local 



tensorial quantities as, e.g., and 9°'^ . In view of the fact that traditional 
nonrelativistic thermodynamics 2871 . 14811 ] intends to describe a many-particle 
system by means of a few macroscopic control parameters, one could argue 
that such a field theoretic approach is already somewhat closer to hydrody- 



namics. 



3.2.2 Juttner gas 

The above results apply to an arbitrary stationary momentum PDF 0(p) in 
the lab frame E. With regard to Brownian motion theory, one is particularly 
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interested in thermal equilibrium distributions. When postulating relativistic 
Langevin equations [ll|, [iT], [l8|, [20], these distributions must be known in 
advance in order to correctly specify the relativistic fluctuation-dissipation 
relation. Similarly, thermal equilibrium distributions are required as an input, 
if one wishes to derive Langevin-type equations from microscopic models, cf. 
Section 12.21 
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At the beginning of the last century it was commonly accepted that a di- 
lute ( quas i-ideal) gas in equilibrium is described by the Maxwellian velocity 
PDF [isj 



iljMiv;m,/3,d) = [/3m/(27r)]'^/2gxp(_/?W/2), 
or, equivalently, by the one-particle momentum distribution 



;il7a) 



MP;m,P,d) = [P/{27rm)f'exp[-Pp'/{2m)] 



(117b) 



where m is the rest mass of a gas particle, v = p/m & the nonrela- 
tivistic velocity, and T = (k^P)"^ the temperature. After Einstein [2I, [sj had 
formulated his theory of special relativity in 1905, Planck and others noted 
immediately that the distribution (11171) is in conflict with the fundamental 
relativistic postulate that velocities cannot exceed the ligh t speed c. In 1911 
a solution to this problem was put forward by Jiittner |l86| , who proposed to 
replace Maxwell's PDF by 



ijj{v;m,(3,d) = Za^m'^-f{v)^^'^ exp[-(3m-f{v)] 6(1 - \v\), (118a) 
yielding for the relativistic momentum p = mv){v) the PDF 



(.j(p; m, /?) = Zt exp[-/?(m2 + p^fl^ 



;il8b) 



Similar to the Maxwell distribution (I117p . Jiittner's distribution (I118p refers 
to a lab frame E where the vessel, enclosing the gas, is at rest. For space 



Knowledge of the relativistic equilibrium distributions is essential for the co rrect 
inte r pret a tion of experimental observations in high energy and astrophysics 318l . 

kill. Examp les include thermalization processes in heavy ion collision 
experiments 318l . |319 ^ and ultra-relativistic plasma beams 309l . 326], or the rela- 
tivistic Sunyaev-Zel'dovich (SZ) effect j41ll ]. describing the distortion of the cosmic 
microwave background (CMB) radiation spectrum due t o th e inte raction of CMB 
photons with hot electrons in clusters of galaxies 4821 . 4831 . 484]. The predicted 
strength of the spectral distortions and the cosmological parameters inferred from 
the S Z eff ect depend on the assumed velocity distribution of the relativistic elec- 
trons 411 ] . 
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dimensions d = 1,2, 3, the normalization constant 



'119a) 



can be expressed as |186 



Zi = 2m KiiPm), 

%2 = 27rm^ exp(— /?m) (1 + /5m) / (/5m) 
Z3 = 47rm3 K2{(3m)/{(3m), 



(119b) 
(119c) 
(119d) 



with Kn{z) denoting modified Bessel functions of the second kind 485| . The 
average energy per particle is obtained by logarithmic differentiation 



yielding 



(120a) 



(e)i=m 



K^iPm) + K2{l3m) 
2Ki{f3m) 
im?l3 



2 

l + m/5' 
3 , Ki(/5m) 
13^ K2{Pm)' 



(120b) 
(120c) 
(120d) 



and exhibiting the limiting behavior 



lim (e) ■ = m, 
lim (e), = d. 

(3^0 ^ 



(121a) 
(121b) 



Moreover, one can show that, for arbitrary space dimensions d the expectation 
value {p ■ v) is independent of the mass m; more precisely 



p-v) = (pVe) = d/P 



(122) 



which allows one to regard {p ■ v) a.s a. statistical thermometer (cf. discussion 
below) . 



53 



Maximum (relative) entropy principle Jiittner 186| orig i nally deri ved 
the distribution (IllSp from a maximum entropy principle 219l . |242| . |474| by 
postulating that, for a sufficiently large particle number iV ^ 1, the one- 
particle equilibrium distribution in phase space (f{x,p) be a maximizer of the 
Boltzmann-Gibbs-Shannon (BGS) entropy I 



&[<p\h-''] = -k^J^d^x J d''p^{x,p) ln[/iV(a^,P)], 



while satisfying the constraints 



:i23a) 



1 = j d'^x j d'^p^(a;,p), 



(123b) 

(123c) 



for a given mean energy value per particle (e) = 8,/N. Pursuing a similar ap- 
proach, Jiittner later also derived the correspondin g exp ressions for relativistic 
quantum gases obeying Bose and Fermi statistics jl94 
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Micro canonical approach An alternative way 1231 . 12271 . |242| of justify- 
ing the Jiittner distribution flllSp is based on the microcanonical ensemble. 
To illustrate this, consider a dilute system of weakly interacting relativistic 
particles enclosed in a vessel V C M°' of volume V = Ad(V) which is at rest in 
the lab frame S. Assume that, with respect to E, this system can - in some 
approximation - be described by a truncated Hamiltonian of the form 



H{xi,...p^) 



N 

El 



+ pI 



(124a) 



complemented by the condition of elastic reflections at the boundaries of the 
vessel. By writing fll24al) one demands that the energy stored in the interac- 
tion is negligible compared with the energy carried by the particles, so that 
£ is approximately conserved. The corresponding microcanonical A^-particle 
distribution is given by the PDF 



^^The BGS entropy &[(j)\h~'^] from Eq. (fT23al) is a relative entropy (isi. lisTl. lisi. 
489l . 14901 ] with respect to the Lebesgue measure (normalized by the constant h ) on 
the 2d-dimensional one-particle phase space Ux, p ) }; i. e.. by writing Eq. (I123P one 
has already fixed a specific reference measure |47i.l490j. 



Alte rnati ve derivations of the Jiittner distribu tion are discussed by Pauli 
Synge [Hi], Matolcsi et al. H, and Debbasch [242] . 



12$, 
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UJ 



0, as, ^ V. 



(124b) 



The normalization constant u = 00(8,, V,N) can be obtained from the inte- 
grated phase volume fl = Q{8,,V,N) by performing a differentiation with 
respect to £, i.e., 



fi:= I dTNQ{^-H{x^,...p^)), 

r dQ 
:=] dTr, 6{e-H{x,,...p^)) = —., 



(124c) 
(124d) 



where, for identical particles, the integration measure in A^-particle phase 
space is defined by 



dr 



N 



{N\h 



dN-,-1 J d a;i ■ ■ ■ d^XArd^i ■ ■ ■ d^AT, 
0, 



Xs^Y. 



(124e) 



The one-particle momentum PDF corresponding to f^'" from Eqs. (1124bll is 
defined by 



N 



MC 



(jy^M){p):=^N-'Y.6{p-p^)j 

N 

= y dr^/]^C(a;i,...,p^) N-'Y.^iP-PN)- 



(125) 



.8 = 1 



As di scuss ed in Pauli's early review |123l | and more recently also by Matolsci 
et al. 2271], in the limit of an infinite particle number —>■ 00, the PDFs 0(Ar) 
converge to the Jiittner distribution 0j from Eq. flllSbl) . 



Quite generally, a microcanonical PDF o f the form f|124bp gives rise to a 
number of rigorous equipartition theorems 287 . 2881 . HqL] such as, e.g. 



56 



MC 



UJ 



d_ 



(£,iv). 



(126) 



for each s = 1, . . . , iV and i = 1, . . . , ci. The rhs. of Eq. ffT26D defines j287l . l288 . 
the usual microcanonical temperature of the gas in the lab frame. In the 
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case of the Hamiltonian f ll24ap . the equipartition theorem ( 1126P is consistent 
with Eq. (11221) upon identifying p = (A;bT^^)"^ in the limit > 1. 



^^In this case, the summation convention is abrogated on the Ihs. of Eq. p26|) : i.e. 
Eq. ()126p holds separately for each s = 1, . . . , N and i = 1, . . . , d. 
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Computer ex periments Although t he J iittner distribution (IllSp became 



widely accepted |l2lll79l . 1188. '242. 12941311 



several modifications have been 
The valid- 



198, 19S, 232, 235, 474 



discussed in the literature |22,, ,1951 . Il96 . 
ity of Jiittner's Eqs. (11181) was recently confirmed in f ully relativistic one- 



202] . see Fig. HI as well 



dimensional (ID) molecular dynamics simulations 
as in (semi-)relativistic two and three-dimensional simulations 201 



492^ . By 



restricting the dynamics to one space dimension, it is possible simulate local- 
ized elastic particle interactions in a relativistically consistent manner without 
needing to introduce fields. This is due to the fact that, in the ID case, the 
outgoing momenta {pa,Pb) of two colliding particles A and B are uniquely 
determined by the momentum values {pa,Pb) before the collision by means of 
the relativistic energy momentum conservation laws 22 



e{mA,PA) + e{mB,PB)=^ijnA,PA) + eimB^Ps): 
Pa + Pb=Pa+Pb, 



(127a) 
(127b) 



where 6(171, p) = {m? +p^)^/^. In higher space dimensions, additional assump- 
tions abou t the spatio-temporal stru cture of the collision processes need to be 



made H, H, H, B, |273 



In order to be fully consistent with the 
requirements of special relativity, relativistic interactions d > 1 space dimen- 
sions must be formulated in terms of fields. Direct simulation of the field dy- 
namics is numerically expensive and, therefore, practically infeasible in most 
cases. Alternatively, one can use simplified semi-relativistic models such as, 
e.g., effective hard-sphere models where the interaction radius is defined with 
respect the center-of-mass frame of the colliding particles. Generally, it can 
be expected th at su ch simplified models yield satisfactory results in the low 



density regime [201|, but they may lead to inconsistencies at high densities. 



e.g., when three-body encounters become relevant. 

In addition to (in)validating theoretically predicted distribution functions, nu- 
merical simulations can be useful for illustrating the meaning of concepts such 
as 'temperature' and 'thermal equilibrium' in special relat ivity . Figure H] shows 
results of a ID relativistic gas simulation as described in fioi . The simulated 
gas consists of A''i 'light' particles (mass mi) and N2 'heavy' particles (mass 
m2 > mi). Particle collisions in the gas are governed by the conservation 
laws (11271) and collisions with the boundaries are elastic in the lab frame (i.e., 
p — >■ —p in S) . As evident from Fig. HI the numerically obtained stationary 
velocity distributions are very well matched by ID Jiittner functions (11181) 
with different masses but same parameter (3. Hence, the Jiittner distribution 
not only agrees with the numerical data, it also yields a well-defined concept 
of 'temperature' in the lab frame: Intuitively, the temperature T is thought to 
be an intensive quantity that equilibrates to a common value if two or more 
systems are brought into contact with each other (i.e., may exchange different 
forms of energy). In the example considered here, it is natural to regard the 
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particle species as two different subsystems that may exchange energy via elas- 
tic collisions. After a certain relaxation time, the combined system approaches 
a 'thermodynamic equilibrium state', where each subsystem is described by 
the same asymptotic, two-parametric velocity PDF ipj{v] rrii, j3), differing only 
via the rest masses m,. The commonly shared distribution parameter j3 may 
thus be used to define a relativistic (statistical) equilibrium temperature in E 
by 



Q 
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heavy particles 
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Fig. 4. Left column: Equilibrium velocity PDFs as measured by an observer with 
velocity w = in the lab frame S. The numerically obtained one-particle velocity 
PDFs (o) are based on a fully relativistic, deterministic ID molecular dynamics 
algorithm 2021 ] with A^i = 1000 light particles of mass mi and = 1000 heavy 
particles with mass m2 = 2mi. Considering elastic interactions between the gas par- 
ticles and elastic reflections at the boundaries, the mean energy per particle in the 
lab frame, e = 'rni'^{vr)vr + Yl^=i ''^2l{vs)vs]/{Ni + N2), remains conserved 

during the simulation. The solid curves in the upper and lower panel correspond to 
Jiittner functions (jll8a|) with different particle masses but same inverse tempera- 
ture parameter /? = 0.701 (mic^)^^. The latter can be determined from the value 
e = 2.5mic used in the simulation by the procedure described in Right col- 



umn: Equilibrium velocity PDFs as measured by a moving observer S' who travels 
at constant velocity w = 0.2c relative to the lab frame S. Parameter values and 
initial conditions are the same as those in the left column. In contrast to the di- 
agrams with w = 0, the PDFs for w = 0.2 are based S'-simultaneously measured 
velocities. The solid curves represent the theoretically expected, Lorentz-boosted 
Jiittner function 'ip'j{v';m, f3,w) from Eq. (jl29p . The very good agreement with the 
simulation data corroborates the validity of the Jiittner distribution. In particular, 
since in the stationary (equilibrium) state both particle species are characterized 
by the same /?- value, the definition 7 := (^b/?)"^ yields a meaningful statistical 
temperature concept. 
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(128) 



in agreement with the interpretation of f3 in the nonrelativistic case. 

However, for this concept of temperature to be meaningful, a restriction of the 
accessible spatial volume is required - be it by means of periodic b oundary 
conditions, by imposing reflecting walls or other types of confinements ISTfJ^ 



Otherwise, it cannot be expected that a many-particle system approaches a 
universal stationary equilibrium state independent of the specific initial con- 
ditions. 

The simulations can also be used to determine the equilibrium velocity distri- 
butions as seen from another frame S', moving with velocity w relative to the 
lab frame S, see Fig. HI In contrast to the diagrams for w = 0, the numerical 
data points for w = 0.2 are obtained by measuring the particle velocities v'^ 
simultaneously with respect to S'. The solid curves in the right column of 
Fig. m correspond to the PDF 



V-K^'; (3, w) = j{w)-^ — exp[-Pj{w)mj{v') (1 + W)], (129) 

obtained by making use of transformation law (I114p . The excellent agreement 
between the simulation data and Eq. (11291) confirms the validity of Eq. (11141) 
and, thereby, also th e sca lar Lorentz transformation behavior of the phase 



space PDF f(t,x,p) |l5ll |. Moreover, for the model under consideration, one 
may state more precisely: Two relativistic gas components are in 'thermody- 
namic equilibrium' for any observer if their one-particle velocity PDFs are 
given by generalized Jiittner functions (I129P with same parameters (3 and w. 
It is only in this case that the net energy-momentum transfer between the 
different gas components in the c onta iner vanishes (a more general version of 
this criterion is discussed in Ref. fl5n| ). 



The numerical simulations illustrate that T = (/cb/?)~^ yields a useful sta- 
tistical temperature concept in the lab frame S, defined by the boundary 
(container box). The question as to which temperature value is measured by 



The critical role of t he bo undary conditions i n relativis tic systems has been em- 
phasized by Sinyukov 493] and van Kampen 150l . 151 1- Loosely speaking, if a 



many-particle system has reached a universal stationary equilibrium state, then 
each particle 'knows' about the presence of the confinement /walls because equili- 
bration typically requires momentum reversal at the walls, e.g., in order to maintain 
a uniform density and a well-balanced average collision frequency. In particular, if 
the walls are considered to be stationary objects then they single out a preferred 
frame of reference. The relevance of the boundary conditions is even more obvious 
in quantum mechanics/statistics d ue t o their direct effect on the energy spectra 



and, thus, on the density of states [194( 1. 
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Fig. 5. Measured mean particle velocity (o, left digram) and estimated temperature 
(o, right diagram) based on Eqs. (|130p as a function of the observer velocity w. Solid 
lines correspond to the theoretically expected values, respectively, using the same 
simulation parameters and initial conditions as in Fig. [H 

a moving observer 0' cannot be uniquely answered - the answer depends on 
the respective definitions (i.e., measurement conventions) for heat, internal 
energy, etc. as employed by 0'. These aspects are discussed in more detail 
in the Appendix [Cl At this point, it suffices to mention that T := (A;b/3j)~^ 
can also be determined from t'-simultaneously sample d vel ocity data (right 



column of Fig. Hj) by means of the expectation value I 141 



kj^7 = m-f{wf (7(t;') {v' + wf)^^ , (130a) 
where 



w 



{v')t' (130b) 



is the mean velocity of the gas measured by the moving observer. One can 
test the validity of Eq. fll30ap explicitly by using simulation data obtained 
for different values of w, see Fig. [51 Hence, in this sense, Eqs. (11301) defines a 
Lorentz invariant statistical gas thermometer - when adopting this thermome- 



ter a moving body appears neither hotter nor colder [202 
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We conclude this part with another general remark: The above discussion 
has focused on the case where an 'ordinary' thermal equilibrium state is ap- 



Equation (jl30ap is obtained by combining the micro canonical equipartition the- 
orem for a Hamiltonian H = "^j^i E{mi,pi) + X^j^i -^("^2;Pj) with the Lorentz 
invariance of the relativistic phase space PDF /. 

The mean value from Eq. (jl30ap can be used to measure the rest temperature, 
which plays a central role in van Kampen's [150] approach to relativistic thermo- 
dynamics. Evidently, upon multiplying Eq. ()130ap by factors ^{w)", a ^ one can 
construct thermometers that measure 'other' temperatures; e.g., a = —1 would cor- 
respond to Planck's [121(] form ulati on of rela tivis tic thermodynamics and a = 1 to 



proposals made by Eddington [l24i | and Ott [l27 1 . 
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proached, corresponding to a stat iona r y 'ex ponential' one-particle momentum 



distribution of the Jiittner type [l86l . [IQJ]. In general, one could also imag- 



ine non-equilibrium scenarios that give rise to (quasi-)stationa ry momenturn 



distributions which differ from the Jiittner function (see, e.g, fi97l . H, liii . 



235l . l236l . 12371 . |239[ |) . When attempting to develop a thermodynamic formalism 
for such non-standa rd d istributions the entropy functional needs to be mod- 
ified appropriately 
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However, our subsequent discussion of relativistic 
Langevin equations will concentrate on the case where the heat bath, which 
surrounds the Brownian particle, is described by a spatially homogeneous 
Jiittner function. 
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4 Relativistic Brownian motion processes in phase space 



The present section describes the generalization of the Langevin theory of 
Brownian motions to the framework o f specia l relativity [ill . Il2l . Il3l . Il4l . Il5 . 
13, M, 0, S H H H H S H, H, [Ili, EJ. More precisely, 



we will consider stochastic differential equations (SDEs) that describe Markov 
processes in relativistic one-particle phase spacej^ Relativistic Langevin equa- 
tions present a useful tool for modeling the dynamics of relativistic particles in 
a random environment . Examples include the analysis of thermalizat ion effect s 



in relativistic plasmas |318l . l319l . l329l . l421| and astrophysical systems |326l . l495 



Here, we would like to provide a comprehensive introduction to the underlying 
mathematical theory. 

From a conceptual perspective, one can distinguish two complementary ap- 
proaches towards modeling relativistic stochastic processes in one-particle 
phase space: The first approach is based on postulating evolution equations for 
the (transition) probability density of the stochastic process, as e.g., integro- 
differential master-type equations or Fokker-Planck-type partial differential 
equations. Within a relativistic framework, the condition of subluminal par- 
ticle propagation imposes stringent constraints on the structure of the evolu- 
tion equatio ns. T hese asp ects were studied in the early pioneering works of 
Lopuszahski 330], Schay 332l | and Hakim 338| . Alternatively, one may start 



by postulating SDKs as phenomenological model equations and subsequently 
derive the evolution equations for the associated (transition) probability den- 
sities. Here, we are going to pursue the latter route for mainly two reasons: We 
will be able to directly compare with the nonrelativistic case discussed in Sec- 
tion [2] and, more importantly, the physical origin of the stochastic dynamics - 
the interaction with a complex background medium - is 'less hidden' within 
the Langevin approach!^ 



According to our knowledge, Debbasch et al. [ll(] were the first to propose 
relativistic Langevin equations as models for relativistic Brownian motions. 
More precisely, these authors constructed a relativistic Ornstein-Uhlenbeck 
process (ROUP) by postulating an additive Gaussian white noise force in 
the rest frame of the bath while adapting the friction such that the stationary 
momentum distribution is given by a Jiittner function. During the past decade. 



various properties of the ROUP were studied by Debbasch and Rivet |12| 
Barbachoux et al. 13, 14], and Zygadlo 2^. As will be illustrated below, the 



^'^A main motivation for considering processes {X{t),P{t)) rather than 'pure' 
position processes X{t) is that in, in the relativistic case, it is impossible to 
define nontrivial relativistic Markov proces ses X (t) in Minkowski spacetir ne; c f. 
Lopuszahski [ssd], Dudley (Theorem 11.3 in fsS^), Hakim (Proposition 2 in 338l |). 
and also the discussion in Section [5l 

For a detailed analysis of the complementary approach we refer to Hakim 338| | . 
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ROUP represents a STJecial limit case of a larger class of relativistic Langevin 

which also includes the relativistic Brownian motion 
Processes belonging to this class share 



models 05, 19, 26, 31, 32 



(RBM) process proposed in [17|, [18 



the common feature that their stationary momentum distribution is given by 
a Jiittner function, but they may exhibit a significantly different relaxation 
behavior, resulting e.g . in a different temperature dependence of the spatial 
diffusion constant 
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The present section is structured as follows. Section 14.11 introduces the ax- 
iomatic Langevin approach, i.e., suitable SDEs and fluctuation-dissipation re- 
lations are postulated in order to provide a simplified model of the complex in- 
teraction between the Brownian particle and its environment (heat bath). Spe- 
cific example processes will be analyzed in Section 14.21 Section 14.2.21 discusses 
the temperature dependence of the asymptotic mean square displacement for 
three different relativistic Brownian motion processes 15|, |25|]. It is demon- 



strated that the diffusion constant is sensitive to variations of the friction 
coefficients. This implies that measurements of the diffusion constant may re- 
veal information about the underlying microscopic interactions. In Section 1^7^ 
relativistic Brownian motion processes will be analyzed from the viewpoint 
of a moving observer. The relativistic generalization of the binary collision 
model from Section 12.2.21 is discussed in Section 14.51 representing a simple 
microscopic model for relativistic Brownian motions. 



4-1 Relativistic Langevin and Fokker-Planck equations 



When considering Langevin equations as models of Brownian motion, one im- 
plicitly assumes that it is possible and reasonable to separate the degrees of 
freedom of the Brownian particle from those of the environment (heat bath). 
Adopting this point of view, one can specify two distinguished frames of refer- 
ence: the lab frame S, defined as the inertial restframe of the heat bath, and 
the inertial frame that is comoving with the Brownian particle at a given 



instant of time|_i| To start with, we consider the question how one can model 



the dynamics of a relativistic Brownian particle in the lab frame S. 



Apart from an irrelevant shift of the origin, the inertial lab frame S is uniquely 
determined by the requirement that the mean velocity of the heat bath particles, 
which is assumed to be constant in any inertial frame, vanishes in S. Similarly, 
the instantaneously comoving frame is determined by the condition that the 
S*-velocity of the Brownian particle is equal to zero at the given instant of time. 
Generally, we assume that the time coordinates t and can be measured, e.g., by 
using atomic clocks that are at rest in S or S*, respectively. 
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4.1.1 Relativistic Langevin equations: general construction principles 

Roughly speaking, a relativistic Brownian motion process is a stochastic pro- 
cess whose absolute velocity | = |dX/dt| does not exceed the light speed 
c = 1 at any time; i.e., the process must satisfy the condition 



I^WI = = (^gi^^ < 1 V«. (131) 

where M > denotes the rest mass of the Brownian particle. Of particular 
interest here are processes can be modeled by means of SDEs similar to the 
nonrelativistic Langevin equations ([5]). The basic idea for constructing such 
processes is to couple the noise sources (e.g., Wiener processes) only to the 
relativistic momentum P = (P*), i = 1, . . . ,d, which can take values in whole 
W^. By doing so, the condition (I13ip is automatically fulfilled. However, before 
one can actually write down specific SDEs for the relativistic momentum com- 



ponents P\ an additional question needs to be addressed, namely |338l . l496l | : 



The choice of the time parameter A fundamental assumption (postu- 
late) of nonrelativistic Galilean physics is the existence of a universal time t. 
Therefore, within the nonrelativistic Langevin theory, it is quite natural to 
identify this universal time t with the time parameter of the stochastic driving 
process, which is often taken to be a mult i- dimensional Wiener process B{t). 
By contrast, in special relativity the notion of time becomes frame-dependent 
and it is important to specify in advance which time parameter is used to 
quantify the fluctuations of the stochastic driving process. 

When considering the stochastic motion of a relativistic Brownian particle, two 
characteristic time parameters can be distinguished: The coordinate time t of 
the inertial lab frame S, which may be interpreted as the proper time of the 
heat bath, and the proper time r of the Brownian particle. In principle, either 
of the two parameters could be used to formulate SDEs for the spatial compo- 
nents of the particle momentum, P = (P*). However, within the conventional 
Langevin picture of Brownian motion, one usually considers friction and noise 
as externally imposed forces that act upon the Brownian particle and reflect 
the fluctuations within the heat bath. Therefore, it seems more natural to 
characterize the statistical properties of the noise source in terms of the lab 
time t. However, as we shall discuss in Section 14.31 one can reparameterize a 
given lab-time Langevin equation in term of the associated proper-timelff. 

Within this lab-time approach, the proper time becomes a stochastic quantity, 
and one could, for example, ask for the probability of finding at lab time t the 
particle's proper time in the interval [r, r -|- dr]. Conversely, if adopting the proper 
time T as the primary deterministic evolution parameter, one could ask for the 
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Lab-time Langevin equation Within tlie lab-time approach Il2l. Il3. 

one aims at constructing t-parameterized 2(i-dimensional 
stochastic processes {X(t), P{t)} with respect to the lab frame S, where 
the position coordinates X = (X*) and the spatial momentum coordinates 
P = (P*) are connected by the standard relativistic differential relation 



dx\t) = (pyp^) dt, 



(132a) 



with P°(t) = E{t) = (M^ + p2^V2 denoting the relativistic energy of the 
particle and V^{t) = P^/P^ its velocity components in S. Stochasticity is 
implemented into the dynamics by coupliji^the momentum components P^{t) 
to an external noise source via the SDEi 
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dP\t) = Tdt - a'j P^ dt + c'j dB^{t), (132b) 
where 

• 3"* denotes a deterministic external force (e.g., Lorentz force), — a*j P^ is a 
phenomenological friction force, and the last term represents noise; 

• G {*, o, •} signals the discretization rule (i.e., the stochastic integral 
definition, cf. Appendix |A|) : 

• in general, the functions IF*, and d j can depend on (t, X, P); 

• the noise source is usually modeled by the rf-dimensional standard Wienej 
process B{t) = {B^(t), . . . , B'^{t)) in a relativistic context. 

Equations fll32p govern the spatial components of the four- vectors (X°, X) 
and (P°, P). One can still add an equation for time component X° by setting 



dX°(t) = (P7P°) dt = dt, (133a) 

The equation for the energy component P° can be derived from Eq. (]132bp 
by applying the (backward) Ito formula, cf. Appendix \^ to the mass-shell 
condition P°(t) = (M^ + P'^f/^, yielding 



probability to find the particle at proper time r within the spacetime interval [t, t + 
dt] X [a;, x+dx\ with re spect to the lab frame. The proper-time approach was studied 
in Refs. j337 . 3381 . 391 1 and most recently in 30]. 

^^E.g., for the post-point rule "0 = •", Eqs. (jl32p constitute a special case of the 
general SDK (1X28]) . upon identifying n = 2d, 1" = (X\ . . . , X^, P\ . . . , P'^) and 



= P'/P^, C'j = 0, A^+'^ 



for i 



In principle, one could also consider other driving process 
processes) but we are not aware of such studies. 



l,...,d. 

(as, e.g., Levy or Poisson 
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dP°(t) = <|(3^.-a,,P^)— + Ae 



CirQdB'{t), 



§ij pi pj 



po (pO) 



>dt + 



(133b) 



where 

3^. := 3^ , 

and, depending on the discretization rule 
A* = 1 , Ao = 



-'ir • r; 



A. = -1. 



The A0-term in Eq. (]133bp reflects the modification of differential calculus for 
Ito and backward-Ito SDEs, cf. Appendix El 



Finally, we still note that the proper-time process associated with Eqs. 
is defined by 



dr(t) = (M/P°) dt. 



(134) 



In Section we shall discuss how the lab-time Langevin equations fll32p can 
be reparameterized in terms of r. 



4-1.2 Fokker-Planck equations 



In the case of the Ito rule "0 = *" in Eq. (]132bp . the corresponding Fokker- 
Planck equation (FPE) is given by [cf. Eq. flA.291) ] 



^ d ^ p"- d \ ^ d 
dt p'^ dx^ I * dp^ 



(135a) 



in the case of the Stratonovich-Fisk mid-point rule "0 = o" by 



' d ^p^ d \ ^ d 
dt p^ dx^ I ° 



I 1^ d ■ \ Id ■ 
-3" + a jp^ - -c r-g^c rj fo + 29^^^ '^'^ 



and in the case of the post-point (backward Ito) rule "0 = •" by 



(135b) 



' d p^ d \ d 
di ^ ]fidx^J •''^dp^ 



-r + a^y - cK^,c\) f. + \Z^{c\.cKf.) 



2dp^ 



(135c) 



with fQ{t,x,p) denoting the phase space PDF of the relativistic Brownian 
particle in E, and p° = (M^ + p^Y^'^, respectively. 
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Deterministic initial data X{0) = Xq and -P(O) = Pg corresponds to the 
localized initial condition 



/0(O,a;,p) = 5(a; -ajo) 5(p-Po)- (135d) 

Due to the linearity of Eqs. (11351) . more general solutions can be obtained by 
integrating the solutions of Eqs. fllSSp over a non-localized initial distribution 
/o(a;o,Po)- However, it should be noted that, within a relativistic framework, 
it is very difficult to determine non-localized lab-isochronous initial data by 
means of experimental measurements. 



4.1.3 Free motions in an isotropic bath and Einstein relations 

Similar to the nonrelativistic case, physical constraints on the coefficient func- 
tions a^j and dj come from symmetry properties of the heat bath and from 
the requirement that Eqs. (11321) must reproduce the correct stationary dis- 
tribution and the correct relaxation behavior. For example, in the absence of 
external force-fields, 5"* = 0, and if the heat bath is stationary, isotropic and 
homogeneous in the lab frame E, the coefficient matrices take the simplified 
diagonal form 



a'j = a 5'j , dj = {20^^^ S'j. (136) 

Under the stated assumptions, the functions a and D depend only on the 
Brownian particles' absolute momentum - or, equivalently, on its relativistic 
energy = (M^ + p'^Y^'^ i.e., a{p) = d;(p°) and D(p) = D{p'^). In this case, 
the relativistic Langevin equations (I132p simplify to 



AX\t) = {P'/P^) dt, 

dP\t) = -aP'dt + {2Dy/^ © dB'{t), 

and the Fokker-Planck equations (11351) take the form 



(137a) 
(137b) 



d d 
7^ + 



Jd 

- + 







pi 


d 


pO 




pi 


d 


pO 


dx^ 



d_ 

dp^ 



dp"^ 
dp^ 



ap'f. + ^iDf, 

apJ. + D^f, 

op^ 



(138a) 
(138b) 
(138c) 



An additional constraint on the functions a{p) = a{p^) and D{p) = D{p^) 
arises from thermostatistical considerations: If the motion of the Brownian 
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particle is restricted to a finite volume V C M'^ and if the heat bath is in 
a thermal equilibrium state at temperature 7 = (/cb/3)~^, then the expected 
stationary so lution f nc(x, p) of Eqs. (11351) is a spatially homogeneous Jiittner 



distribution 186 
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U{x,p) = N exp[-/?(M2 +p2)i/2] V), 



(139) 



with J(£c; V) being the indicator function of the accessible volume V as defined 
in Eq. (1103bp . By inserting Eq. (I139p into the Fokker-Planck equations (11351) 
one finds that, depending on the discretization rule, the functions a{p) = 
d(p°) and D{p) = D{jP) must satisfy the generalized fluctuation-dissipation 



relations [ll|, [l7|, [l8| 



= *: = «(/)/-/?£)(/) + £)'(/), (140a) 
= o: {) = a{jP)p^ - P b{jP) + b'{p^)/2, (140b) 
= .: = a(/)p°-/5 £)(/). (140c) 

where £>'(/) ■= dD(/)/dp°. Equations ffHOj) are also referred to as the rel- 
ativistic Einstein relations. In particular, by comparing Eq. (ll40cp with the 
nonrelativistic Einstein relation fl20l) . we note that the mass has been replaced 
with energy in the relativistic case. 



4-2 One- dimensional examples and mean square displacement 



In this part, we will consider one-dimensional example processes [d = 1). 
Their generalization to higher space dimensions is straightforward but analytic 
calculations become more tedious. We first summarize the SDEs for the energy 
process P^it) and the velocity process V(t) := P{t)/P^{t) in Section 14.2.11 
Subsequently, analytical and numerical results for the asymptotic diffusion 



constants of specific example processes will be discussed [15|, l25 



4-2.1 Discretization rules, energy and velocity equations 

Considering an isotropic thermalized heat bath m d = 1 space dimensions and 
the post-point discretization, Eq. (11371) takes form 



dX{t) = {P/P°) dt, 

dP{t) = -a.{P) Pdt+ [2D(P)]i/2 . ^141g^) 
Cf. the discussion in Section 13.2.21 
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where B{t) is characterized by the Gaussian distribution (l5dl) and a,{p) = 
d,(p°) and D{p) = D{p^). Similar to the nonrelativistic case, one can replace 
Eq. (114 lap by the equivalent Stratonovich-Fisk SDB^^ 



dP{t) = -ao{P)Pdt + [2D(P)]^/2 o dB{t) (141b) 
ao{p):=a,{p)-D\p)/{2p), 

or by the equivalent Ito SDE 



dP{t) = -a,{P)Pdt + [2D(P)]^/2 ^ (141p) 
a^{p) ■.= a,{p) - D'{p)/p, 

where D'{p) := dD{p)/dp. Compared with Eqs. (I141ap and fll41bp . the Ito 
form (I141cp is most convenient for numerical simulations. 

Imposing that the solution of corresponding FPE (I138p be given by a one- 
dimensional Jiittner function of the form (11390 . the friction and noise co- 
efficients must satisfy the Einstein relations (I140p . In terms of the friction 
coefficients ^©(p), the Einstein relations can also be rewritten as 



= a,{p) p° -pD{p) 
O = aoip) p'^-/3D{p) 
= a4p) p°-pD{p) 



D\p) p'/{2p) 
D'{p) p'/p. 



From Eq. (]133bp one obtains for the relativistic energy process P^{t) :- 
p2)i/2 the following SDE 



(142a) 
(142b) 
(142c) 

(M2 + 



dP°(t) 



-do(P°)P° 

2Z)(P°) 



+ Xp 



P)(P0) /M\ 



0" 



po VP"/ 



dt + 



1- 



My 



1/2 



&dB{t), 



(143) 



where A* = 1, Ao = 0, and A, = —1. Furthermore, by defining for P{V) 
MV (1 - 1/2) -1/2 new coefficients 



a{V) := a{P{V)) 



D{V) := D{P{V)) 



and applying the (backward) Ito formula, cf. Appendix |Al to the relativistic 
velocity formula V{t) = P/{M^ + P'^f/'^, one finds the following SDE for the 
velocity process 



See also the corresponding discussion by Hanggi 498l |. and Hanggi and Thomas 
(page 293 of Ref. [77,]). 
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dV{t) 



Vdt + 



2ir 

M2 



.1/2 



'1-V 



2\3 



QdB{t). 



(144) 



As discussed in the next section, this equation can be used to calculate the 
asymptotic mean square displacement. 



4-2.2 Asymptotic mean square displacement 

A primary objective within any Brownian theory is to determine the asymp- 
totic diffusion constant Dqo, corresponding to the plateau values in Fig. El For 
a one-dimensional diffusion process X{t) with velocity V{t), the asymptotic 
diffusion constant !Doo is defined by 



Doo= lim([X(t)-X(0)]n/(2t) 



t—*co 



(145a) 



where the spatial displacement is given by 



X(t)-X{0) = f'ds V(s). 
Jo 



(145b) 



The asymptotic diffusion constant Dqo may be expressed in terms of the ve- 
locity correlation function {V{t)V{s)) by virtue of 



Doo = lim - 

t^<^ 2 dt 

1 d /•* 



[x{t)-x{oy^' 

lim ^^J'ds fds' (y{s)V{s')) 
t^oo 2 dt Jo Jo 

lim / ds {V(t)V(s)). 



(146) 



Assuming that the velocity process V{t) is (approximately) stationary, which 
means that {V{t)V{s)) = {V(t — s)V{0)) holds (at least in good approxima- 
tion), and substituting u = t — s, we recover Kubo's formula 



2)oo= lim / du {V(u)V(0)). 



(147) 



As recently discussed by Lindner |25| , for a one-dimensional Langevin equation 
of the form 



dV{t) = -a,{V) Vdt+[2 b{V)Y/^ • dB{t) 



(148) 
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with symmetric coefficient functions, a,{v) = a,{—v) and b{v) = b{—v), the 
Kubo formula fll47p gives rise to the following integral representation for the 
asymptotic diffusion constant: 



Jo+ dy e^^y^ [/;+ dx e-^(^) 



/o+ dz e-ui-)/h{z) 

Here, f + G [0, oo] represents the upper bound for the velocity range, and 



(149) 



V 



U{v) := / dw fi,{w)/b{w) (150a) 
Jo 

is an effective velocity potential with Ito drift 

/i*(f) := a*(f ) V = a,(f ) v — h'{y). (150b) 

In general, the formula ( 1149^ has to be integrated numerically, but for the first 
two models from Section l4:.2.3l the integrals may also be evaluated analytically. 
The (i-dimensional generalization of Eq. (11491) was recently derived by Angst 
and Franchi 1151. 



4.2.3 Examples 

We discuss three specific one-dimensional relativistic Langevin models whose 
stationary momentum distributions are Jiittner functions 0j(p) oc exp[— /?(p^+ 
M^)"^/^] with heat bath temperature T = {k-Q(3)~^ . In this case, the relativistic 
Einstein relation (I140p implies that only one of the two functions a,{p) and 
D{p) can be chosen arbitrarily. 



Constant noise amplitude As a ffist example, we consider the so-called 
'Relativistic Ornstein-Uhlenbeck process' (ROUP), proposed by Debbasch et 



al. Ill, I12I] and also studied by Zygadlo [20]. The ROUP is defined by the 



choice 



a,(p) = etc M/p°, (151a) 

where etc > is a constant friction parameter. From the relativistic Einstein 
relation (11401) . one then finds 



D{p) = a.{p)pT^ = acMp-' =: D^, (151b) 
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i.e., the ROUP corresponds to the hmit case of constant noise amphtude. The 
associated Langevin equation reads 



dP{t) = -a, — Pdt 

M , 
= -a, — Pdt + 



1/2 




dB{t) 



*dB{t). 



(151c) 
(151d) 



The discretization rule is irrelevant here, because the noise amplitude Dc = acM/ j3 
does not depend on the momentum P for this particular case. However, the 
rules of stochastic calculus have to be specified, if one wishes to write down 
the SDE for the associated velocity process V{t) := P/P^. Defining the useful 
abbreviation 



X:=(3M = M/{kB7), 
one finds 



(151e) 



dV{t) = ~a, 



X 

1/2 



Vdt + 



X 



:i-v 



2\3 



dB{t) 



X 

'^{l-V^fT^dBit) 



Vdt + 



(151f) 



:i5ig) 



Using the Ito form, as required in Eqs. (1150bp . we see that the ROUP (1151dp 
corresponds to 



l_^2^3/2^ O 



2\2 



X 



b(v) = — a-v 

X 



2\3 



(152a) 
(152b) 



with an upper velocity bound = c = 1 in Eq. fll49p . In this case the general 
integral formula fll49p for the asymptotic diffusion constant can be evaluated 
analytically by making use of the identity ( 1153^ 



:-ii 



^Ko{x) = £dv exp| 



X 



(l_^2)-(.+2)/2^ 



(153) 
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where, for v = , 1, 2 . . ., Ky{z) denotes the modified Bessel function of the 
second kind 485l |. Remarkably, one recovers for the ROUP the 'classical' re- 
sult [15i, cf. Eq. dM}, 



2) 



ROUP 



(154) 



for all parameter values (etc, 'J', M). 



Constant friction coefficient in the backward-Ito SDE An alterna- 



tive relativistic Brownian motion (REM) model [17|, ll8| corresponds to the 
special case of a constant friction function a,{p) = in the backward-Ito 
SDE fll4ip . In this case, the relativistic Einstein relation ( 11401) yields the mo- 
mentum dependent noise amplitude 



D{p) = a'fp^ (3-\ (155a) 

The relativistic (backward) Ito Langevin equations of this model read explic- 
itly 



dP{t) = -a^P dt+ (^^^^^^ •dB{t) (155b) 

= -«t (^9^) ^ + * dBit). (155c) 

The corresponding SDEs for the velocity process V{t) = P/P^ are given by 



dV{t) = -a^ 



X 



2at 
X 



1 - \/2)5/2 



1/2 



X 

1/2 

*dB{t) 



Vdt + 

dBit) 

V dt + 



^ (1 - v^f/^ 



X 



(155d) 



(155e) 



With regard to numerical simulations, the Ito form is more convenient, cf. 
Appendix Rl 



Various properties of the REM process (11551) have been analyzed by Fa 24 
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Lindner 25|, Fingerla I 395| . and Angst and Franchi 15j. In particular, from 
Eq. fll55ep we see that this model corresponds to 



^(1 

X 



6(^)= ^2)5/2^ 

X 



(156a) 
(156b) 



with an upper velocity bound ?;_|_ = c = 1 in Eq. (11491) . As for the ROUP, 
the integral formula (I149p can be calculated analytically by making use of the 
identity (11531) and one then finds that 
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RBM 



("tx) 



-1 Mx) 

KiixY 



(157) 



At low temperatures f3 := (^b'^) ^ oo, Eq. (11570 reduces to the well-known 
classical result D'^^ = kBl/i^Ma^), cf. Eq. ( !33|) . In the opposite limit of very 



high temperatures, i.e., for (3M -C 1, one finds a logarithmic dependence [15 



= («tM)-i {-7. + ln(2/x) + ©[(/^M)^]} , (158) 

where 7^ — 0.577216 is the Euler constant. However, it should be kept in mind 
that, due to the increasing importance of particle annihilation/creation at high 
energies, classical non-quantum theories become invalid in the high temper- 
ature limit X = l3M -C 1, and, therefore, the asymptotic expansion (11581) 
appears to be of limited practical use. 

By comparing with the ROUP, we observe that T)^^ < T)^^^ holds true 
for same values of the friction coefficients = a^, cf. Fig. [71 Intuitively, this 
can be explained by the fact that, for the ROUP, the absolute value of the 
friction force is bounded by acM, cf. Eq. f ll51dl) . whereas the friction force is 
unbounded for the RBM model fll55bp . thereby suppressing spatial diffusion 
more strongly in the latter case. 



Constant friction coefficient in the Ito SDE The RBM process de- 
fined by Eq. fll55bl) is characterized by a constant friction coefficient a-^, when 
adopting the post-point discretization rule (•). Another model, referred to as 



Fingerle [395l | discusses a fluctuation theorem for this process; see also Cleuren et 



al. [24l[ | 



Equation (|157l) is an equivalent, more compact representation of Lindner's result 
Eq. (10) in 
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RBM(I) hereafter, is obtained by considering a constant friction coefficient 
in the Ito Langevin equation 



dP(t) 



-a, P dt + 



2a, 
/32 



1/2 



(1 + /?P°) 



*dP(t), 



(159a) 



where the noise amphtude is chosen such that the Einstein relation fll42cl) is 
satisfied. The Ito SDE of the associated velocity process Vif) '■= P/P^ reads 



dV{t) = -a, 



1 _ v^) + -(1 - ^2)3/2 + 3 _ ^2^2 
X X 



Vdt + 



I X 

In this case, we have 



1/2 



*dP(t). (159b) 



a* 



^1 _ ^2) ^ 3 2)3/2 ^ 3 2 



6(^;) = ^ (1 _ ^2)5/2 ^^-1(^_ ^2)3 

X ^ 

From these equations the velocity potential is obtained as 



(160a) 
(160b) 



U{v) = In 



X + 1 



Xi'^-v^) + (1 -t;2)3/2 
yielding for the asymptotic diffusion constant: 



-X[l-{l-v 



2^-l/2^ 



^T""^'^ = [a.K.ixT' £ dv — 



,-X(l-^2)-l/2 



X(l-t;2) + (l-t;2)3/2- 



(161) 



The remaining integral can be evaluated numerically. As illustrated in Fig. [TJ 
the theoretical predictions from Eqs. 01571) . fll54p and fll6ip are in good agree- 
ment with the numerically obtained estimates of the asymptotic diffusion con- 
stant. 



The three model processes (1151d|) . (1155b|) and fll59ap give rise to the same 
stationary momentum PDF 0j(p), but their respective relaxation behavior 
differs strongly. This is illustrated in Fig. [6l which depicts the time evolution 
of the spatial mean square displacement divided by time. 



Dr.= {[x{t)-xm')/m, 



(162) 
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Fig. 6. Time evolution of the spatial mean square displacement 
Dt := ([^(*L- ^(0)]^> /{2t) for the ROUP [Tl^ model (solid line) from Eq. (fTsTdl) . 
the RBM fl?!] model (dashed) from Eq. (fTB5b]) and the RBM(I) model (dotted) 
from Eq. (1159ap at same temperature-mass ratio ■= k^J/iMc^) = 1. The plots 
are based on a simulation with N = 1000 trajectories, initial conditions X{Q) = 0, 



P(0) = for each trajectory, and discretization time step At = 10 a 



for all three models at same temperature 7. The curves in Fig. [6] were calcu- 
lated numerically from Eqs. fllSldl) . fll55c| and (I159ap . respectively, using an 
algorithm similar to those described in l25|, see also Appendix IA.5[ 



The ROUP and the two RBM models considered in this part represent spe- 
cial limit cases of the general Langevin equation fll4ip with arbitrarily chosen 
friction coefficient functions. Nonetheless, they yield useful insights: As evi- 
dent from Fig. [3, at moderate-to-high temperatures the diffusion constant can 
vary significantly for different friction models. For realistic systems, the exact 
functional shape (i.e., energy dependence) of the friction function a is deter- 
mined by the microscopic interactions. This result implies that simultaneous 
measurements of the temperature and the diffusion constants can reveal infor- 
mation about the underlying microscopic forces. Below, in Section 14.51 we will 
outline a general procedure for deducing more realistic friction coefficients a 
from microscopic models. 



4-3 Proper-time reparameterization 



The relativistic Langevin equations considered thus far are parameterized in 
terms of the time-coordinate t of the lab frame S. We shall now discuss how 
these equations can be re-parameterized in terms of the proper-time r. For 
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0.01 0.1 1 10 



Fig. 7. Temperature dependence of the asymptotic diffusion constant Dqo for the 
ROUP HH from Eq. (fT^Tdj) . the RBM model [17] from Eq. (fT55bl) . and the RBM(I) 
model from Eq. (|159a|) . Symbols 'o', '+' and 'x' represent the results of computer 
simulations [N = 100 trajectories, initial conditions X(0) = 0, -P(O) = for each 
trajectory, discretization time step At = 10~^ in units of a~J^i^^ obtained by av- 
eraging over the numerically determined values {Ifioo! ^iiO) ^^120> • • • > 2)500}. Solid, 
dashed and dotted lines correspond to the theoretical predictions from Eqs. (jl54p . 
(|157|) and p61|) . respectively. 

this purpose, we start from the Ito Langevin equatioij^ 



dX"(t) = (P°/P°) dt, 
dP*(t) = A* dt + C*j *dP^(t), 



(163a) 
(163b) 



In general, the fun ction s A} and may depend on (t, X*, P*). Upon applying 
Ito's formula [tI, ^ the mass-shell condition P°(t) = (M^ + P'^fl'^^ 
Eq. (]163bl) yields the following equation for the relativistic energy 



dP°(t) 




po 



pipj 
(P0)3 



>dt + |^Q,*dP^(t), (163c) 



where Ai 



A\ D 



J2rCl.cl and 



cV- The Fokker- 



Planck equation for the associated phase space PDF f{t,x,p) reads 



' d d 
dt p° (9x* 



1/ 



_d_ 



(164) 



'^^ For simplicity, we assume that B(t) = {B^ (t)) is d-dimensional, implying that C*^- 
is a square matrix. However, all results still hold if B(t) has a different dimension. 
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As before, we consider deterministic initial data X{0) = Xq and -P(O) = Pq, 
which translates into /(O, x,p) = S{x — Xq) 6{p — p^). 



We are interested in rewriting Eqs. (11631) in terms of the proper-time r, defined 
by 



dr{t) = (1 - VY^Mt , r(0) = 0. 

The proper-time differential may also be expressed as 



(165) 



dr(t) = (M/P°) dt. 



(166a) 



The inverse of the function r(t) will be denoted by X°(r) = t(r) and represents 
the time coordinate of the particle in the frame S parameterized by the proper 
time T. Our goal is to find SDEs for the reparameterized processes X"(r) : = 
X°(t(r)) and P"(r) = P"(t(r)) in S. Heuristically, the SDEs can be derived 
from the relation 



dB^{t) ~ Vdt = (P7M)^^^yd7 ~ (p7M)'^%dS-'(r 



,1/2 



(166b) 



where B^{t) is a standard Wiener process with time-parameter r. Equa- 
tion (fT66b|) can be justified rigorous ly by a pplying theorems for the time- 
change of (local) martingale processes |30l . l499t |. Inserting Eqs. fll66p in Eqs. fll32p 
one finds 



dX"(r) = (P°/M) dr, 
dP*(r) =A' dT + &j * dP^ (r), 

where the transformed coefficients are given by 



(167a) 
(167b) 



A':={P"/M) A\X^,X,P), 
c\.:=(P7M)1/2^\.(X°,X,P) 



(167c) 
(167d) 



Analogous to Eq. fll63cp . the reparameterized energy equation is obtained as 



dP°(r) = \A~ + ^ 
^ ' 1 po 2 



§ij pi pj 



po (p 



o^3 



^dt + Ci^*dP^(r), (167e) 



where now D^^ := c*c^. The FPE for the associated probability density /(r, x^, x, 
reads 
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' d p°' d \ - d 



- - 1 (9 / - 
-A'f + -—{D'^f 



(168) 



We note that now / dx^d'^'xd'^p gives probability of finding the particle at 
proper-time r in the {2d + l)-dimensional interval [t, t + dt] x [x, x + dx] x 
[p,p + dp] in the inert ial frame E. 

An interesting consequence of this reparameterization can be illustrated by 
considering the free motion in a thermalized heat bath, which is stationary, 
isotropic and position independent in the lab frame S. As discussed above, in 
this case a plausible ansatz reads 



A' 



-a{p^)p' 



31 



(169a) 



where the friction and noise coefficients a and D depend on the energy p° 
(M^ +p^)^/^ and satisfy the relativistic Einstein relation 



= «(/)/ + D'(/) -/3Z}(p"). 



(169b) 



Then, the stationary mome ntum dis tribution of the lab-time FPE ( 11641) is 



given by a Jiittner function [l86l . |202 



I.e., 



lim / oc 



t— ►oo 



exp(— /3p°) 



0J(P) 



(170a) 



Remarkably, however, in this case the stationary solution f^o of the correspond- 
ing proper-time FPE (11681) is given by the modified Jiittner function 474 



/oo := lim / oc 



exp(— /5p°) 



(170b) 



This can be confirmed by direct numerical simulation of Eqs. (11631) . An exam- 
ple is given in Fig. [HI which depicts the numerically obtained PDF of the abso- 
lute momentum \P\ for the one-dimensional ROUP ll|]. The numerical PDFs 
were calculated from 10000 sample trajectories, by measuring the momentum 
at constant lab-time t (x) and constant proper-time r (o), respectively. The 
solid and dashed lines show the corresponding standard and modified Jiittner 
functions, respectively. The physical explanation for the difference between 
and /oo is that measurements at t = const and r = const are non-equivalent 
even if both r, t — oo. 

The modified Jiittner distribution 0m j, representing the stationary solution of 
the proper-time FPE, differs from the standard Jiittner function 0j (1118bp by 
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Fig. 8. 'Stationary' probability density function (PDF) of the absolute momen- 
tum |P| measured at time t = 15 (x) and proper-time r = 15 (o), respectively, 
from 10000 sample trajectories of the one-dimensional {d = 1) relativistic Orn- 
stein-Uhlenbeck process [HJ, corresponding to coefficients D{p^) = Dc = const 
and a{p^) = fiD^jP in Eqs. (jl63p and (|169p . Simulation parameters: dt = 0.001, 
M = c = /3 = Dc = 1. The solid and dashed lines show the corresponding stan- 
dard and modified Jiittner functions cpj oc exp(— and (pMj oc ex.p{—(3p^)/p^, 
respectively. 



a prefactor proportional to the inverse energy. As discussed in Ref. [47J], the 
modified distribution 0mj can be derived from a relative entropy principle by 
using a Lorentz invariant reference measure in momentum space, whereas the 
Jiittner function 0j is obtained when adopting a constant reference density 
(Lebesgue measure), cf. Eq. fll23al) . 

We may summarize: If a relativistic Langevin-Ito process has been specified 
in the inertial frame S and is parameterized by the associated S-coordinate 
time t, then this process can be reparameterized by its proper-time r and the 
resulting process is again of the Langevin-Ito type but exhib its a modified 
stationary distribution. With regard to applications 3161 . l319l | the latter fact 
means that the correct form of the fluctuation-dissipation relation depends on 
the choice of the time-parameter in the relativistic Langevin equation. In prin- 
ciple, however, one can - similar to the case of purely deterministic relativistic 
equations of motions - choose freely between different time parameterizations 
by taking into account that the noise part needs to be transformed differently 
than the deterministic part. 



4-4 Moving observers 



The preceding discussion has focused on Langevin equations that describe the 
random motion of a relativistic Brownian particle in the lab frame, defined 
as the rest frame of the heat bath. In this part we would like to address 
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the following question [26 



30| : Assuming that a Langevin equation of the 
type (11411) holds in the lab frame S - how does the corresponding process 
appear to a moving observer? 

There exist (at least) two different ways to tackle this problem: Either one 
uses the Lorentz invariance of the phase space density, or one applies a Lorentz 
transformation directly to the Langevin equation 26, 3ol . 



4-4- 1 Lorentz transformation of the phase space density 



An observer at rest in the lab frame would measure the PDF / governed by 
the Fokker-Planck equations (I135p or (11381) . As discussed in Sec tion [3.1.21 the 



one-particle phase space density transforms as a Lorentz scalar [151 



I.e. 



f\t',x',p') = f{t{t',x'),x{t',x'),p{p')), 
and, conversely. 



:i71a) 



fit,x,p) = f'it'{t,x),x'{t,x),p'{p)), (171b) 
where {t',x',p') and {t,x,p) are related by the Lorentz transformation 



x'^(t,a;) = AV + AV, (172a) 
p"{p)=A\{m^ + p2)i/2 + _ (172b) 

For an observer moving at constant velocity w through the lab frame, the 
corresponding Lorentz boost matrix elements are given by Eq. fl83l) . Therefore, 
in order to find /', it suffices to solve the Fokker-Planck equations (11351) in the 
lab frame S, e.g., for a given t-simultaneous initial condition f{0,x,p), and 
to insert the solution into (llTlap . 



4.4-2 Lorentz transformation of the Langevin equation 

Alternatively, in order to obtain explicit SDEs for the Brownian motion in S', 
one has to apply a Lorentz transformation to the Langevin equations in the 
lab frame E [261 . 3o| . This can be done by employing a similar heuristics as in 



the case of the proper-time reparameterization, cf. Section 14. 3[ 

To illustrate this in more detail, we start again from the Langevin equa- 
tions (11631) and consider a Lorentz transformation Q from the lab frame S to 
another inertial frame S', mediated by a constant matrix A"^ that leaves the 
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metric tensor rjajs invariant. For convenience, we exclude time reversal trans- 
formations, i.e., we restrict ourselves to proper Lorentz transformations with 
A°o > 0. It is convenient to proceed in two steps: First we define 

¥""{1) := A%X^(t) , G"'{t) := A%P^{t). 

Then we replace t by the coordinate time t' of S' to obtain processes 

X'"(t') = , = G"°(t(t')). 

In order to find the SDKs for and P'"{t'), we note that dY'%t) = A%dX^'{t) 

and, therefore. 



dt'(t) = dY'\t) = \ dt = —dt= ' \ dt, (173a) 

where is the inverse Lorentz transformation. A similar heuristics as in 
Eq. (]166bp then gives 



dB^{t) 



■ pO \ 1/2 



^dt' 



(A- 



-1\0 pl^^ 

) fl-^ 



p/0 



1/2 



*dB'^{t'), 



;i73b) 



where B'^{t') is a Wiener process with time parameter t'. Defining 



2 



ptpj 



(174a) 
(174b) 



we can introduce Lorentz-transformed coefficient functions in S' by 



A'^(x'°,a;',p') 



c'Ux",x',p') 



■(A-^)V^" 



p. 



/o 



X 



A\A^ ((A-^)V^ (A-^)Va:'^ (A-^)^^' 



(175a) 



p. 



/o 



A^ ((A-i) (A-') V^'", (A-^) Vp"^) . (175b) 

Substituting the formulas fll73p - (ll75l) into the Langevin equations fll63p . one 
finds that the particle's trajectory (X'(t'), P'{t')) in E' is governed by an SDE 
of the form (ll|, [l3, [3 
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dX'°(t') = (P'"/P'°) dt', (176a) 
dP'\t') = A" dt' + c"j * dB'^it'). (176b) 

It should be noted that, due to the [ ■ ]-prefactors in Eqs. (11751) . the coefficients 
A'^ and c^^ do not form Lorentz tensors. This is essentially a consequence of 
the coordinate-time parameterization in Eqs. (I163P and (I176P Furthermore, 
it is straightforward now to obtain the corresponding SDE for the energy 
P'° = (M^ + p'^)i/2 in as well as the Fokker-Planck equation by replacing 
unprimed with primed quantities in Eqs. (I163cl) and (11641) . respectively!^ 



4-5 Relativistic binary collision model 



The preceding subsections have focused on general aspects of relativistic Langevin 
and Fokker-Planck equations. Similar to the nonrelativistic case, relativistic 
SDEs present a useful tool for analytical and numerical studies of relaxation 
processes in relativistic systems. Stochastic models of this type provide a 
simplified picture of the underlying microscopic dynamics. In order for the 
Langevin approach to be successful, one must know in advance which fric- 
tion coefficient function a{P) and noise amplitude D{P) are appropriate for 
the system under consideration. Realistic friction models can be obtained, 
e.g., by d erivi ng Fokker-Planck equations from relativistic Boltzmann equa- 
tions 16l. |299| . In this context, howev er, it should be noted that the validity 
of the relativistic Boltzmann equation 312l | is less understood than that of its 
nonrelativistic counterpart. In the remainder of this section, we will sketch an 
alternative procedure for obtaining friction coefficients and noise amplitudes 



from a simple microscopic interaction model j21j. The latter can be viewed as 



the direct relativistic generalization of the elastic binary collision model from 
Section [2221 More precisely, we consider a one-dimensional system consisting 
of a heavy Brownian particle (mass M) which is embedded into a heat bath 
of smaller particles (mass m <C M, total number ^ 1). Our model assumes 
that the stochastic motion of a Brownian particle arises due to frequent elastic 
interactions with the surrounding heat bath particles. Similar to Section [5221 
we are interested in finding the 'best' approximation of the 'exact' dynamics 
within the class of SDEs defined by Eq. (I14ip . 



The coefficients A* of the corresponding r-parameterized process are the spatial 
components of a Lorentz four- vector A^; similarly, the quantities = c^^c^^ 
are spatial components of a second rank Lorentz tensor. 

''^ In particular, if the coefficients of the lab frame Langevin equation satisfy the 
relativistic Einstein relation ()169bp . then the stationary momentum distribution 
in the moving frame S' is given by a boosted Jiittner distribution of the form 
4>'{p') oc exp(/3?7^p'"), where is the four-velocity of the heat bath in S'. 
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Relativistic collision kinematics To begin with, we consider a single 
collision of the Brownian particle (momentum P, energy E = pO) with a heat 
bath particle (momentum p, energy e). The relativistic energy, momentum and 
velocity of the two particles before the collision are given by 

P = MV-f{V), E{P) = [m^ + P^y^\ (177a) 

p = mv'y{v), e{p) = (m'^ + p"^^ ^ . (177b) 

where 7(f) := (1 — v'^)~^^'^. Considering elastic interactions, the collision kine- 
matics is governed by the relativistic mass-energy-momentum conservation 
laws 



M = M, m = m, E + e = E + e, P + p = P + p, (178) 

where hat-symbols refer to the state after the collision. Inserting Eqs. (11771) 
into the conservation laws fll78p . and sol ving for the momentum of Brownian 



particle after the collision, P, we obtain [21 



P = -f{u)^[2uE - {l + u^)P], (179a) 
where the collision-invariant center-of-mass velocity u is given by 



u{p,P) = ^. (179b) 

Accordingly, the momentum change AP^ := P — P of the Brownian particle 
in a single collision with the heat bath particle 'r' is given by 



APr = -2j{Urf P + 2j{Ur)' Pr, (180) 

where Ur '■= u{pr,P) and e,. := e{pr). In the non-relativistic limit case, where 
< 1, E ~ M and ~ m, Eq. fflSOjl reduces to Eq. fliSj) . 



Furthermore, by making the same assumptions as in Section 12.2.21 we find 
that the momentum change SP{t) := P{t + 6t) — P{t) of the Brownian particle 
during a small-but-sufficiently-long time interval [t, t+6t] can be approximated 
by 
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dP{t) ^-2Y: liUrY Pit) Irit, 6t) + 

r=l ''" 

^ E 

2 E 7(Wr.)' TT— Pr Ir{t, St). (181a) 
r=l ^ + er 

Formally, the collision indicator Ir(t,6t) is again given by [cf. Eq. f l53al) 

Ir{t,5t)^^\Vr-V\5{Xr-X), (I81b) 

but now one has to use the relativistic velocities V = P/ (M^ + P^)-*^/^ and Vr = 
Pr/iM"^ +Pr)^^'^, respectively. Equation fll81al) is the relativistic counterpart 
of Eq. (]53bl) . Heuristically, the first term on the rhs. of Eq. fll81ap can again 
be interpreted as 'friction', while the second contribution may be viewed as 
'noise'. 

Bath distribution and drift Similar to the nonrelativistic case, Eqs. (11811) 
can be used to calculate the statistical moments of the momentum increments 
- provided one specifies the phase space distribution of the heat bath particles. 
We will assume here that the heat bath is in a thermal equilibrium state, so 
that the one-particle phase space PDF is given by a spatially homogeneous 
Jiittner function on V = [0, L], i.e.. 



fl{Xr,Pr) = {ZiLY 



exp 



-P{m'+pl)'/' 



6(L — Xr)Q{Xr 



[182) 



where T = (/^/cb)"^ is the temperature, and Zi = 2m Ki{l3m) with Ki{z) de- 
noting the modified Bessel function. With regard to our subsequent discussion, 
we are interested in calculating the mean drift force 3C, defined bj- 
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P{t)=p 



(183) 



Inserting SP{t) from Eq. (118 lal) . we find 



X{p) = -nb ( 7(Mr) 



nh ( l{Ur)" 



2 '^r 



E + er 
E 



L \Vr - V\ 6{Xr - X) 



E + e, 



Pr L \Vr — V\ 6{Xr — X) 



P{t)=p) P + 



Pit) = p 



:i84) 



In principle, higher moments can be calculated in a similar manner, but then one 
has to specify the corresponding many-particle heat bath PDFs. 
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Fig. 9. Relativistic binary collision model. Mean drift force ^{P) '■= {SP{t)/5t)^ 
numerically evaluated from Eq. (|184p for different values kB^, with nb = N/L 
denoting the number density of the heat bath particles. 

where nb = N/L is the number density of the bath particles. In order to 
determine 3C(p), we note that for some arbitrary function G{p, P), we have 



G{pr,P)L\Vr-V\S{Xr-X) P{t) = p^ 

f dpr G{pr,P) exp[-/5(m2+p2)i/2 

Pr P 



X 



(m2 + p2)l/2 (]Vf2 + p2)l/2 

The rhs. of Eq. (11841) involves the functions 



(185) 



e{pr 



Pr 



(186a) 
(186b) 



Unfortunately, it is very difficult or perhaps even impossible to analytically 
evaluate the integral (11851) for the functions Gi/2- Figure [9] depicts the mean 
drift force %, obtained by numerically integrating the formula (11851) for dif- 
ferent values of P. 



Langevin approximation We conclude this section by discussing how one 
could, in principle, approximate Eqs. (I18ip by, e.g., a nonlinear backward Ito 
SDK of the form (I14ip . reading 
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dpit) = 



-a{P) Pdt + [2D(P)]i/2 . 



(187) 



The considerations for the two-component Jiittner gas from Section [3] imply 
that the stationary momentum distribution of the Brownian particle in the 
binary collision model is given by the Jiittner function 



Mp) = exp[-/?(p2 + m2)1/2]/[2M K,{(3M)]. (188) 

Hence, in order for Eq. f ll87p to yield the expected Jiittner distribution, the 
functions a and D must be coupled by the relativistic Einstein relation f ll40p . 



Dip) = P~'aip)Eip), (189) 

where E = {p'^ + M'^y/'^. In order to determine the function a, we demand 
that the Langevin equation yields the same mean drift force % as the collision 
model, i.e.. 



\ dt 



P{t) = p 



l5P{t) 



6t 



P{t) = p 



(190) 



For the post-point (backward-Ito) Langevin equation (I187p we know that [cf. 
Eq. (lAjsj) ] 



/dm 

\ dt 



P{t)=pj = -a{p)p + —D{p). 



(191) 



Thus, by means of the Einstein relation fll89p . the Ihs. of Eq. fll90p is given 
by 



/dP(t) 



dt 



P{t)=p\ = -a{p)p + p-'-^[a{p)E{p)], 



(192) 



and the condition fll90p becomes equivalent to the differential equation [cf. 
Eq. dSSD] 



-aip)p + p-'-^[a{p)E{p)] = X{p). (193) 

In the case of the collision model, where the function %{p) is not exactly 
known, one could, e.g., try to fit %{p) by a simple analytic expression and, 
subsequently, use this approximation in Eq. (11930 . 
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5 Non-Markovian diffusion processes in Minkowski spacetime 



The preceding section has focused on relativistic Brownian motions in phase 
space. In the remainder we will discuss relativistic diffusion models in Minkowski 
spacetime, i.e., continuous relativistic stochastic processes that do not explic- 
itly depend on the momentum coordinate. On the one hand, such spacetime 
processes may be constructed, e.g., from a Brownian motion process in phase 
space by integrating out the momentum coordinates. As a result of this av- 
eraging procedure, the reduced process for the position coordinate will be 
non-Markovian. Alternatively, one can try to derive or postulate a relativis- 
tic diffusion equation and /or di ffusi o n pr opagators in spacetime on the basis 



of microscopic models |10l . l28l . 1 1531 . l366l | or plausibility considerations 23 



Regardless of the approach adopted, in order to comply with the principles 
of special relativity, the resulting spacetime process must b e no n-Markovian, 
in accordance wit h th e results of Dudley (Theorem 11.3 in 333| ) and Hakim 
(Proposition 2 in 338[ |). Roughly speaking, this means that any relativistically 
acceptable generalization of the classical diffusion equation ([1]) should be of 
at least second order in the time coordinate. 



The construction and analysis of relativistic diffusion models in Minkowski 
spacetime poses an interesting problem in its own right. Additionally, the 
investigation of these processes becom es re l evan t in view of potential analo- 
gies with relativistic quantum theory 



383, 391 



similar to the analogy be- 
tween Schro dinge r 's eq uation and the diffusion equation ([1]) in the nonrela- 



tivistic case 



500, 



501 



The present section intends to provide an overview 
over classical relativistic d iffus i on models that have been d iscussed in the lit- 
erature 0, M, M, S, H H, H, H, H, H, H, [Mil • For this purpose, 
we first recall basic properties of the Wiener (Gaussian) process, which con- 
stitutes the standard paradigm for nonrelativistic diffusions in position space 
(Section 15. ip . Subsequently, relativistic generalizations of the nonrelativistic 
diffusion equation fll and/or the nonrelativistic Gaussian diffusion propagator 
will be discussed j23l] . 



5.1 Reminder: nonrelativistic diffusion equation 



We start by briefly summarizi ng a few r elevant fa cts about the standard non 



relativistic diffusion equation [2871 . l339l . |422| . |501 



t > t, 



0) 



(194) 



87 



where D > denotes the spatial diffusion constant, and g{t, x) > the one- 
particle PDF for the particle positions as G M"^ at time t. Within classical 
diffusion theory, Eq. (11941) is postulated to describe the (overdamped) random 
motion of a representative particle in a fluctuating environment (heat bath). 
In particular, Eq. fll94p refers to the rest frame of the bath. 

There exist several well-known ways to motivate or derive the phenomeno- 
logical di f fusion eq uation (11941) by means of microscopic models (see, e.g.. 



2871 . l339l . I422I . ISOli p. With regard to our subsequent discussion of relativistic 



alternatives, it is useful to briefly consider a 'hydrodynamic' derivation [502 
which starts from the continuity equation 



^^g(t,x) = -V-jit,x), (195) 

where j{t,x) denotes the current density vector. In order obtain a closed 
equation for the density g, the current j has to be expressed in terms of g. 
One way of doing this is to postulate the following rather general ansatz {cf. 



Eq.(2.81) in [50^} 



j{t, x) = -V /* dt' K{t - t') g{t', x), (196) 



where, in general, K may be a memory kernel. However, considering for the 
moment the memory-less kernel functiorf^ 



KF{t-t') -.= 21) 6{t-t'), (197) 
one finds 

j^{t,x) = -'DVgit,x). (198) 

Upon inserting this expression into the continuity equation (11951) . we recover 
the classical diffusion equation fll94p . 

Now, it has been well-known for a long time that the diffusion equation fll94p is 
in conflict with the postulates of special relativity. To briefly illustrate this, we 
specialize to simplest case of ci = 1 space dimensions, where = d"^ /dx^. In 
this case, the propagator of Eq. (I194p at times t > to is given by the Gaussian 



The factor '2' in Eq. (jl97p appears because of the convention 
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p{t,x\to,xo) 



1/2 



exp 



"4D(t-to) 



(199) 



The propagator (11991) represents the solution of Eq. (11941) for the initial con- 
dition 

g(to, x) = 6{x — Xq). 

That is, if X{t) denotes the random path of a particle with fixed initial position 
X(to) = Xq, then p{t, x\to, Xo)dx gives the probability that the particle is found 
in the infinitesimal volume element [x, x + dx] at time t > to- As evident from 
Eq. (11991) . for each t > there is a small, but non- vanishing probability that 
the particle may be observed at distances |a;— Xq] > c(t—to), where c = 1 is the 
speed of light in natural units. The evolution of the nonrelativistic Gaussian 
PDF from Eq. ffT99D is depicted in Fig. M 



It is worthwhile to summarize a few essential properties of Eqs. (I194p and 
(11991) : Equation (11941) is a linear parabolic partial differential equation. Due to 
the linearity, more general solutions may be constructed by superpositioning, 
i.e., by integrating the solution (11991) over some given initial PDF go{xo)- 
Equation (I194p describes a Markov process which means that the transition 
PDF (I199p satisfies the Chapman-Kolmogoroff criterion 



p{t,x\to,xo) = / dxi p{t,x\ti,xi) p{ti,xi\to,xo) 



(200) 



for all ti G {to,t). The corresponding diffusion process X{t) can be character- 
ized in terms of the following SDE: 



(201) 
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Fig. 10. Spreading of the Gaussian PDF Q{t,x) = p(t,x|0,0) from Eq. p99p at 
different times t, where t is measured in units of Djc? . At initial time t = to = 0, 
the PDF corresponds to a 5-function centered at the origin. 
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where B{t) is a standard Wiener process as defined in Section [2.1.2[ Formally, 
Eq. (12011) may be obtained from the Langevin equations (!22|) of the classical 
Ornstein-Uhlenbeck process with 3" = as follows: First we rewrite Eq. fl22bp 



as 



dV{t) 
Ma 



-Vdt 



2Do_ 



1/2 



(202) 



Upon letting (Ma) —>■ oo and Dq ^ oo such that T) = Do/laMY remains 
constant, the Ihs. of Eq. (12021) should become negligible. Then, by making use 



of dX = Vdt, Eq. (I20ip is recovered] I This limiting procedure defines the 



so-called overdamped regime of the Ornstein-Uhlenbeck process. The mean 
square displacement of the overdamped process (12011) is given by |287 



[X{t)-X{to)]'' 



dx {x — xo)^ p{t, x\to, Xq) 
2D (t-to), 



(203) 



qualitatively similar to the asymptotic behavior of the classical Ornstein- 
Uhlenbeck process; cf. Eq. ( I25l) . Finally, we note that the solution of Eq. (11941) 
with initial condition 



Q{to,x) = 0o{x), 



(204) 



can be expressed in terms of the Feynman-Kac formula [T^, |76 



g{t,x) = (go{x + {2Vy/^B{t) 



(205) 



where ( ■ ) indicates an average with respect to the Wiener measure of the 
standard Wiener process B(t) with initial condition B(to) = 0. Equation (12051) 
yields an efficient Monte-Carlo simulation scheme for computing the solutions 
of the diffusion equation (I194p for a broad class of initial distributions qq. 



5.2 Telegraph equation 



The problem of constructing continuous diffusion models which, in contrast 
to the classical nonrelativistic equations (11941) and (11991) . avoid super luminal 
velociti es, has attracted considerable inte r est o ver the past years |lOl. l23l. |27. 
28|, 0, H H, H, H, H, H, H, |39l|. Nonetheless, it seems fair to 



Debbasch and Rivet 1^ discuss the difficulties that arise when attempting a 
similar reduction for the relativistic Ornstein-Uhlenbeck process. 
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say that a commonly accepted solution is still outstanding. Apart from the 
theoretical challenge of developing a consistent relativistic diffusion theory, 
there exist several practical applicatio ns in c luding, e .g., the analysis of data 
from high energy colli sion experime nts 320 . 321 ^ 322| or the diffusion of light 
through turbid media |372l . l503l . |504| | and foams I374l375j . In this context, 
a frequently considered alternative to the cl assical diffusion equation (|194D is 
give n by the telegraph equation |10l. |27. 
380l 



29|, 13201, 13211, 13661, 1370, \S77\, 13781, 1379| 



Here, T) > plays again the role of a diffusion constant, while > is 
an additional relaxation time parameter. Equation fl206p can also be obtained 
from the continuity equation (11951) . provided one uses the exponential memory 
kernel (Section 2.10 in Ref. |502] ) 



KT{t - t') := - exp[-(t - t')/r.], (207) 

instead of Kp from Eq. (I197p . Similar to Eq. fll94p . the telegraph equa- 
tion fl206p refers to a special frame where the background medium, causing 
the random motion of the diffusing test particle, is at rest (on average). The 
'nonrelativistic limit' corresponds to letting r^, ^ in Eq. ( I206p . which leads 
back to Eq. fll94p . For > 0, Eq. fl206p is a linear hyperbolic partial differ- 
ential equation. Because of the second order time derivative in Eq. fl206p . one 
now also has to specify the first order time derivative of the initial distribution 
at time to- 

Considering particular initial conditions 



d 

g{to,x) = S{x - xo), ■^Qito,x)= 0, (208) 



one finds that the corresponding solution of Eq. (12060 is given by lOl, 1366 



g-(t-to)/(2r„) r 

p{t, x\to, Xo) = < S[\x -xo\-v{t- to)] + 



2r„t> 



2n, e 



(209a) 



Here, we have abbreviated 
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X — Xq 



1/2 



V := {D/r, f\ 



(209b) 



and the modified Bessel functions of tlie first kind, Iu{z), are defined by 



2k+u 



with T{z) denoting the Euler gamma function. According to our knowledge, 
the solution fl209p was first obtained by Goldstein in 1938/1939. Actually, 
Goldstein derived the result (1209 ) by considering the continuum limit of a 
persistent random walk model 3691 ]: subsequently, he proved that this function 



satisfies the telegraph equation (12061) . cf. Section 8 of his paper [366 
The propagator (I209p is characterized by two salient features: 



(1) As evident from the 5-function term, the solution exhibits two singular 
diffusion fronts traveling at absolute velocity v := {Ti/t^Y^'^ to the left 
and right, respectively. 

(2) Due to the appearance of the Heaviside G-function, the solution is non- 
zero only within the region |x — Xo| < '^(t — to); i-e., upon fixing r„ such 
that V = c = 1 the solution vanishes outside the light cone. 



Thus, in contrast to the nonrelativistic propagator (I199p . Eqs. (12090 define a 
relativistically acceptable diffusion model. Because of the second order time 
derivative, the telegraph equation (I206P describes a non- Mar kovian proce ss, in 
accordance with the aforementioned theorems of Dudley 333[ | and Hakim 338 . 
The non-Markovian character of the propagator (12091) can also be proven di- 



rectly by verifying that this solution does not fulfill the condition 



The linearity of Eq. (I206p implies that more general solutions can be ob- 
tained by integrating the propagator (12091) over some given initial distri- 
bution ^o(a;o)- In principle, one may also construct other classes of solu- 
tions wi th dp(to, x)/dt ^ 0, e.g., by applying a Laplace- Fourier transforma- 
tion [l^, ISOS , 506 ] to Eq. (12060 . We note, however, that in order for the solu- 
tion g{t,x) to remain normalized and positive at all times t > to, additional 
constraints on the initial conditions must be imposed, cf. Eq. (12081) . Various 



solutions and extensions of the telegraph equation (12060 . including different 
types of boundary conditi ons, additional extern al sources , etc., ha ve b een dis- 
cussed, e.g., by G oldstein |366| . Masoliver et al. |505l. l506l| . Foong |507|, Foong 



and Kanno [508|, Renardy [509j, and Dorogovtsev [510 



Similar to the nonrelativistic diffusion equation (I194p . the telegraph equa- 
tion (12060 may be derived and/or motivated in many different ways. A detailed 
overview is given by Masoliver and Weiss [lO| , who discuss four different possi- 
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bilities of deducing Eq. (I206P from underlying models; see also Koide [27|, 1291 . 
During the past decades, the telegraph equation fl206p has been used to de- 
scribe a number of different phenomena. The applications include: 



Transmission of electrical signals. According to Masoliver and Weiss (I206p . 
the earliest deriva tion of the telegraph equation is based on a paper by 
William Thomson 365| (who later became Lord Kelvin), published in 1855. 
He considered the problem of how to transmit electrical signals without 
distortion, a question closely related to the design of the first transatlantic 
cable. 

Continuum models of persistent diffusion. It seems that the c oncept of 
persistent diffusion was first proposed in 1917/1922 by Fiirth 3681 . 1369 
who aimed at describing the random motion of biolo gical objects. Indepen- 
dently, a similar approac h wa s suggested by Taylor [370 1 in an attempt to 
treat turbulent diffusion [363]. Fiirth and Taylor considered discrete mod- 
els, assuming that a particle moves with constant absolute velocity between 
neighboring lattice points. At each lattice point, the particle is either back- 
scattered or transmitted, with the transmission probability being larger 
than the back-scattering probability {persistence). A few decades later, in 



1950, Goldstein |366| demonstrated for the ID case I that, for a suitable 
choice of the transition probabilities, the continuum limit of this model 
leads to the telegraph equation. Hence, in contrast to the ordinary diffusion 
equation fll94p . the telegraph equation (12060 relies on asymmetric transi- 
tion probabilities, causing the non-vanishing probability concentration at 
the diffusion fronts. In more recent years, persistent diffusion models have 
been e mployed to d e scribe th e propagation of photons in thin slabs and 



foams 371, 372, 373, 374, 375 



Heat transport and propagation of heat waves. In this case, the function 
^(t, x) in Eq. (12061) is interpreted as a temperature field and the normaliza- 



tion condition is dropped. For a detailed account of the vast ht erature on 



heat waves we refer to the review articles of Joseph and Preziosi [3771 . |378 



High energy ion collision experiments. In recent years, the telegraph equa- 
tion has been used to estimate the dissipation of net charge fluctuations, 
which may obh terate signals of QCD phase transitions in nuclear colli- 
In this context, however, the coordinate x in Eq. 



sions 



32G, 321 



is interpreted as a rapidity variable. 



Anot her i nteresting aspect of the telegraph equation is elucidated in a paper by 
Kac |367l | . He observed that the solutions of Eq. (I206p with initial conditions 



d 

Q{to,x) = Qo{x), -Q^Q{to,x) = 0, (210) 



76 



Boguna et al. ^] discuss persistent random walks in higher space dimensions. 
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may be expressed in the fornj^ 



+ 



(211) 



where v = (D/r^)^/^, and (■) indicates an average with respect to the t^- 
parameterized Poisson process N{t); i.e., for any given time t > to have 



Prob{A^(t) = A;} 



-(t-to)/(2r„) _ + N A: 



and for any finite sequence to < ti < . . . < tn the increments 



A; = 0,1,2,...; (212a) 



iV(ti) - N{to), N{t2) - iV(ti), . . . , iV(t„) - iV(t„_i) 



(212b) 



are independent. Equation (121 ip is the direct counterpart of Feynman-Kac 
formula (12051) for the classical diffusion equation. Equations (12111) and (12121) 
together provide an efficient Monte-Carlo simulation scheme for computing 
solutions of the telegraph equation (12061) . Moreover, the Poisson path inte- 
gral representation (12 lip di sclose s an interesting correspondence between the 
free-particle Dirac equation 382| ar id th e telegraph equation (12060 . which was 



first pointed out by Gaveau et al. |383| in 1984: The solutions of both equa- 



tions may be linked by means of an analytic continuation. This connection is 
similar to the relation between the diffusion equatio n (|194l) and the free par- 
ticle Schrodinger equation in the nonrelativistic casell^ The crucial difference 
is given by the fact that the measures of the functional integration refer to 
different underlying processes, respectively. 



However, the telegraph equation 



is not the only possible relativistic 



generalization of the nonrelativistic diffusion equation ([T]) and, recently, ther e 



380, 381 



has been some controversy about its applicability and validity [379 
An early critical discussion of Eq. (12061) in the context of relativistic heat 
transport was given by van Kampen [153| in 1970. Starting from a simple 
microscopic model, consisting of a cloud of material particles that exchange 



'^'^ The result (j21ip may be gene ralized to an arbitrary number of space dimensions; 
cf. pp. 500 in Kac's paper [3671 ]. 

78 Pq]^ further reading about path integral r epresentations of the Dirac pro pagator 
we refer to the papers of Ichinose 384 . 385], Jacobson and Schulman SSd ]. Barut 
and Duru [3871], and Gaveau and Schulman [388l ]: see also f ootn ote 7 in Gaveau et 
al. [3831 ] and problem 2-6, pp. 34-36 in Feynman and Hibbs [3891 ]. 
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electromagnetic radiation, van Kampen derived an integral equation for the 
temperature of the particles as function of time and space. He then showed 
how the telegraph equation fl206p can be recovered as an approximation to 



the more precise integral equation, but that the validity of this approximation 
breaks down in the vicinity of the diffusion fronts. 



Similarly, the singular diffusion fronts predicted by Eq. fl209p represent a source 
of concern if one wishes to adopt the telegraph equation fl2U6p as a model for 
particle transport in a random medium. While these singularit ies may be 



acceptable in the case of photon diffusion [3711 . l372l . l373l . 1374 l375l |. they seem 



unrealistic for massive particles, because such fronts would imply that a finite 
fraction of particles carries a huge amount of kinetic energy (much larger than 
mc^). In view of these shortcomings, it appears re ason able to explore other 



constructions of relativistic diffusion processes [23|, l390l ] . In the next part we 
will discuss a different approach [i^ that may provide a viable alternative to 
the solutions of the telegraph equation. 



5.3 Relativistic diffusion propagator 



In principle, one can distinguish two different routes towards constructing rela- 
tivistic diffusion processes: One can either try to find an acceptable relativistic 
diffusion equation, or one can focus directly on the structure of the diffusion 
propagator. In the present part we shall consider the latter approach [2^. The 
basic idea is to rewrite the nonrelativistic diffusion propagator fll99p in such 
a form that its relativistic generalization follows in a straightforward manner. 
This can be achieved be reexpressing Eq. (I199p in terms of an integral-over- 
actions. 

For this purpose, we consider a nonrelativistic particle traveling from the event 
= (^0, xq) to X = {t, x) and assume that the particle can experience multiple 
scatterings on its way, and that the velocity is approximately constant between 
two successive scattering events. Then the total action (per mass) required 
along the path is given by 

a{x\xo) = 1 f' dt' v{t'f, (213) 
2 Jto 

where the velocity v{t') is a piecewise constant function, satisfying 

x = xo+ f &t' v{t'). (214) 
Clearly, the nonrelativistic action fl213p becomes minimal for the deterministic 
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(direct) path, i.e., if the particle does not coUide at all. In this case, it moves 
with constant velocity v(t') = (x — Xo)/(t — to) for all t' G [to,t], yielding the 
smallest possible action value 



On the other hand, to match the boundary conditions it is merely required that 
the mean velocity equals {x — xq)/ (t — to). Consequently, in the nonrelativistic 
case, the absolute velocity of a particle may become arbitrarily large during 
some intermediate time interval [t',t"] C [to,t]. Hence, the largest possible 
action value is a_|_(a;, Xq) = +00. These considerations put us in the position 
to rewrite the Wiener propagator fll99p as an integral-over-actions: 



a+{x\xo) / Q \ 

p{x\xo)(x / da exp -— , (216a) 



supplemented by the normalization condition 



dxp{x\xo). (216b) 

The representation (12161) may be generalized to the relativistic case in a 
straightforward manner: One merely needs to insert the corresponding rela- 
tivistic expressions into the boundaries of the integral (I216ap . A commonly 
considered relativistic generalization of Eq. (12131) . based on the particle's 
proper time, reads ^ 



a(a;|a;o) = — / dt' 

Jto 



i\2 



1 - v{t') 



1/2 



(217) 



Analogous to the nonrelativistic case, the relativistic action (12171) assumes its 
minimum a_ for the deterministic (direct) path from xq to x, characterized 
by a constant velocity v{t') = [x — Xo)/{t — to)- One explicitly obtains 



x\xo) = - (t - tof - (x - Xo)^ 



1/2 



(218a) 



i.e., a_ is the negative Minkowski distance of the two spacetime events Xq 
and X. The maximum action value a+ = is realized for particles that move 
at light speed[l3 Hence, the transition PDF for the relativistic generalization 
of the Wiener process reads 



In general, particles must undergo reflections in order to match the spatial bound- 
ary conditions. 
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p{x\xo) = 1^ ^ < exp 



22) 



- 1 , (218b) 



if (x — xoY < (t — toY, and p{x\xo) = otherwise, with a_ determined by 
Eq. (l2T8i!l . 

The relativistic diffusion process described by Eqs. fl218p is non-Markovian, 
i.e., it does not fulfill Chapman-Kolmogoroff criterion (12001) . The functional 
form of the propagator (]218bp remains the same for higher space dimensions 
d > 1. The normalization constants 'Nd for d = 1,2,3 read 



u 

y^d = K-^Od, (219) 

where u := t — to, Od = 27r'^/^/r(c?/2) is surface area of the d-dimensional unit 
sphere, and 'N'^ can be expressed in terms of m odifi ed Bessel functions of the 



first kind !„ and modified Struve functions Lu 485 



as 



y{[ = u7r[h{x) + L^i{x)], (220a) 
7{'^ = u^^[l + {x- 1) exp(x)] , (220b) 
27r2 

^3 = — {X [hix) + Loix)] - 2L,{x)} , (220c) 
X 

with X = m/(2D). 

In contrast to the solution (12091) of the telegraph equation, the propaga- 
tor (1218bp vanishes continuously at the diffusion fronts. Figure [TT] depicts the 
PDF g{t, x) = p{t, x|0, 0) of the diffusion process (12181) for the one-dimensional 
case (i = 1 at different times t. The corresponding mean square displacement 
is plotted in Fig. [12] (dashed curve). 

It is also interesting to note that the PDF (I216p is a special case of a larger 
class of diffusion processes, defined by 



i-a+{x\xo) 

Pu,ix\xo) = y^[w]'^ daw{a), (221) 

Ja-{x\xo) 

where w(a) > is a weighting function, and J^lw] the time-dependent nor- 
malization constant. In particular, Eq. (122 ip may be viewed as a path integral 
definition in the following sense: Physically permissible paths from xo to x have 
action values (per mass) a in the range [a_, a+]. Grouping the different paths 
together according to their action values, one may assign to each such class of 
paths, denoted by C{a;x,Xo), the statistical weight w{a). The integral (12211) 
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Fig. 11. Transition PDF Q{t,x) = p{t, x\0,0) for the one-dimensional {d = 1) rela- 
tivistic diffusion process (I218P at different times t (measured in units oi D/c^). At 
time t = to = 0, the function g{t, x) reduces to a (^-function centered at xq = 0. 
In contrast to the nonrelativistic diffusion propagator, cf. Fig. [TOl the PDF (I218[) 
vanishes outside of the hght cone. 

can then be read as an integral over the equivalence classes C(a; x, xo) and their 
respective weights w{a). The nonrelativistic Wiener process corresponds to the 
specific choice w{a) = exp[— a/(2D)]; hence, it appears natural to define the 
relativistic generalization by using the same weighting function. It is, however, 
worth mentioning that a very large class of functions w{a) yields an asymp- 
totic growth of the spatial mean square displacement that is proportional to t, 
corresponding to 'ordinary' diffusion. Moreover, Eq. (122111 c ould also be used 



to describe super-diffusion or sub-diffusion processes 
asymptotic mean square displacements grow as t", a ^ 1 



lialllTl. I51l| . whose 



°^ This can be achieved, e.g., by choosing the integral boundaries as 
a_ = (x — xqY /{t — to)") a / 1 and a+ = oo, but then the variable a may not 
be interpreted as a conventional action anymore. 
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Fig. 12. Comparison of the mean square displacements divided by 2Dt, for 

the one-dimensional [d = 1) nonrelativistic Wiener process (jl99p and its relativistic 
generalization from Eq. (j218p with initial condition (I^q^xq) = (0,0). 
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6 Outlook 



We would like to conclude our discussion by summarizing problems which, in 
our opinion, deserve further consideration in the future: 

• Microscopic models. The one-dimensional binary collision model from Sec- 
tion 14.51 appears to be the simplest example for motivating relativistic 
Langevin-type equations by means of an underlying microscopic model. Fu- 
ture investigations should focus on constructing relativistic Langevin equa- 
tions from more precise particle-field interaction models; e.g., one could 
consider the motion of a classical relativistic point particle in quasi-static 
external random fields. If successful, this approach would yield more precise 
noise and friction models for relativistic systems] I In particular, this would 
contribute to clarifying under which circumstances stochastic differential 
equations may provide a feasible approximation to complex relativistic sys- 
tems. A useful starting p oint for future research in this d i recti on could be 
the w ork o f Blanco et al. 416], Johnson and Hu 518 . 519 . 52d | and Galley 
et al. 



521 



who proposed to model the interaction between quantum test 
particles and vacuum fluctuations of quantum fields by means of effective 
Langevin equations. 

Relativistic processes driven by non-Gaussian noise. The most frequeritly 
studied relativistic stoch astic diff erential equations 11, 12, 13, [l3, [13, El, 



20|, l2l|, M, M, l32|, l33|, M, M, mM are driven by Brownian motion (Wiener) 



processes which couple to the momentum coordinates. It would be inter- 
esting to also consider other driving processes (e.g., Poisson or Levy noise) 
and to co mpare w ith the results of the corresponding nonrelativistic equa- 



tions 



83, 467, 470 



On the one hand, the properties of such generalized rel- 
ativistic processes deserve to be studied from a purely mathematical point 
of view, e.g., with regard to potential modifications of normal or anomalous 
diffusion effects. On the other hand, it would be desirable to identify classes 
of physical systems that can be described by these processe s. Fo r example. 



it seems likely that the quasi-random particle acceleration |522l | in various 
astrophysical systems may be efficiently modeled by non-Gaussian driving 
processes. 

Relativistic fluctuation theorems. Fluctuation theorems are mathematical 
relations that may be used to deduce certain thermodynamic properties of 
equilibrium systems by measuring s uitable averages fro r n dif ferent realiza- 
tions of a non-equilibrium process Isii . [sii . B . Isiel . [528]. During 
the past two decades, numerous fluctuation theorems have been established 
for various closed and open nonrelativistic systems. Both from a theoretical 



In this context, technical and conceptual challenges are likely to arise when con- 
sidering the energy loss of the particle due to radiation emission 
516l.[517j. 
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and practical point of view, it would be most interesting to generalize these 
results to the rela tivist ic case. Recently, a first step in this direction was 



made by Fingerle (395 



who derived a relativistic fluctuation theorem for 

It will be 



the special relativistic Brownian motion process proposed in [17 



worthwhile to extend these investigations to other, more general relativistic 



processes (see also Cleuren et al. 24 1| ) 



Relativistic diffusion models in Minkowski space. The above problems refer 
to stochastic processes in relativistic phase space. Alternatively, one may 
focus on constructing and analyzing novel types of relativistic diffusion pro- 
cesses in Minkowski spacetime, similar to those discussed in Section O One 
particularly important issue in this context concerns the existence of rea- 
sonable path integral representations for the propagators of such processes. 
The latter question is closely related to the problem of find i ng path in tegral 



formulations of relativistic quantum propagators ^384j . l385l . 13871 . |389 



Relativistic quantum Brownian motions. Another potential generalization of 
the above ideas concerns the construction of relativistic quantum Brownian 
motion processes. Recent eff orts in th i s dir ection include the afore mentioned 
papers by Johnson and Hu [51 8l . l519l . and Galley et al. [521|, who con- 
sidered the motion of quasi-particles in fluctuating quantum fields. Similar 
concepts have also been studied within the theory of stochastic semiclas- 
sical gravity, where the gravitational field plays the role of a 'Brownian 
particle' with the vacuum fluctuations of quantum fields forming a stochas- 
tic environmen t ('b ath'); for details we refer to the recent review by Hu 
and Verd aguer 529[ |. Another promising star ting point may be the work of 
Pechukas [lQ5|, and Tsonchev and Pechukas 106|, who developed a simple 
elastic collision model of nonrelativistic quantum Brownian n iotions. In thi s 
context, we also mention the work of Breuer and Petruccione 530. 531. 532 



who proposed a Langevin equation approach to describe qu antu m state dif- 



fusion in the framework of special relativity, see also Diosi [533 



Extensio ns to general relativity. Last but not least, notwithstanding recent 



progress |395|, |494, |497|, l53J, l535|, l536|, |537|, l537|, l538|, [SSg , the generalization 



of stochastic concepts and their applications within the framework of general 
relativity offers many interesting challenges for the future. 
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T temperature 

f3 inverse thermal energy f3 := (/cbT)"^ 

U; internal energy; thermodynamic energy-momentum four- vector 

§ thermodynamic entropy 

& relative entropy 

T pressure 

V volume (number) 
a friction coefficient 
D noise coefficient 

V spatial diffusion constant 

B{s) d-dimensional standard Wiener process with time parameter s 

P probability measure of the Wiener process 

placeholder for discretization rule in stochastic integrals, e {*, o, •} 

* Ito (pre-point) interpretation of the stochastic integral 

o Stratonovich-Fisk (mid-point) interpretation of the stochastic integral 

• backward Ito (post-point) interpretation of the stochastic integral 
R set of real numbers 

V C M'' finite d-dimensional spatial volume set 

S) d-dimensional hyperplane in (1 -|- d) -dimensional Minkowski spacetime 

€{S') backward lightcone of the event 

T(^°) isochronous hyperplane {{t,x) \ t = ^"^j 

A Lebesgue measure, e.g., V = A(V) 

II, p measures 

J indicator function with values in {0, 1} 

{X) expected value of the random variable or process X 
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A Stochastic integrals and calculus 



This appendix summarizes the most commonly considered stochastic integral 
definitions and the corresponding rules of stochastic calculus. For a mor e rig - 



orous and comprehensive introduction, we refer to, e.g., Refs. |73l. |74|. l76l . |422 



We consider a Wiener process (standard Brownian motion) B{t) as defined in 
Section [2. 1.2[ i.e., the increments dB(t) := B(t + dt) — B{t) are stochastically 
independent \l% l76| and characterized by the Gaussian distribution 



¥{dB(t) G [y.y + dy]} = (27rdt)-^/'exp[-yV(2dt)]d?/. 
We are interested in defining integrals of the form 



(A.l) 



f{Y{s))QdB{s) 



(A.2) 



where /(y) is some real-valued function, Y{s) a real-valued time-dependent 
process, and signals a discretization rule discussed below. If B{s) were some 
ordinary differentiable function of s G [0,t], then the integral in Eq. (1A.2P 
would simply be given b} 
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f{Y{s)) B{s) ds. 



(A.3) 



where B = dB/ds. Unfortunately, B{s) is not well-defined for the Wiener 
process [tI, 422| . but it is possible to generalize the concept of integrati on t o 



also include the Wiener process and other stochastic processes [T^, iTg, |422 



However, in contrast to the standard Riemann-Stieltjes integral (1A.3|) . the 
integral with respect to a stochastic process may depend on the choice of 
the discretization scheme and, in particular, also require modifications of 
differential calculus. 

To illustrate these aspects for the most commonly considered stochastic in- 
tegral definitions, we will always consider the following equidistant partition 
{to, ti - ■ ■ 1 ^n} of the time interval [0, t\: 



At = ti 



tk- 



k-l 



t/N, k 



t 







0, 



t 



N 



t. 



(A.4) 



By writing Eq. (jA.Sh . it is implictly assumed that f{y), Y and B are sufficiently 
smooth functions so that this integral exists in the sense of Riemann-Stieltjes; in this 
case, the value of the integral (1A.3|) is independent of the underlying discretization 
scheme 76| . 
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A.l Ito integral 



We first summarize the properties of Ito's stocliastic integral 60|, l6l|]- Its re- 



lationsliip to other stochastic integrals is discussed in Section IA.4I 



A. 1.1 One- dimensional case 

The Ito stochastic integral of some real-valued function f{Y{t)) with respect 
to a standard Brownian motion process B{t) over the time-interval [0,t] can 
be defined by 



/ f{Y{s)) * dB{s) := hm ^ f{Y{h)) - i?(t,)] , (A.5) 

where the partition {tg, • • • ,tN} is given by (lA.4p . The peculiar, defining fea- 
ture of this integral is that, on the rhs. of Eq. (lA.Sp . the argument of the 
function / must be evaluated at the lower boundary points tk of the discrete 
intervals [tk,tk+i]; i.e., the definition of the Ito integral is non- anticipating. 
Accordingly, the Ito discretization scheme is also known as the pre-point rule. 

Now consider a stochastic process Y{t) which, for two given functions A{y) 
and C{y), is defined by 

Y{t) = F(0) + r A{Y{s)) ds + f C{Y{s)) * dB{s), (A.6) 
JO Jo 

and where the last term is interpreted as an Ito integral (lA.Sp . Stochastic 
integral equations like Eq. flA.6p are usually abbreviated by rewriting them as 
an Ito stochastic differential equation (I-SDE) 



dY{t) = A{Y) dt + C{Y)*dB{t), (A.7) 

complemented by the initial condition Y{0). From the non-anticipating def- 
inition (lA.Sp of t he Ito integral and the properties of the Wiener process it 



follows that 422 
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{C{Y)*dB{t)\Y{t)=y) = 0. (A.8) 

The Fokker-Planck equation for the PDF f(t,y) of the stochastic process 
defined by Eq. flA.Tp reads 



( • I Y(t) = y) denotes the conditional expectation with respect to the Gaussian 
measure of the Wiener process B{t). 
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where A = A{y) and C = C{y). A deterministic initial condition 1^(0) = yo 
translates into f{0,y) = 6{y — yo). 



Finally, an important peculiarity arises when one considers nonlinear trans- 
formations G of the stochastic process Y{t). More precisley, assuming that 
Y is defined by the I-SDE (1A.7|) . then the differential cha nge of the process 
Z{t) := G{Y{t)) is given by (see, e.g.. Section 4.3.2 in |42^ ) 



dZ{t) = G\Y) * dr + i G(Y f G"{Y) dt 



A{Y) G'{Y) + - G{Yf G"{Y) 

2 



dt + G{Y) G'{Y)*dB{t), 



(A.IO) 



where G'{y) = dG{y)/dy and G"{y) = d^G{y)/dy'^. Within ordinary differ- 
ential calculus, the term containing G" is absent. Equation ( lA.lOl) is usually 
referred to as Ito formula. 



A. 1.2 The n- dimensional case 

Consider the n-dimensional stochastic process Y{t) = {Y^{t), . . . , F"(t)), de- 
fined by the following n-dimensional generalization of Eq. (I A. 70 : 



dY\t) = A\Y) dt + GMY) * dB'{t), 



(A.ll) 



where i = 1, . . . ,n and r = 1, . . . , i^. In Eq. ( lA.llI) . the Wiener processes B^(t) 
represent K independent noise sources, and each term C'^r(Y)*dB^{t) symbol- 
izes an Ito integral. The Fokker-Planck equation for the PDF f{t,y^, . . . ,?/"■) 
reads 



dt dy^ 



(A.12) 



The generalized Ito-formula reads (see, e.g.. Section 4.3.2 in 422| ) 



dG[Y{t)] 
where di := d/dy''. 



A' diG + - G\G\ did.G 
2 



dt + CV 9,G*d5"(t), (A. 13) 
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A. 2 Stratonovich-Fisk integral 



Next we summarize the properties of an alternative stochastic integral defi- 



nition which was proposed by Stratonovich |67l . l68l . |69| and Fisk |65l . |66| . In 
contrast to the non-anticipating Ito integral, the Stratonovich-Fisk (SF) inte- 
gral is semi-anticipating, but satisfies the rules of ordinary stochastic calculus. 



A. 2.1 One- dimensional case 

The SF stochastic integral of some real-valued function f{Y{t)) with respect 
to a standard Brownian (Wiener) motion process B[t) over the time- interval 
[0, t] can be defined by 



/ fiYis)) o dB{s) := hm ^ - [/(F(t.+i)) + /(F(4))] x 

[5(tfc+i) - 5(4)] (A.14) 



where the partition {tg, • • • jt^} is given by flA.4p . In contrast to Ito's inte- 
gral flA.51) . the SF definition flA.14p uses the mean of the boundary values 



of / on the intervals [tk,tk+i]', i.e., the definition of the SF integral is semi- 
anticipating. This discretization scheme is also known as the mid-point rule. 

Similar to Eq. ( ]A.6[) . we may consider a stochastic process Y{t) defined by 



Y{t) = Y{0) + /* A{Y{s)) ds + f C{Y{s)) o dB{s), (A.15) 
Jo Jo 

where now the last term is interpreted as an SF integral (1A.14I) . The integral 
equation (1A.15I) can be abbreviated in terms of the equivalent SF stochastic 
differential equation (SF-SDE) 



dY{t) = A(Y) dt + C(Y)odB{t), (A.16) 



with initial condition 1^(0). From the semi-anticipating definition flA.14p o f 



the SF integral and the properties of the Wiener process it follows that 422 



(C(r) o dBit) I Y{t) =y) = l Civ) C'iy) dt, (A. 17) 

where C'{y) = dC{y)/dy. The Fokker-Planck equation for the PDF f(t,y) of 
the stochastic process (1A.16P reads 
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dt dy 



(A. 18) 



where A = A{y), and C = C{y). The deterministic initial condition Y{0) = y^ 
translates into /(O, y) = 6{y — yo). 



It can be shown (tgI, that the SF integral definition preserves the rules 
of ordinary stochastic calculus; i.e., if Y(t) is defined by the SF-SDE (lA.16p . 
then the differential chan ge of the process Z{t) := G{Y{t)) is given by (see, 
e.g.. Section 4.3.2 in |422l ]) 



dZ{t)=G'{Y)odY 

= A{Y) G'{Y) dt + C{Y) G'{Y) o dB{t), 

(A.19) 

where G'{y) = dG{y)/dy. 

However, as will be discussed in Section \AA\ for a given SF-SDE with suffi- 
cently smooth coefficient functions A and G one can always find an I-SDE, 
which yields the same Fokker-Planck equation. Hence, in order to describe a 
certain physical process, one may choose that integral definition which is most 
convenient for the problem under consideration. 



A. 2. 2 The n- dimensional case 

Consider the n-dimensional stochastic process Y{t) = (F^(t), . . . , de- 
fined by the following n-dimensional generalization of Eq. flA.16p : 



dY\t) = A'(Y) dt + G\{Y) o dB'it), (A.20) 

where i = 1, . . . ,n and r = 1, . . . , i^. In Eq. f lA.20p . the Wiener processes B^{t) 
represent K independent noise sources, and each term G^r(Y)odB^{t) symbol- 
izes an SF integral. The Fokker-Planck equation for the PDF f{t,y^, . . . ,?/") 
reads 



dt dy''- 



(A.21) 



and the transformation rules of ordinary differential calculus apply. 
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A . 3 Backward Ito integral 



We still consider a third stochastic integral definition which is also known as 
the backward Ito (BI) integral 7^, 74 1. Its relationship to the other stochastic 



integrals is discussed in Section IA.4[ 



A. 3.1 One- dimensional case 



The BI stochastic integral of some real- valued function f{Y{t)) with respect 
to B{t) over the time-interval [0,t] can be defined by 



nt N-i 

/ f{Y{s)) . dB{s) := lim ^ f{Y{tk+i)) - , (A.22) 

Jo JV-oo 



where the partition {to, • • • , ^iv} is given by (lA.4p . On the rhs. of Eq. (1A.22I) . in 



contrast to the Ito and SF integrals, the argument of the function / must be 
evaluated at the upper boundary points tk+i of the discrete intervals [tk, tfc+i]; 
i.e., the definition of this integral is anticipating. This discretization scheme is 
also known as the post-point rule. 

Similar to above, we may consider a stochastic process Y{t) which, for two 
given functions A{y) and C{y), is defined by 



Y{t) = Y{0) + f AiY{s)) As + f CiY{s)) • dB{s), (A.23) 
Jo Jo 



and where the last term is now interpreted BI integral (1X22]) . Equa- 



tion (1A.6P can be abbreviated by rewriting it as a backward Ito stochastic 
differential equation (BI-SDE) 



dY{t) = A{Y)dt + C(Y)»dB{t), (A.24) 

complemented by the deterministic initial condition 1^(0). From the antici- 
pating definition (1A.22D of the BI integral and the properties of the Wiener 
process it follows that [42 2| 



{C{Y) . dB{t) I Y{t) =y) = C{y) C\y) dt. (A.25) 

The Fokker-Planck equation for the PDF f{t,y) of the stochastic process 
defined by Eq. flA.24j) reads 
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(A.26) 



where A = A{y) and C = C{y). The deterministic initial condition F(0) = y^ 
translates into /(O, y) = 6{y — yo). 

It can be shown that, similar to the Ito integral, also the BI integral requires a 
modification of differential calculus. More precisley, assuming that Y is defined 
by the BI-SDE ([Oil) , the differential change of the process Z{t) := G{Y{t)) 
is given by 



A. 3. 2 The n- dimensional case 

Consider the n-dimensional stochastic process Y{t) = (Y-^lt), . . . , de- 
fined by the following n-dimensional generalization of Eq. flA.24p : 



where i = 1, . . . ,n and r = 1, . . . ,K. In Eq. (IA.28I) . the Wiener processes 
B'^{t) represent K independent noise sources, and each term CV(^) • dB^{t) 
symbolizes a BI integral. The Fokker-Planck equation for the associated PDF 
fit,y^, ■■■^y'^) reads 




(A.27) 



where G\y) = dG{y)/dy and G"{y) = d^G{y)/dy\ 



dY\t) = A'{Y) dt + G\{Y) • dB'Xt), 



(A.28) 



dt 



l_d_ 
2'dyi 



{G\Girf) 



(A.29) 



The generalized backward Ito-formula reads 




where di := d/dy\ 
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A. 4 Comparison of stochastic integrals 



As anticipated in the preceding sections, the three different stochastic inte- 
grals/SDEs may be transformed into each other. In particular, a given Fokker- 
Planck equation can usually be realized by any of three SDE types, upon 
choosing the coefficient functions appropriately. To illustrate this by example, 
we reconsider the n-dimensional SDEs from above, assuming identical noise 
coefficients CV but different drift coefficients Al^^^^{Y), respectively, i.e. 



dY\t) = Al{Y) dt + C\{Y) * dB'it), 
dY'{t) = Ai{Y) dt + C\{Y) o dB'-it), 
dY\t) = Al{Y) dt + C\{Y) • dB'Xt), 



(A.31a) 
(A.31b) 
(A.aic) 



where i = 1, . . . ,n and r = 1, . . . ,K. We would like to determine the drift 
coefficients such that these three different types of SDEs describe the same 
n-dimensional stochastic process Y{t) = {Y^{t), . . . , Y"'{t)) on the level of the 



Fokker-Planck equationa I. which can be compactly summarized as follows 



dt 



(A.32) 



where di := d/dy\ and A* = 0, Ao = 1/2, and A, = 1. We distinguish three 
cases. 



Eq. flA.31ap is given: Eq. ([021) implies that Eqs. flA.31bp and flA.31cD 
describe the same process like Eq. (IA.31al) . if we fix 



Al = A:--C^d,C' 



— M _ r<3 pj.r<i 



(A.33) 



Eq. (]A.31b|) is given: Eqs. (]A.31aP and (]A.31cP describe the same process 
like Eq. (IA.31bD . if we fix 



M — M r<j B-r"^ 

/I, — /1q 2 r'-']^ r- 



(A.34) 



° For most practical purposes, two Markovian stochastic processes can be consid- 
ered as physically equivalent if their PDFs are governed by the same Fokker-Planck 
equation. 
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Eq. (lA.31cP is given: Eqs. f lA.SlaP and (lA.31bp describe the same process 
like Eq. (lA.Slcp . if we fix 



K = K + C\d,C\ , a: = Al + ^ C\dfi\. (A.35) 



To summarize, by means of Eqs. (lA.SSp . (]A.34p and (1A.35P one can change 



between the different forms of stochastic integration and stochastic differ- 
ential calculus, respectively. Each SDE type has advantages and disadvan- 
tages: The Ito formalism is well suited for numerical simulations 73, 89, 420l 



and yields a vanishing noise contribution to conditional expectations of the 
form flA.81) . The Stratonovich-Fisk approach is more difficult to implement nu- 
merically, but preserves the rules of ordinary differential calculus (in contrast 
to Ito/backward Ito integration). Finally, within the backward Ito scheme, 
fluctuation dissipation relatio ns ma y take a particularly elegant form (cf. Sec- 
tion 6.2 in Ref. ^72], and Ref. [m^). 



A. 5 Numerical integration 

A detailed introduction to the numerical simulation of SDEs can be found 
in 



73l . 1891 . 1420j . A simple Monte-Carlo algorithm for numerically integrating 
Eqs. (lA.3ip follows directly from the deflnition of the stochastic integrals. The 
corresponding discretization scheme, which works sufficiently well for many 
purposes, reads 



Y\t + At) - Y\t) = KiY{t)) At + C\iYit)) AB'\t), (A.36a) 
Y\t + At) - Y\t) = Ai{Y{t)) At + 

^[C\{Y{t + At)) + C\{Y{t))] AB\t), (A.36b) 

Y\t + At) - Y\t) = Ai{Y{t)) At + C\{Y{t + At)) A5"(t). (A.36c) 

Here, the AB'^{t) are random numbers, sampled from a Gaussian normal dis- 
tribution with density 



As evident from Eqs. (lA.36p . for given functions Al and CV, the discretized 
I-SDE (]A.36aP allows for calculating the values + At) directly from the 



preceding values Y'^{t). By contrast, the discretized SF-SDEs (]A.36b|) and BI- 
SDEs (IA.36cl) are implicit equations, which must be solved for Y^{t + At). 
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The latter difficulty can be avoided by transforming a given SF/BI-SDE to 
the corresponding I-SDE by means of Eqs. (]A.33|) . (1A.34I) and (]A.35I) . 
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B Surface integrals in Minkowski spacetime 



We would like to integrate a tensorial quantity 6^'^^'''{t,x) over a fixed d 
dimensional hyperplane ^ in (1 + (i)-dimensional Minkowski space, e.g. 



T^-[^]:= I da^r"«-(t,x). (B.l) 

Relevant examples are integrals over the energy-momentum tensor 9°"^ or the 
current density as discussed in Appendix O We outline the general proce- 
dure for the case of c? = 3 space dimensions, considering Cartesian coordinates 
with metric tensor {rjap) = diag(— 1, 1, 1, 1). In this case, the surfa ce el ement 
da^ may be expressed in terms of the alternating differential form 2201 



da^ = ^e^ap^dx" A dx'^ A dx^, (B.2) 



where e^afs-y is the Levi-Cevita tenson I and 'A' denotes the antisymmetric 
product 



dx" A dx^ = -dx^ A dx", (B.3a) 
implying that 

dx° A dx° = dx^ A dx^ = . . . = 0. (B.3b) 



Inserting flB.2p and ordering differentials with the help of Eq. flB.3ap . the 



integral fIB.ip can be rewritten as 



g"-[f,] := f dx^ A dx^ A dx^ - / dx° A dx^ A dx^ O^''- + 

I dx° A dx^ A dx^ - - / dx° A dx^ A dx^ e^""- . (B.4) 



85 j£ gfj.ai3... tensor of rank n then the quantity 9"^'"[^] from Eq. (jB.ip is tensor 
of rank (n — 1). 

The total antisymmetric covariant Levi-Cevita tensor Ea/S'yS is if two or more 
indicies are equal, -|-1 for even permutations of the indices (0123) and —1 for odd 
permutations. Similar to the metric tensor, the Levi-Cevita tensor is numerically 
invariant under Lorentz transformations with determinant -|-1. For a general dis- 
cussion of the properties of Levi-Cevita tensors we refer to Section 5.5 of Sexl and 
Urbantke M. 
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With regard to thermodynamics, we are particularly interested in integrating 
over space-like or time-like surfaces Sj given in the form 



= t = g{x) = g{x\...,x'^). 



(B.5) 



Typical examples are the isochronous hyperplane T(^'') of an inertial frame 
S, defined by 



:={{t,x)\t = e =:9ix) }, 



fB.6a) 



or the backward lightcone of some spacetime event with coordinates 
(^°, ^) in S, which is given by 



:={{t,x)\t = e-\x-^\=: g{x) }. 



(B.6b) 



Given such an explicit representation of the hyperplane, we may express the 
differential dx" in Eq. flB.4p in terms of dx* by using 



dx° = dig dx'. 



(B.7) 



Inserting this expression into Eq. (IB.4P and taking into account Eqs. (IB.Sp . 
we obtain 



= / dx^ Adx' Adx' 



d'^x 



e'^-ig{x),x)-{d,g) e^'^-{g{x),x) 



(B.8) 



Hence, in the case, we may write for the covariant surface element four- vector 



(dcTa) = nQ.d'^x 



rir 



'X-dig)- 



(B.9) 



In particular, for the isochronous hyperplane T(^°) from Eq. (lB.6ap . we have 
g{x) = and dig = in E, yielding 



(B.IO) 



For comparison, when integrating over the lightcone from Eq. flB.6bp . 

one has to use 



dig 



X* — 



fB.lll 
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such that n^n" = 0, yielding exphcitly 



^^"■■■(^° - I a;- II, a;) I. (B.12) 



^ pia...(cO 



Example 1: Gas in a cubic container. We calculate S"[<^] for the energy- 
momentum tensor from Eq. (11061) . describing an isotropic, stationary, 
spatially homogeneous gas in the lab frame S. Considering a cubic vessel 
V = [-L/2, L/2Y in S, the components of 9'^^ in S read 



87 



r (e) , « = /? = 0, 

a;) = V-^ ^{x- N)\{p.v)/d, a = P = l,...,d, (B.13) 

where V = L*^. Since is diagonal and time-independent in S, we find 

gO[^] = I d^x ^°°(c/(a;), a;) = V^-^ (e) 1 d^x J(a;; V) = (e) . (B.14) 

The spatial components g*[^] can be calculated by means of a partial inte- 
gration 

5'W] = -J d'x id,g) e'\g,x) 

= -iVd)-'{p-v) I d^x id,g) J(a;;V) 
= (Vd)-^ (p-v) J d'^x g{x) di3{x;Y) 

= {Vd)-'{p-v) Jd^'xie-lx-^Dd.Jix-Y), (B.15) 
where, cf. Eq. ffTOTbl) . 



d,J{x; V) = [5(L/2 + x') - 5(L/2 - x')\ [] e(L/2 + x-'') e(L/2 



87 fj-Qj]^ ()B.13P is the energy-momentum tensor per particle and has to be 
multiplied by the particle number N to obtain the total energy and momentum of 
an A^-particle system. 



117 



is antisymmtric with respect to x*; hence 



g^[^] = -{Vd)-^ {p ■ v) J d'^x \x - II did{x; V). 



In the one-dimensional case, where V = L and d = 1, this integral can be 
solved exactly, and one finds 

g^[^] = -L-' (A^) i\L/2 + e\ - \L/2 - e\)- (B.16) 

In particular, for a lab observer located in the origin ^* = we have S'^[(S'] = 0. 
In this example, the origin may be interpreted as the photographic center-of- 
mass, which can be defined by the condition S^[«f] = in the lab frame. 

If the picture is taken from outside the volume V = [—L/2, L/2], corresponding 
to the condition L/\^^\ < 1/2, then 

9'[^] = -{pv) {e/\e\y, m?) 

i.e., if the mean momentum is naively sampled from a photo taken outside of 
the photographic center-of-mass, then even the lab observer (who is not mov- 
ing with respect to the system) will find a non-vanishing momentum average 
value always pointing away from herself. In the case of a Jiittner gas with 
{p ■ v) = d/f3, the magnitude of this photographic drift effect is proportional 
to the temperature (and particle number). 

In the d-dimensional case, performing the dx*-integration first, we find 

g^[^] = J d'/^x [{x' - e? + ... {L/2 - ef + . . - 

I d'I'x \{x^ - + ... {L/2 + Cf + ■ ■ .]'/^ (B.18) 

with the remaining {d— l)-dimensional integration (denoted by d'^/'x) ranging 
over [-L/2,L/2]'^-\ Again, S'W] vanishes if = 0. In the limit L/\^\ < 1/2, 
corresponding to an observer position far outside the volume, one may expand 
the integrands, yielding 

-{vd)-' (p-v) I d''/'x L e [{x' - + ■■■ {cf + ■ ■ 
^-(rrf)-i (p-v) {c/m) Jd'/^xL 

= -d-' {p-v) (e/lll). (B.19) 
In the one-dimensional case this reduces to Eq. (1B.17p . 
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Example 2: Arbitrary confinement. In the preceding example, the spa- 
tial density was taken to be piecewise constant with a singular drop-off at the 
boundaries (walls of the vessel). As a slightly more general example, consider 
an equihbrated gas in a smooth confinement, e.g., due to some external po- 
tential. Moreover, assume that there exists a preferred inertial frame S, where 
the gas is stationary with an isotropic momentum density, so that it can be 
described by a Lorentz-scalar one-particle density of the form 

ip{x,p) = g{x)(Pip), (B.20) 

with normalized marginal distributions g and (p, and being spherically sym- 
metric. In this case, the (kinetic) energy- momentum tensor 



e''^(t,x) ■= I -^pyP ^(x,p) (B.21) 
J 

has components 



r (e) , a = /? = 0, 

e^^^it, x) = q{x) (p^v^) , a = (3 = t = l,...,d, (B.22) 

Then, from Eq. (IB. 121) . one obtains 

g°[^] = (e) J d'x 0{x) = (e) , (B.23) 



and 



= (pV) J d'x ^-^ g{x) =: 9\^) (B.24) 
In this case, the photographic center-of-mass position in S is defined by 



SXC) = 0, t = l,...,d; (B.25) 

e.g., if g is rotationally invariant with respect to some space point z in S, 
then = z. On the other hand, assuming that the density g vanishes outside 
a finite region V G M'^ in E and that a stationary observer is located at a 
position ^ far away from V, then we may approximate 



^ {p'v') / d'^a; — ^ g{x) = ip'v') ' , . (B.26 
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C Relativistic thermodynamics 



In this appendix, we shall refer to the Jiittner gas 188l | from Section 13.2.21 in 
order to compare different formulations of relativistic thermodynamics 
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1221 . 11271 . Il50l . Il52| . This serves to illustrate how different definitions and con- 
ventions imply different Lorentz transformation formulas for, e.g., the temper- 
ature. 



C.l Reminder: nonrelativistic thermodynamics 



Traditionally, nonrelativistic thermodynamics intends to describe a many- 
particle systemsby means of a few macroscopic control parameters. Accord- 



ing to Callen Hsil, Isiol . the most 'natural candidates for thermodynamic 



variables are either conserved quantities (such as total internal energy IX in 
the rest frame, particle number N, etc.) or parameters that characterize the 
breaking of symmetries (e.g., the volume parameter V characterizes the break- 
ing of translation invariance, external magnetic fields may break iso t ropy , 



etc.). Within the axiomatic formulation of thermodynamics |48ll . l540l . |541 



one further postulates the existence of an extensive thermodynamic potential 
§(11, V, . . .) called entropy, required to be a monotonous function of U. Given 
the entropy S, intensive thermodynamic quantities (temperature T, pressure 
7, etc.) are defined by a differentiation of S with respect to the extensive 
control parameters (IX, V, . . .), yielding the first law of thermodynamics 



dS = 7-\dU + 7dV + ...). 



(C.la) 



Upon introducing a heat differential form dQ and a work contribution d^l 

by [liol-liii) 



dQ:=TdS, dA:=-7dV - 

the first law fIC.laP can be rewritten in the form 



(C.lb) 



dQ = du-dyi. 



(C.lc) 



In the nonrelativistic framework, it is rath er s traightforward to generalize 
Eqs. ( IC.ll) to also describe moving systems 150|. To illustrate this, consider 



a nonrelativistic thermodynamic system of constant total mass M, moving at 
constant mean velocity w. In this case, one can interpret the mean veloci ty w ' 
as an additional intensive variable, and define the internal energy U by jisj 



120 



£' = U + 

where £' is the total energy of the moving system; hence, 



(C.2) 



dU = d£' - w' ■ d{Mw') = d£' - w' ■ dG', 



(C.3a) 



where G' := Mw' is the total (mean) momentum of the system. Combining 
this with Eq. (IC.ip . the first law can be expressed as 



dQ = d£' - w' ■ dG' - dyi, (C.3b) 

where the second term on the rhs. may be interpreted as acceleration work. 
C.2 Relativistic case 



The relativistic g ener alization of Eqs. (IC.ip - flC.3|) is less trivial due to the 



following reasons [150 



One cannot decompose the relativistic energy £ of a thermodynamic system 
into the sum of a kinetic part and a term IX depending only on the internal 
state. 

In general, the rest mass (energy) of a thermodynamic system is not constant 
in thermodynamic processes, since, e.g., any heat transfer dQ represents a 
change of energy and, therefore, of mass. 

The transfer of energy and momentum between moving systems implies the 
transfer of mass and, hence, of momentum. 

There exist several different, reasonable ways to define heat and work in 
relativistic systems. 



When adopting a statistical approach! I towards relativistic thermodynamics 



i.e., when identifying thermodynamic variables with statistical averages, one 



can add one more source of difficulty to this list [178 



Equations (jClh and (IC.3|) represent relations between macroscopic observables. 
The microscopic, statistical justification of these equations rests on the assumption 
that thermodynamic variables can be related to we ll-de f ined expectation values with 
respect to a phase space probability distribution [2871 . l49l[ |. A statistical distribu- 
tion (ensemble) with density / provides a "good" model of thermodynamics if the 
expectation values satisfy the differential rela tion (IC.3P for a suitable entropy func- 
tional §[/]. Recently, it has been shown [49ll | that thermodynamic relations of the 
type (1C.1|1 hold for a rather wide class of distributions (not only micro canonical and 
canonical ones), provided the entropy functional is chosen appropriately. 
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• Expectation values such as the mean energy or the mean momentum of 
a many-particle system are nonlocal quantities whose definition requires 
choosing (i.e., fixing) a specific hyperplane in spacetime. The choice of the 
hyperplane may single out a preferred frame of reference. 

In particular, the last two aspects have led to considerable debate and confu- 



sion during the past century, e.g., coi icerning the transforma t ion o : 
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In the remainder of this section we briefiy summarize three diffe rent formalisms 

Ott 



of relativistic thermody namics, as proposed by a Planck [121 



127 and 



van Kampen 15Cll . 152 |. respectively. This serves to illustrate how different 
definitions and conventions yield, e.g., different transformation formulas for 
the temperature. For this purpose, we consider a confined system with fixed 
particle number that, in the lab frame S, can be described by a spatially 
homogeneous, stationary Jiittner distribution flll8p . Aiming to identify ther- 
modynamic quantities with statistical averages, we distinguish three different 
types of hyperplanes in Minkowski spacetime, as illustrated in Fig. IC.ll : 

We define the S-isochronous hyperplane corresponding to events {t,x) 

with t = ^'^=constant in the lab frame S, by 



%{e):={{t,x)\t = e, xeR''}. (C.4a) 

Similarly, we define for a frame S' moving at velocity w = {w, 0, . . . , 0) relative 
to the lab frame S, the S'-isochronous hyperplane by 

T'(^'O) := { (t', x') I t' = x' G }. (C.4b) 

We denote the backward lightcone of some spacetime event (f, having 

coordinates ^(^) = (^°, ^) in S and coordinates (^'°) = (A°^^^) in E', by 

) := {{t,x)\t = ^^ -\x-^\, xeR'^ }. (C.4c) 

We next introduce relativistic thermodynamic variables as integrals over these 
hyperplanes. 



Energy and momentum We distinguish energy-momentum mean values 
defined along the different hyperplanes] 



Although the energy-momentum tensor was introduced above with the help of 
the one-particle phase space PDF, Eqs. (IC.SP can also be used to develop more 



122 



t = e 

t' = f 



-L/2 



S 



(0,0) 




L/2 



Fig. C.l. The different liyperplanes as defined in Eq. ()C.4p . Tlie worldlines of two 
particles are labeled by "PI" and "P2"; the worldlines of the container walls corre- 
spond to vertical lines at x = —L/2 and x = L/2 (thick dashed lines), respectively. 
Assume a lab observer, resting at position x = in S, takes a photograph of the sys- 
tem at the spacetime event S' with coordinates {t, x) = {^^ , 0) in S. This photograph 
will reflect the state of the system along the lightcone ^((^). 

(C.5a) 
(C.5b) 
(C.5c) 

In the presence of an arbitrary external confinement, the energy-momentum 

cf. Eq. (11071) . In this case, each 
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nonlocal four- 



tensor is not divergence-free everywhere, 
of the three quantities from Eq. (1C.5I) represents a differentL 
vector, since they refer to different hyperplanes, respectively. In principle, 
either of the three four- vectors could be used as the rmodynamic variable. 



giving rise to different thermodynamic formalisms [178 



The lab-isochronous four- vector 'U^[^°] from Eq. (IC.SaP is most easily calcu- 
lated in the lab frame S; cf. Appendix [Bl Using Eq. (llOScI) . one finds 




/3 = 0, 



(C.6a) 



general thermodynamic theories based on the energy-moment um tensor (and other 
conserved currents) constructed, e.g., from field theories 

Similar problems arise ifone wants to construct a relativistic continuum model 
for the electron [lii. l475l.[478|. 

The three four-vectors from Eq. (jC.Sh would coincide if the energy-momentum 
tensor were divergence- free everywhere |178i] . 
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In the case of a Jiittner distribution, the one-particle expectation value (e) = 
(e)^ is given by Eqs. (]120b|) - (]120dp . In another frame S', moving with velocity 
w along the x^-axis of S, one finds the transformed energy-momentum vector 
lX'^[e°] = Aiwy^Wi^] as 




/3 = 0, 
f3 = l, 
P>1, 



(C.6b) 



where 7(w) = (1— w^) and A(w )^^, denotes the Lo rentz transformation ma- 
trix. As discussed below, Ott 127| and van Kampen 150 . 152 1 use Eqs. (IC.6P 



as starting point for their formulations of relativistic thermodynamics, but 
consider different definitions of heat and temperature, respectively. 

For comparison, the S'-isochronous energy-momentum four-vector, defined in 
Eq. (]C.5bl) by integration along the hyperplane is most conveniently 

calculated in E', yielding, by virtue of Eq. (lllSbp . 



fOl 



Nip' 



h{{^)+w^p-v) /d), (3 = 0, 

Nl-^w{{e) + {p-v)/d), (3 = 1, (C.7a) 

[o, (3>1, 

where {p ■ v) = d/(3 for the Jiittner distribution, cf. Eg. 01220. Eq uation (ICTaP 
presents the basis of the Planck-Einstein formulation 12ll . Il22( | of relativistic 
thermodynamics, cf. discussion below. Applying an inverse Lorentz transfor- 
mation A{—w) to Eq. flC.7al) gives 



■(e), (3 = 0, 

U^[e] = N{-w{p-v)/d, (3 = 1, 
0, (3>1. 



(C.7b) 



The non-vanishing component with (3 = 1 reflects the fact that the integra- 
tion was performed along the hyperplane "t'=constant" , which introduces an 
apparent asymmetry in the lab frame, cf. Fig. IC.ll and results in a spurious 
mean momentum 
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borne authors |l23, [l5l|, [iHJ] have interpreted the w^-term in Eq. ()C.7p as an 
energy co ntrib ution due to presence of the walls. This view, although objected to 
by others 178| . seems to be at least partially correct. The above derivation shows 
that the difference between U'^f^"] and U'^f^"^] can be attributed to the different 
underlying hyperplanes X(^'') and X'(^'°), respectively. However, Eqs. (jC.SaP and 
(fCl5al) would give the same result if the energy-momentum tensor were divergence- 
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In fact, th e most fr e quently discussed versions o f relativistic thermodynamics 



see, e.g.. (Til l , [lii. [Tirl. [Tisl. jiiol . CT . [Tsi . [m, 



Tii, CT. Il4ll. [l42ril43 
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MaCTlTHliii iTsa iTsi . [159 



160 



138 



165 



lee, lea. Iml . I17I Il78l . l214j . ]215l . I217I . I234I . kTsj Vadopt - imphcitly or exphcitly 



- either Eqs. (IC.ep or Eqs. (IC.7P as the starting point. From an experimental 
point of view, these equations refer to different, nonequivalent measurements 
procedures: 

In order to determine 'U'^[.^°] a moving observer (at rest S') had to reconstruct 
the velocities and/or momentum values along the hj^erplane t = whereas 
to obtain U'^ [^'^] velocities must be sampled S^-sim ult aneously along t' = 



On the other hand, as pointed out by Gamba |475 |. even for a lab observer it 



will be very difficult (if not impossible in practice) to reconstruct the velocity 
data along the hyperplane t = due to the finiteness of the speed of light. 



Isochronous vs. lightcone averages A more natural way of performing 
measurements (e.g., in astronomy) is to 'take a photograph' of a given system. 
A photograph recorded by an observer at the event reflects the state of 
the system along the lightcone ^[S'], see Fig. IC.li Let us assume that an 



idealized photograph encodes both position and velocitM I of the particles in 



the system. Then, the empirical energy and momentum averages which can be 
sampled from this data correspond to the lightcone average 'U^[^] as defined 
in Eq. (lU^cj) . 

Calculating 'U''[(f] for a stationary lab distribution of the form (p{x,p) = 
g{x)(f){p), with isotropic momentum PDF 0, gives in the lab frame (cf. Ap- 
pendix [B]) 



1X°[ 



N{e), 

d'^X f ^ 

\x — £1 



QK.X) 



(C.Sa) 
(C.8b) 



where (^°, ^) are the spacetime coordinates in S. By comparing with Eqs. (lC.ea|) 
and (lC.7bp . we see that the 0-components are identical, 



'Ol. 



(C.9) 



free everywhere |l78l ]. In the above example of a homogeneous gas, the divergence 
is non-zero only on the boundary due to the appearance of the G-function in the 
phase space density, cf. Eq. ()107p : i.e., the difference between U'^[^''] and U'^[^"^] is 
indeed related to the presence o f the boundary. It seems that Einstein was aware of 
this problem, cf. his remarks in 122l |. but at that time did not follow up this issue. 

In practice, velocities can be reconstructed from color variations due the Doppler 
shift of spectral lines. 
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i.e. when sampling energy values it does not matter whether this is done from 
a photograph or from simultaneously collected (i.e., reconstructed) data. 



The situation is different, when estimating the mean momentum. As evident 
from Eq. flC.8b|) . even for a lab observer at rest in S, the lightcone average 



depends on the position ^ of the observer J ^ I A distinguished "photographic 



center-of-mass" position in E can defined by 



U'[^] = , i = l,...,d. (C.IO) 

For example, if g is symmetric with respect to the origin of S, then = 0; 
this would correspond to a lightcone as drawn in Fig. IC.ll 

To illustrate how WIS"] depends on the observer position, assume that the 
density g vanishes outside a finite region V G in S and that a stationary 
observer is located at a position ^ far away from V. In this case, one can 
approximate |a; — ^| ~ |^| in the integrand of Eq. (IC.Sbp . yielding 



Nd-' {p ■ v) 



(C.ll) 



In particular, when considering a homogeneous Jiittner gas with position mean 
value (x) = and {p ■ v) = d/f], ther. 
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C N 



(C.12) 



i.e., an outside observer gas, who naively estimates W'lS'] from his photo- 
graphic data, could erroneously conclude that the gas is moving away from 



him (at a velocity proportional to the temperature). This effect! I should be 



taken into consideration when estimating the velocities of astrophysical ob- 
jects from photographs. 



From the conceptual point of view, it is worthwhile to note that the hyperplane 
"backward lightcone C(<^)" is a relativistically invariant object which is equally 
accessible for any inertial observer. Put differently, if a second observer, moving 
relative to the first one, takes a snapshot at the same event then her picture 



Clearly, averages in the lab frame do not depend on the specific value ^ of the 
time coordinate if the PDF is stationary in this frame. 

Note that Eqs. (|C.8bp and ()C.12p are consistent with Eq. (|C.7bp . as can be seen 
by letting > 0, = _ _ = = q ^nd -u; ^ 1 in Eq. (lafb]) . 

The effect becomes neglible if j3m(? ^ 1. 
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will reflect the same state of the system - even though the colors will be 
different due to the Doppler effect caused by the observer's relative motion 



Another advantage of lightcone averages as defined in Eq. (]C.5cP lies in the 
fact that they can be generalized to general relativity [s], 0| in a straightforward 
manner, whereas it becomes very difficult to single out a globally acceptable 
"simultaneous" hyperplane in curved spacetimeL^ In view of these benefits, 
hghtcone averages appear to be the better - if not the best - suited candidates 
if one wishes to characterize a many-particle system by means of nonlocally 
defined, macroscopic variables within a relativistic framework. 

However, since historically most authors considered either lab-simultaneously 
or observer-simultaneously defined quantities, we will restrict ourselves in the 
remainder to discussing the implications of these two choices. 



Entropy Having identified the potential candidates for the thermodynamic 
state variables 'energy' and 'momentum', one still needs to specify 'entropy'. 
In the case of a Juttner gas, one can define an entropy density four-current 



per particle by |242l . |312 



s-(t,x) = - J^p" f{t,x,p)\n[h'f{t,x,p)]. (C.13) 



The 'logarithmic' form of this entropy current is specifically adapted |491 
to the exponential form of the Jiittner distribution and/or vice versa, cf. 
Eq. (^2M>- Inserting the Jiittner function ffTTSD into Eq. flClSl) . one finds 
that the entropy current is stationary in the lab frame E and given by 



."(«..) = V'-'5(.;V)^'"(^^^/"') " = (C.14) 

0, a = 1, . . . ,d. 



Hence, the current (IC.14p satisfies the conservation law 



= 0. (C.15) 

The thermodynamic entropy § is obtained by integrating Eq. flC.14p over some 
space-like or light-like hyperplane S), yielding the Lorentz invariant quantity 

§[^] ■=N f daa = N f da'^ s"'{t\x') =: S'[Sj]. (C.16) 



^^In the nonrelativistic limit case c — > oo, the lightcone "opens up" so that photo- 
graphic measurements become isochronous in any frame in this limit. 
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The conservation law (IC.ISP implies that the integral (1C.16|) is the same for 
the hyperplanes T(^°), T'(^'°) and i.e., 



(C.17) 



The integral (]C.16P is most conveniently calculated along Sj = T[^'^] in S, 
yielding 



§ = N Jd'^x x) = N \n{VZd/h'^) + (3 N (e) . 



(C.18) 



Heat and temperature By means of the preceding considerations, we 
can now summarize and cor apare the most cornmonly d i scus s ed versions of 
rela t ivist i c thermodynamics 

lii 
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137 . 



160 



138 . 



165 



139, 



166 



140 



169 



121 



141 



171 



122 



142 



172 



127 



143 



178 



128 



144 



214 



129 



145 



215 



1301. 



14f 



211 



131 



147 



234 



132, 



148 



475 



133, 



149 



134 
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. As men- 
tioned above, these versions are based on simultaneously defined averages, 
respectively. In spite of this common feature, differences exist regarding 

• the choice of the underlying hyperplane and/or 

• the adopted definition of heat, 

leading, e.g., to different temperature transformation laws. Taking into ac- 
coun t the hist orical order, we begin by recalling the Planck- Einstein formula- 
tion 



1211 . |122| | of relativistic thermodynamics. 



Guided by Eqs. (lC.3p . Planck and Einstein opted in 1907/08 for the following 
definition of heat {cf. Eq. (23) in [l2^ } 



dQ'(w') := d£' - w'dS' + 7'dV' =: 7'dS' 



(C.19a) 



where w' = —w is the constant velocity I of the thermodynamic system in 
the frame S', anc 



For simplicity, we again consider a thermodynamic system that moves with ve- 
locity w' along the x'^-axis of S', so that w' = {w' ,0, . . . , 0) = {—w, 0, . . . , 0) = —w 
with w denoting the velocity of S' relative to the restframe S of the system. 

The scalar transformation (IC.lQcl) of the pressure T is implied by the transfor- 
mation laws of force and area 



121 



im 



^'^''After reinstating constants c in Eq. (lC.19fjl . it becomes evident that S'{w' 
w'M in the nonrelativistic limit case c 



oo. 
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V'{w 
T'{w 

8.'{w 
9'{w 



s, 

7 (£ + w'^ TV 



(C.19b) 
(C.19c) 
(C.19d) 

(C.19e) 

(C.19f) 



with 7 = 'y{w') = 7(w) denoting the Lorentz-f actor. The choice flC.19bp - 
(]C.19fp corresponds to defining thermodynamic energy and momentum observer- 
simultaneously as in Eqs. (lC.7aP and identifying the pressure as 



9V = N{p-v) /d, 
so that for a Jiittner gas we have 



(C.19g) 



-PV = N/3-\ 

Furthermore, substituting Eqs. (lC.19bp - (lC.19f[) into Eq. (lC.19aP gives 



dQ'(w') = 7-1 (d£ + 7dV)= 7-^ dQ = 7-1 TdS, 



(C.19h) 



and comparing Eq. (]C.19hD w ith (lC.19ap leads to the temperature transfor- 
mation formula of Planck [121,] I 



r = 7-17= (l-«;'2)i/2 7, 



(C.19i) 



stating that a moving body appears cool er (for the Jiittner gas 7 = l3~^). 
This formalis m was criticized by O tt 127l | in 1963 and later also by van Kam- 
pen [l50 . 152| and Landsberg 139 . 141 | - mostly because the quantities (£', S') 
and (£, S) = (£'(0), S'(0)) are not related by a Lorentz transformation. As dis- 
cussed above, this drawback can be traced back to the fact that £' and S', 
if take n at different values of w', correspond to different hyperplanes, respec- 



tively [178 



To overcome this deficiency, van Kampen [I50l . ll52l | proposed to define, instead 
of Eq. (lC.19ap . a thermal energy-momentum transfer four- vector by 



dQ^ 



d£"-d^^ 



(C.20a) 



seems that, in the later stages of his life, Einstein changed his opinion about 
the correct transformation laws of thermodynainic quantities, favo ring formulas 



which were later independently deriv ed b y Ott 1271] and Arzelies [131 



ipond] 



the corres pond ing discussion by Liu 18ll . Il82l |. Schroder and Treder 
Requardt 



Ip; cf. 
175( 1 . and 
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where (d/l") := ( ^PdV^ 0) in the lab frame E, and £" is chosen to be the 
lab-simultaneous energy-momentum vector from Eqs. (lC.6p . i.e., 



(C.20b) 



The particular choice (lC.20bp singles out the lab-isochronous hyperplane 
and in a moving frame one then has 



dQ'" = d£'"-dyi'", 

(de'°) = (7£°,7^«'£°), 



(C.20c) 



with unprimed quantities referring to the lab frame E, defined hj w' = 0. 
While essentially agreeing on Eqs. (lC.20ap - (lC.20cp and also on the scalar 
nature of entropy §' = S, van Kampen and Ott differ in how to defi ne t he 
temperature of a moving system, i.e., how to formulate the (first partu^ of 
the) second law. More precisely, Ott opted for the definition 



rdS' :=dQ'° =7(d£° + Tdl^) =7dQ° =7a-dS, 



yielding the modified temperature transformation law 
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(C.21a) 



r = 77= (1-u;'2)-i/2T; 



(C.21b) 



i.e., according to Ott's definition of heat and temperature a moving body ap- 
pears hotter. 



However, van Kampen was able to demonstrate by example 15d | that the 
Eqs. (IC.2ip are not well-suited if one wants to describe heat and energy- 
momentum exchange between systems that move at different velocities (hetero- 
tachic processes). To obtain a more convenient description, he propos ed to 
characterize the heat transfer by means of the scalar quantity 150l . 1521 



dQ' := 



-WadQ" = dQ = dQ° 



(C.22a) 



with {w'^) = {—'-f,'~fw') denoting the velocity four-vector of the system in E', 
and reducing to (wa) = (—1,0) in the lab frame E. He then defined the 
temperature by 



^o^The "second part" 
^'^^ee also Eddington 



542 1 of the second law states that dS > in a closed system. 



12|1 and Arzelies [13l[ |. 
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T'dS' :=dQ' =dQ = TdS, 



(C.22b) 



so that, in view of S' = §, van Kampen's temperature is a scalar 



(C.22c) 



i.e., according to this definition a moving body neither hotter nor colder. 
Adopting this temperature definition, one can define an inverse temperature 
four-vector by 



PL 



r-'w'=7-'w. 



1 



(C.23a) 



which allows us to rewrite the second law in the compact covariant form 



dS' = -/3;dQ' 

Thus , the the rmodynamic formalisms proposed by Ott 



127 



150l . Il52l | are based on the same lab-isochronous hyperplane 



(C.23b) 

and v an K am- 
TtsI but 



pen 

both formulations differ only with regard to their res pective te mperature defi- 
nition. By contrast, the Einstein-Planck formulation 12 ll Il22l | is based on an 
observer-dependent hyperplane. 

To summarize, the above discussion shows that whether a moving body ap- 
pears hotter or not depends on how one defines thermodynamic quantities. 
Of course, before comparing the results of experimental observations with 
theoretical predictions one should make sure that the applied measurement 
procedures are consistent with the definitions employed in the theory (choice 
of the hyperplane, definitions of heat and work, etc.). 
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